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The papers presented here describe a new roinput 
erized system developed to project forest growi h and 
mortality with or without tree rutting for any nn\ of 
species or stand structure encountered in tin 1 field 

Growth and mortality coefficients have heendevel 
oped for the forests of the Lake States Region, and the 
system has boon used to update forest inventories, 
identify silviculturul treatment opportunities, and 
project long-run timber output and growth response. 
The system is designed to he used JIM the core of a 
larger Forest Resources Evaluation Program (KRKP) 
being developed by the North Central Station, hut H. 
lias potential for much wider use. 

Work to enhance and refine the system's capabili- 
ties is continuing, and further progress has been 
made since these papers worts prepared lor publica- 
tion. Interim modules to aecount for regeneration 
and for tree cutting to moot timber product removal 
objectives have been added. Although the rapnbili 
ties of the projection system are changing rapidly we 
thought it appropriate to publish these papers to give 
a snapshot picture of its development as of September 
1977. 

AKNK K.KEMP 

Assistant Director 

North Central Forest Experiment Station 
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The Forest and Rangeland Renewable Re- 
sources Planning Act of 1974 requires the Secre- 
tary of Agriculture to prepare a Renewable 
Resource Assessment for the Nation every 10 
years, and a long-range Renewable Resource Pro- 
gram for the Forest Service every 5 years. The first 
assessment and program were published in 1976 
(USDA Forest Service 1976a, 1976b). The second 
assessment is scheduled to be made during 1979, 
with a second program document based on this 
assessment due in 1980. These resource assess- 
ments provide the basis for national programs, but 
also can serve as guides for framing regional re- 
source policies and programs. 

State land management and planning agencies 
require resource assessments at the state and local 
level for long-range resource planning and devel- 
opment. River basin and other special-area plan- 
ners need resource inventory information for 
political subdivisions updated to a common base 
year. 

Large individual private and public forest prop- 
erties and management units likewise need re- 
source evaluation systems designed to project 
prospective supply and demand from their prop- 
erties, and to identify and evaluate potential man- 
agement opportunities. 



NORTH CENTRAL'S 
RESPONSIBILITY 

The North Central Forest Experiment Station of 
the Forest Service, U.S. Department of Agricul- 
ture, has the responsibility for conducting state- 
wide forest resource inventories and regional 
resource assessments within the North Central 
Region of the United States. To meet these and 



other needs expressed above, the Station must 
have the capability to acquire, process, and evalu- 
ate information about the current and future ex- 
tent, condition, and use of the renewable resources 
within the Region. An evaluation system must be 
designed specifically to fit the existing data base 
and data analysis capabilities already existing for 
this Region. This system must be usable for re- 
gional, State, and local assessments as the need 
arises, and be compatible with any national as- 
sessment system developed, 



WHAT SHOULD A 
RESOURCES EVALUATION 

SYSTEM DO? 

Obviously, we need an evaluation system that 
can, using basic forest resource inventory and 
other data available for a resource area, project 
future forest resource supply and demand for the 
area, and evaluate potential alternative resource 
management programs designed to change supply 
or demand. Such a system or series of systems 
should be able to evaluate timber, wildlife, recrea- 
tion, water, forage, and special uses of the forest 
resource and their interactions from an ecologic- 
environmental, economic, and socio-cultural point 
of view. There is little doubt that such a resource 
evaluation system must be computer-based if it is 
to handle the large amount of information and 
analytical procedures necessary for regional 
assessment, 

The dimensions of this job are enormous. They 
include economic, ecologic-environmental, and so- 
cio-cultural evaluations of the following resource 
use systems; timber, wildlife, recreation, water, 
forage, and special uses (fig. 1). For each resource 




Figure l.-Potential dimensions of a forest resources evaluation program. 



use system we must obtain and organize the physi- 
cal and biological information needed for evalua- 
tion, and develop methods for doing this within 
reasonable cost and time constraints. We also 
must develop and apply methods for evaluating 
economic and social costs and benefits of defined 
alternatives. Substantial work is needed to deter- 
mine the nature and extent of interrelations 
among the various resource use systems so that 
these can be considered in evaluating opportuni- 
ties and making regional assessments. These 
analyses should be made for the different resource 
components, resource owners, owner organization 
levels, and resource users (both direct and indi- 
rect), and for different spatial and temporal ex- 
tents. They also should recognize the degree of 
uncertainty regarding our knowledge. 

Because the existing forest resource is usually 
in a condition far from ideal, such an evaluation 
system must be capable of dealing with the re- 
source situation as it exists on the ground. It must 
be able to handle the entire range of sites and 
stand conditions, all levels of tree stocking, cull 
and sound trees, and almost every conceivable 
mixture of tree species and sizes. It must be able to 
produce both short- and long-run projections of 
resource conditions for a wide range of potential 
management treatments within an acceptable lev- 
el of accuracy. Further, it should be able to handle 
;he other components of forest ecosystems besides 
,rees in order to cover all forest uses. 



EXISTING SYSTEMS 
ARE INADEQUATE 

A forest growth projection system capable of 
handling the complexities of the real world in 
enough detail to evaluate alternative stand treat- 
ment programs has not been available in the past. 
Individual growth and yield models have been de- 
veloped, but mainly for pure stands of single spe- 
cies or for a few types under limited conditions. 
These, however, have proven to be inadequate in 
dealing with the mixed ages, sizes, and species, the 
non-normal stocking, or the high proportion of cull 
trees that often typify existing stands. Because 
each model was developed independently, using 
different kinds of inputs and producing different 
kinds of output, attempts to incorporate them into 
a comprehensive analytical system have not been 
successful. Also, since no growth and yield models 
exist for many forest types and stand conditions, 
any system developed from existing models would 
be, at best, incomplete. 

It is doubtful that a comprehensive system capa- 
ble of conducting integrated multi-resource, 
multi-use, multi-objective evaluations simultane- 
ously will be successfully developed in the near 
future. Rather than wait until a complete system 
is developed, the North Central Station has de- 
cided to develop parts of the system now. Since 
much work already has been done on the timber 
resource, developing a computerized system to 



evaluate the timber resource in economic terms 
was chosen as a logical first task. Later, other 
resource components will be linked to timber. 



As a First step in meeting the need for a regional 
forest resource assessment system, we are devel- 
oping a computerized timber resource projection 
and evaluation system that will update and project 
timber inventories and product yields, determine 
the extent of potential stand treatment opportuni- 
ties, and evaluate from an economic standpoint 
potential programs to increase future timber sup- 
plies. Figure 2 shows the major components of this 
system that require development. Each of these 
components in turn includes subcomponents that 
will require considerable research. 



Basic to this timber resource evaluation system 
is a subsystem to project future stand development 



in response to proposed management activities. To 
be useful for regional analyses, such a forest 
growth projection system must handle all forest 
conditions encountered in forest inventories. In 
short, it must be a generalized forest growth pro- 
jection system. 

This growth projection system must project fu- 
ture forest growth on forest inventory plots so they 
can be updated to any given base year. To do this 
adequately it must account for growth on all trees 
inventoried, mortality, recruitment or ingrowth of 
new trees, and removals of trees through cutting. 
It must cover the entire range of situ and stand 
conditions encountered. 

This growth projection system also must project 
the response of individual stands or classes of 
stands to prescribed treatments, such as precom- 
mercial thinning, partial cutting, and cull-true 
removal, so treatment opportunities can bo evalu- 
ated. When fully developed, it will also be used to 
estimate potential yields from fully-stocked 
stands to establish the potential productivity of 
specified sites. Such a system is basic to many 
types of timber resource evaluations. 
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Figure 2,Major components of a timber resource evaluation system for econom- 
ic analysis. 



NORTH CENTRAL'S 
APPROACH 

In recognizing the need for a new approach to 
this problem, researchers at the North Central 
Forest Experiment Station have developed a com- 
puter-based, generalized forest growth projection 
system that can be used as the basis for a timber 
resource evaluation system. The details of this 
system are summarized in the papers that follow. 
Although the system is still under development, it 
already can project short-run growth responses in 
enough detail to be useful in evaluating treatment 



alternatives in many stands. With further devel- 
opment it will provide the basis for a timber re- 
source evaluation system for regional forest 
resource assessments. 



USDA Forest Service 1976a. RPA: The nation's 
renewable resources an assessment, 1975. 345 
p. Washington, B.C. 

USDA Forest Service 1976b. RPA: A recom- 
mended Renewable Resources Program. 648 p. 
plus Appendices. Washington, D.C. 



Rolfe A. Leary, Principal Mensurationist 



A major new direction is surfacing in the devel- 
opment of forest growth models. The trend is away 
from single-species models toward models suffi- 
ciently general to handle any pure or mixed stand 
of any structure and density on any site. This pa- 
per explains the design features of a generalized 
forest growth projection system that can be used 
for short- to medium-range projections of existing 
stands. The details are found in other chapters of 
this report. 

Development of mathematical equation models 
of forest growth has been an active area of mensu- 
ration research for the past 40 years. One of the 
early studies reported was that of MacKinney et al. 
(1937). In the early work, goodness of fit of model 
to data was the primary discriminating variable in 
model selection. Grosenbaugh (1958) developed 
this to a highly detailed state. Little consideration 
was given to generality or theorification potential, 
Most models applied only to pure stands of a single 
species with relatively simple stand structures. 
The last 5 to 10 years have seen many efforts to 
develop growth models general enough to handle a 
variety of conditions (Pries 1974). The models for 
some species are quite detailed, giving a high reso- 
lution level, but other models, especially for eco- 
nomically unimportant species and mixed stands, 
are still at low resolution levels. Several approach- 
es (e.g., Stage 1973, Ek and Monserud 1974, Hegyi 
1975, and probably others 1 ) have been developed 
and are being improved and extended to cover 
additional forest conditions. 



I had two fundamental objectives in designing 
the model part of this growth projection system. 

l Arney, J. S. 1977. Personal communication. 
Centralia, Washington, 



The first was to capture, to the degree possible, the 
"essence" of the pi'ocess of stand growth and devel- 
opment. Thus, I sought a structural similarity be- 
tween the model and the forest stands. The second 
objective was to accomplish the first objective as 
efficiently as possible. 

From the standpoint of modeling, a forest stand 
can be viewed as a set of interactive, differential 
energy transformers and accumulators {Leary 
1970b) . Thus, to model forest stands, one must deal 
effectively with these five key concepts: (1) set, (2) 
interactive, (3) differential, (4) transformer, and 
(5) accumulator. Nearly all models address these 
concepts, although they may be given different 
names or not be identified explicitly. 

Prom "set" we know we are dealing with an 
aggregate of individuals; from "interactive" we 
know the individuals in the aggregate are inter- 
acting; from "differential" we know that not every 
individual in the set is the same; from "transform- 
er" and "accumulator" we know that the individu- 
als are both factory and warehouse. If we can 
develop a mathematics that will address each of 
the concepts and synthesize them into a whole, we 
may 'have some structural similarity between 
model and nature. 

My preference was to look first to "transformer" 
and "accumulator" when shopping for a mathe- 
matics to use. This nature-based constraint on our 
shopping list directs us to differential equations. 
For the sake of simplicity we go to first-order ordi- 
nary differential or difference equations. Here we 
see the symbolic expression of this ideal; 
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Verbally these symbols are saying: "that which 
results from growth is itself typically capable of 
growing". It is important to understand that Y is 
the dependent variable in this equation and t the 
independent. Persons accustomed to using regres- 
sion methods sometimes have difficulty, especially 
with the differential equation form, because the 
dependent variable, Y, is on both sides of the 
equality. 

Once having made the critically important 
choice of a mathematics to use, we can proceed to 
examine the possibilities for dealing with "set", 
"interactive", and "differential". The concepts of 
set and differential are complementary in the fol- 
lowing way: from set we easily consider subsets, 
but not just any subsets. They should be structured 
in a way so the individuals within a subset are 
similar and those between subsets are different. 
Subsetting or stratifying an aggregate is a natural 
way to formalize differentia. Thus, the union of 
concepts "set-subset" and "differential" suggest 
that we may need as many Y's as we have subsets. 
We turn, of course, from scalars to vectors to ac- 
commodate this need. The result, in terms of a 
symbolic expression of our ideal, is now only 
slightly altered: 
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These equations must be solved together be- 
cause they are coupled or simultaneous equations. 
We see that the mathematics of differential/differ- 
ence equations is so well-suited to our problem 
that in the handling of "bet-subset" and "differ- 
ential" we have also dealt with the concept of 
"interactive". 

The practicing forester often considers "stand" 
through more operational concepts: (1) structure, 
(2) density, (3) species composition, (4) site, and (6) 
stand history. How does the mathematics we have 
selected allow us to account for these important 
aspects of "stand"? The later sections of the paper 
go into the details of this question, so here I will 
only indicate their general occurrence in the above 
equations. 
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Thus, we see that density, site, and stand history 
are dealt with by making them arguments in the 
right-hand side of a difference equation. Stand 
structure and species composition are accommo- 
dated by changing the number of simultaneous 
equations. 

An important operational advantage made pos- 
sible by our choice of the mathematics of difference 
equations is the ease with which resolution level 
may be handled. The next few sections of the paper 
are organized around the concept of resolution 
level. They reflect the current design state, and 
are in part chronological because a conscious 
method used was one suggested by Polya (1945): 
take an insoluble (but apparently simple) problem 
and embed it into a larger problem so as to create a 
class of problems. The reasoning is that the solu- 
tion to the larger claas of problems may be more 
obvious than the solution to the one simple 
problem. 

For example, the potential choices of a model for 
a single-species stand with unimodal diameter dis- 
tribution are numerous. Generalization to handle 
mixed stands with varying stand structures and 
densities is possible with some and not with oth- 
ers. This eliminates a large fraction of the can- 
didate models. Other forms are eliminated by 
applying the reverse process, and assessing the 
ability to reduce or collapse the complicated model 
to a simple form for simple stands and leas precise 
estimates. The result of this process is a tremen- 
dous reduction in the number of feasible models 
that can make up the core of a generalized pro- 
jection system capable of variable resolution 
levels. 



HOW THE PROJECTION 

SYSTEM WORKS 

Like any system, the workings of the projection 
system may be broken up into three parts: input, 
output, and transformation or internal workings. 



Part 
input 



output 



transformation 



Form 

a list of interacting trees show- 
ing species, d.b.h., and crown 
ratio 2 for each tree. 
a list of trees showing species, 
new d.b.h. and new status if tree 
is a mortality tree, 
a stand component change 
equation, an allocation rule, 
and a mortality function. 



Input and Output 

The input list must identify each tree's species 
and d.b.h. and hopefully its crown ratio. Although 
many temporary inventory procedures do not re- 
cord trees in this detail, it is possible to fabricate a 
tree list from a stand table that would, if summa- 
rized, produce the stand table. Crown ratio is col- 
lected by the national Forest Survey, but by few 
others. This important, although usually missing, 
input variable has been successfully approxi- 
mated through a simple relation between mean 
crown ratio and 10-year mean basal area density 
(Holdaway et al. 1979). We use crown ratio as an 
indicator of stand history. It helps to tell us 
whether a stand with 80 square feet of basal area 
grew up from 40 square feet or was cut back from 
160 square feet. 

Output from the projection system is a new set of 
d.b.h. values for the trees which survived the pro- 
jection interval and a status change from live to 
dead for all trees projected to die. 

Although the input tree list is aggregated to 
form initial conditions for the system, the tree list 
is not discarded. This would be throwing away 
information. Instead, it is carried along and up- 
dated every year or periodically. This greatly facil- 
itates a flexible output of information. The job of 
summarizing the output tree list ia thereby sepa- 
rated, as it should be, from the job of summarizing 
the input tree list and growing the stand. Output 
tree list diameter class breakdowns should not be 
forced on the projection system. They are normally 
based on technology and economics, and are of no 
value in projection. 



^Crown ratio is the ratio of full live crown length 
to total tree height. 



Transforming Input to Output 

The mathematical model for transforming the 
input tree list to an updated output tree list has 
four parts: potential function, modifier function, 
allocation rule, and mortality function. 

The first two_parts fittogether as shown below: 

[mimbefl fraction of lliy 
of trees [ ' potential growth 
actually 
occurring 



In symbolic equation form the relationship is: 
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(1) 



The potential function is designed to estimate 
how rapidly the mean tree would be growing in 
d.b.h. if it were not interacting with any other 
trees (Hahn and Leary 1979). The potential 
growth is multiplied times the number of trees, 
giving a potential change in sum of diameters for 
NT trees. The modifier function reduces the poten- 
tial change to what has been observed from perma- 
nent growth plots (Leary and Holdaway 1979). In 
somewhat more detail, and with a slightly 
changed notation, equation (1) is: 
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where: 

fi is the potential function, 
j^_ is the modifier function, 
D is mean d.b.h. 
CR is mean crown ratio, 
SI is site index, 

Y is sum of tree diameters, and 
NT is number of trees. 

The allocation rule disaggregates the projected 
growth put on by NT trees, AY/At, to the individu- 
al trees on the input tree list. The simple rule: 



[percent of AY/dt 

that should be given 
\Jx> tree i 



[percent of growing 
= f stock represented 
[by tree i 



(3) 



has proven to be a useful point of beginning for this 
task (Leary, et al. 1979). 

The mortality rule takes the allocated d.b.h. 
growth for each tree and, using an appropriate 
function, computes a probability of death for a tree 



of that species growing at the allocated rate (Buch- 
man 1979) . Equation (4) shows this basic relation: 



Probability of 1 f/allocatetl d.b h \ 1 

death for i lh =f I growth for previous! , (species) -(4) 
tree of species jj [Vperiod / J 



The computed probability is compared with a 
pseudo-random number drawn from a uniform dis- 
tribution on the interval [0,1]. If the random num- 
ber is less than the computed probability, the tree 
is called a mortality tree, and its status changed to 
reflect this. It does not enter further in the pro- 
jection computations. 

Once the mortality has been specified and the 
tree list otherwise updated, it is again summarized 
for input, and equations (2), (3), and (4) are exe- 
cuted to project another year's growth. The entire 
process is repeated over and over again as many 
times as the user desires. 

Throughout this discussion it has been implicit 
that all trees in the stand would be grouped togeth- 
er and projected with a single equation. Thus, at 
this level of resolution the model is at the stand 
model end of the stand model-tree model contin- 
uum of growth models. 

In moving away from the stand model we go to 
stand component models, wherein groups of trees 
in the stand are projected with different equations. 
The key to this flexibility lies in the use of simulta- 
leous equation models and in the particular form 
>f the modifier function. The potential function, 
sllocation rule, and mortality functions are little 
affected as resolution level is changed. 



LOW RESOLUTION MODEL 
STRUCTURE 

When the projection system is operating at a low 
resolution level, the input tree list is summarized 
into two or more groups called components. For 
bhis discussion we assume that tree species is the 
Dreakdown criterion. If one has, for example, a 
bwo-species mixture, the input tree list is summa- 
'ized into two components. Each component is 
arojected by one of a system of two simultaneous 
equations as follows: 



(5) 
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AYl/At is disaggregated to trees of the first spe- 
cies using a species -specific allocation rule, and 
AY2/At to trees of the second species. Mortality 
probabilities are computed with species-specific 
coefficients in the mortality equation. An updated 
tree list is then ready to be summarized again and 
projected for another period. 

Let us go back and take a closer look at the 
modifier function, since it, of all parts of the model, 
is changed the most. The functional form of the 
modifier used to date has been: 



[l_ e /3/f(Y,NT)]. 



(6) 



The form in (6) is used when the model is for a 
single species, hence has a single equation the 
lowest resolution level possible. Here jS represents 
how the species interacts with itself in a pure 
stand. When going to a mixed stand, /3 is not 
changed, but the denominator of the exponent is, 
as follows: 
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represent how species one and two, re- 
spectively, affect themselves when grow- 
ing in pure stands, 

f b f 2 are functions of sum of diameters (Yi) and 
number of trees (NT ( ) of the species being 
projected, and 

(Hi are the coefficients that tell how species 
two affects species one (first equation) and 
species one affects species two (second 
equation). 

If 0)!=!, this indicates that species two affects 
species one the same as species one affects itself. A 



similar interpretation may be given to co 2 =l. The 
specific forms of the functions fi and f z in equation 
(7) are discussed in detail elsewhere (Leary and 
Holdaway 1979). 



INTERMEDIATE STRUCTURE 

Intermediate level models differentiate not only 
by species but also by size within species. Thus, the 
input tree list is partitioned by size class within 
species. How do we incorporate size of tree into the 
model? 

There are three separate but related aspects to 
'this question. One has to do with the mechanism 
by which the model handles size in addition to 
species. The second aspect deals with the class 
movement problem the movement of trees from 
one size class to another. The third deals with 
choosing class bounds for aggregating the input 
tree list. 

The first aspect is easily handled because the 
model form is such that equations may be added to 
the system to project size classes within species. 
When this is done the denominator of the exponent 
in the modifier function is more complex, but the 
potential function, allocation function, and mor- 
tality function remain essentially unchanged. The 
modifier function for the first equation would take 
on a form such as: 



1-e 
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The denominator of the exponent of e thus has 
the form of a linear combination of interaction 
coefficients and functional relations between the 
component being projected and all other compo- 
nents. This pattern of relations (8) is repeated for 
each size class within each species, that is, for each 
equation in the system. 

The second aspect, how to move trees from one 
class to another, is dealt with by keeping all trees 
in the class to which they were initially assigned. 
Since the trees in a class are growing, the attrib- 
utes of the class change over time, rather than the 
trees changing classes. Two substantial benefits 
result from this approach: (1) Because trees do not 



change classes, summaries of permanent growth 
plot data do not have discontinuous trends caused 
by trees growing into and/or out of a class. The 
summaries, such as basal area in a class at several 
points in time, show fairly continuous changes. 
This allows modern calibration methods to be used 
to estimate the /3 and to values (Leary 1970a, 
1970b, 1972, 1975; Leary and Skog 1972). (2) Com- 
putational efficiencies are gained in the projection 
process itself since it is not necessary to check the 
correctness of class assignment as the trees' diam- 
eters are increased. An added benefit is that the 
user has the choice of not applying the allocation 
rule-mortality submodels every year, rather only 
when output is desired. 

The third aspect is treated by partitioning the 
list of trees of each species into relative size group- 
ings rather than into size groupings set on a cardi- 
nal scale of numbers. If, for example, we use two 
size equations for each of two species, the d.b.h. 
class limit might be such that a predetermined 
percent of the trees are larger than the limit and 
the rest smaller. By using this relativistic ap- 
proach when calibrating the modifier function, 
each equation of the simultaneous system is as- 
sured of having trees present in the class it repre- 
sents. In using the model for projection, we avoid 
the computational inefficiencies incurred when 
many equations (size classes) have no trees in 
them. Of course, mortality or partial cutting could 
cause some of the relativistic classes to have few, if 
any, trees, but this can always be rectified by real- 
locating the trees to the classes. 

The increased precision that occures from add- 
ing more size equations is not known at this time. 
It is anticipated that three equations per species 
will meet many needs. 

The allocation rule works as before except that 
the percentages are based on each component-that 
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given to tree i 



percent of growing 

stock in component 1,2 
represented by tree i 



(9) 



Again, the workings of the mortality submodel are 
unaffected by resolution level, It functions the 
same as for the simplest model. 



HIGHEST RESOLUTION 
MODEL STRUCTURE 

Let's skip now to the highest resolution level 
that this model can produce. It is attained when 
the number of simultaneous equations is in- 
creased until each tree is in a class by itself, hence 
projected by a separate equation. In this way it is 
clear that there is a continuum of growth models 
from stand (all trees in one class) to tree (one tree 
in each class). 

In going to the highest resolution level the basic 
model structure is the same: 



At 



1-B 



[1] 



(Y t ) + 



, Y Jf x k - 



.(10) 



Here the denominator of the exponent of e has the 
same formulation as many individual tree growth 
models. Y A is now the diameter of the subject tree 
with spatial coordinates x A) y^. The neighborhood 
trees are of diameter Y with spatial coordinates 
xj,y ; . The function f 2 characterizes the competition 
factor common to these formulations. 

The potential function now predicts the maxi- 
mum d.b.h. growth of a single tree with diameter 
Y A , crown ratio CR A , on site, SI. The allocation rule 
is no longer needed because in equation (10) tree i 
is 100 percent of change as well as growing stock. It 
can be retained, however, by making the func- 
tional relation in (9) the identity relation. The 
mortality function and relation operate normally, 
giving an estimate of the probability of death 
based on AY fc /At rather than the disaggregated 
growth given by the allocation rule. 



EFFICIENCIES OF THE 
DESIGN 

Efficiency of the design may be viewed from a 
number of perspectives. Consider the following: 

(a) Efficiencies in modeling are permitted by 
breaking the overall problem into smaller, 
semi-independent parts. 

When the overall model is decomposed into the 
potential, mortality, and modifier functions, and 
an allocation rule, efficiencies result because all 
but one can be separately completed prior to a 



synthesis of the parts. Team efforts need this par- 
allelism. The relation of the jobs to each other and 
to the data base is shown in figure 1 for the present, 
nearly parallel effort, and for a hypothetical serial 
effort. Figure la shows that only the modifier func- 
tion development is dependent on completion of 
another part before it can be totally completed. 
This allows a rapid updating of coefficients in all 
models once a new or updated calibration data 
base is available. It also minimizes the amount of 
work needed to recover from errors. Figure Ib 
shows a serial or cumulative type of model devel- 
opment. This approach has two major disadvan- 
tages: submodel 4 cannot be completed until sub- 
model 3, and so on, and an error in submodel 1 
requires changes in submodels 2, 3, and 4 as well 
as in itself. Both of these disadvantages relate to 
the speed with which work can be accomplished. 
Models bf a serial or cumulative nature are best 
constructed by one person. 



ALLOCATION RULE 




(CALIBRATION! 

DATA I - SUBMODEL 2 

BASE 




MORTALITY FUNCTION 



POTENTIAL FUNCTION 



MODIFIER FUNCTION 



SUBMODEL 1 



SUBMODELS 



SUBMODEL 4 



Figure I. (a) Relation between data base and sub- 
models for current design showing a nearly par- 
allel structure, (b) Relation between data base 
and submodels for a hypothetical design show- 
ing a serial structure, 

(b) Efficiency is related to use of information in 
the calibration data base. 

Recall that a design feature of the model calls for 
each tree to be assigned to a particular class and to 
remain in that class, and for the size distribution of 
trees to be partitioned so that about the same num- 
ber of trees is in each class. When these two fea- 
tures are combined it is possible to use a multi- 
point boundary value problem approach to the es- 
timation of the modifier function. The source of 
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efficiencies here is that a boundary value approach 
does not require the taking of differences, often 
over uneven time intervals (see table 2, Christen- 
sen et al. 1979). Numerical analysts caution 
against taking differences of a noisy signal be- 
cause doing so accentuates the noise relative to the 
true signal. The regression approach, wherein dif- 
ferences are taken, admittedly used in the poten- 
tial, mortality, and allocation rule, ignores this 
caution, Keeping a small noise-to-signal ratio is 
very important in regions of the country where 
forest growth rates (i.e. , signal strengths) are rela- 
tively low. 

When the input tree list is partitioned by abso- 
lute size class limits, for example, 4.0-5.9, 6.0-7.9, 
etc., in a stand table projection, numerous stands 
to be projected may have no trees in the upper or 
lower classes. When the input tree list is parti- 
tioned on a percentile basis, each size class is cer- 
tain of having a number of trees present. This 
ensures that the calibration data base information 
available to develop the projection equation coeffi- 
cients is used on every plot. Figure 2 compares the 
traditional stand table and relativistic approach to 
tree list partitioning. 

(c) Efficiency is related to simplicity of the 
mathematical equations. 

When using differential/difference equations 
such as dY/dt or A Y/At=f(Y) to explain a pattern of 
Y over time, also given by Y=g(t),the rule is that f 
is algebraically simpler than g. In fact, quite sim- 
ple functions, f, can produce an extremely varied 
pattern of Y values over time. The use of systems of 
difference equations to handle the process com- 
plexity introduced by mixed tree sizes and species 
helps to keep the right-hand side of each equation 
as simple as possible. 



(d) Efficiency is gained by building into the 
model the response to treatment. 

The particular form of the right-hand side of the 
difference equations helps answer the question of 
stand response to a treatment calling for tree re- 
movals. This is made possible by having an upper 
limit of component change, given by the product of 
potential growth and number of trees, and a modi- 
fier of the upper limit, expressed as a function of 
stand density. The modifier has a form such that as 
density increases by growth, the value of the modi- 
fier decreases, but if density is abruptly decreased 
by cutting or mortality, the modifier is increased 
(fig. 3). 

(e) Efficiency of computations must also be con- 
sidered. 

Here the concern is the number of arithmetic 
operations that must be performed per unit of out- 
put information; for example, an updated tree list. 
The multiple resolution level design feature is re- 
sponsible for most of the gains here. If one wishes, 
the number of equations may be reduced to one (for 
a pure stand), and the tree list diameters updated 
every 10 years rather than every year. Doing this 
would have an effect on the precision of estimates 
of future forest conditions, but it would materially 
reduce the computational load to produce an up- 
dated tree list. The trade-offs between computa- 
tional load and precision of estimates have not 
been thoroughly examined at this time. 

(f) Efficiency in computer hardware require- 
ments should be a part of the design. 

The design of this projection system is such that 
one is not restricted to either a large storage capac- 
ity or extremely fast digital computer. Storage 
requirements can be reduced by such measures as: 
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1 2 3 4 56 78 9 10 11 12 13 14 15 16 1? 18 19 
D.B.H, (Inches) 

Figure 2, (a) The traditional (metric) approach to partitioning 
a sample distribution of trees. Notice that the 12. 9-16. 9-inch 
class has no trees, (b) The relativistic approach to tree list 
partitioning used in this Design, For illustration the group- 
ings are nearly *&, 1 A, *A 
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PROPORTION OF 
POTENTIAL 



MODIFIER INCREASED 
BECAUSE OF DECREASED 
DENSITY _[ 

MODIFIER REDUCED 
BECAUSE OF INCREASED 
DENSITY 



DECREASED DENSITY 
DY MORTAUTY AND 



INCREASED DENSITY 
BY GROWTH 



FigureS. Relation between changes in the propor- 
tion of the potential change, increased dejisity 
due to growth, and decreased density due to mor- 
tality and cutting. 



(1) projecting shorter tree lists (i.e., smaller plots), 
and (2) using lower resolution levels. Execution 
times can be reduced by employing both of these 
measures and by updating tree diameters and esti- 
mating mortality every 10 or so years (rather than 
every year). By a judicious combination of mea- 
sures to reduce both storage and execution time 
requirements, it is anticipated that the projection 
system can be made to function on both full-size 
and mini-computers. 



DIMENSIONALITY 

The fundamental premise in calling the higher 
dimensional models higher resolution models is 
based on the reduction in heterogeneity within 
class as class "width" is reduced. Theoretically, as 
heterogeneity is reduced, ability to predict growth 
is increased. In a practical sense, however, the 
increased predictability may not be a linear func- 
tion of number of classes (equations). In other 
words, going from a single equation for a two- 
species mixed stand to two equations may increase 
precision 20 percent, but adding two more equa- 
tions may not increase precision another 20 
percent. We are confronted with the costs of cali- 
brating and using whatever resolution level is 
developed. The determination of an optimum reso- 
lution level for growth modeling seems analogous 



to finding optimum sampling intensities wherein 
sampling intensities have precisions and coats as- 
sociated with each. 

The use for which a growth projection system is 
intended, it has been argued, determines the reso- 
lution level at which one should operate. Another 
consideration is the following relation between 
dimensionality of the model (equivalent to my res- 
olution level) and the time horizon over which 
projections are to be made: 



(11) 



where h is the time horizon over which pro- 
jections are to be made, 

d is the dimension of the model state 
vector, 

An is a constant for a given accuracy or 
precision level, and 

f is a function (Kahne 1976). 

For simplicity's sake, if f is assumed the identity 
function the relation is: 



d= 



k, 



h 



(12) 



Thus, the model dimension should be reduced 
for long-range projections and increased for short- 
range projections, to give the same level of 
accuracy. 
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The relation in equation (11) is apparently well- 
accepted by scientists involved in large-scale mod- 
eling. In general, however, it is not accepted by 
scientists working in forest projection. Too often 
projections are made for any time horizon of inter- 
est without regard to the dimensionality of the 
model. Sometimes single tree models have given 
unreasonable values after projections for 100 or 
200 years, and this fact has been used to discredit 
the entire approach. On the other hand, stand (sin- 
gle equation) models cannot be used for short-term 
evaluations of tree spacing relationships, and this 
fact has been used to criticize this well-established 
approach. Obviously there are needs requiring 
growth projection systems with a spectrum of reso- 
lution levels. 

In theory, it seems reasonable that one should be 
able to go from one resolution level to another 
without having to calibrate the denominator in 
the modifier function for every level. One would, of 
course, have to calibrate at several, say three, lev- 
els, but from then on one would simply input to the 
system the resolution level a user desires or can 
afford. The modifier function, as well as the other 
parts, would be computed by interpolation or ex- 
trapolation from the three calibrated levels. 



Once the initial calibration of the projection sys- 
tem has been completed there are at least three 
types of environments in which it may be used. 
One environment is where the objective is simply 
to project plot tree lists to a future time. This 
input-transformation-output schema is shown in 
figure 4. The algebraic forms come from the design 
phase, coefficients come from the calibration 
phase, and the input tree list may come either from 
temporary or permanent forest plots. No pro- 
jection system should be used in this manner for 
prolonged periods of time without also being used 
in the second environment. 

In the second environment accuracy and preci- 
sion of the projection system are evaluated. The 
schema for this environment is shown in figure 6. 
Here, in addition to input tree list, projection sys- 
tem, and output tree list, I have added the "experi- 
ence", detector and controller parts to the schema. 
This environment requires that permanent forest 
growth plot data be used for the input tree list and 




TREE LISTS . (TRANSFORMATION! PROJECTED TREE ^mi-miiT 
IITHE PROJECTION! usfs "UU I KU I 



Figure 4. Schema of environment for unverified 
projections. 




Figure 5.- Schema showing processor in a cyber- 
netic system involving use of permanent forest 
growth plots to correct the processor. 

that, after a number of years haa elapsed, a re- 
measurement of the plots be made and used as 
"experience". A comparison of the difference is the 
function of the detector, which passes the distribu- 
tion of differences back to a controller, which uses 
it to modify the algebraic forms and coefficients in 
the projection system. 

Reason suggests that with an untested pro- 
jection system one should not wait too long for 
another measurement an "experience" reading. 
It also suggests that when little deviation is de- 
tected after such a reading one may want to project 
for a longer period before the next check is made 
with "experience". Under conditions where exoge- 
nous variables are thought important, for exam- 
ple, climatic variations, it may be desirable to both 
fix an upper limit to the projection period and 
automatically use the most recent "experience", in 
place of the most distant, to recalibrate the model. 
This would produce a running average type model, 
as opposed to a cumulative model wherein all past 
measurements are included in the calibration data 
base. 
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A third use environment is shown schematically 
in what I have called the scenario environment 
(fig. 6). Here, the projection system is instructed to 
execute various management recommendations 
and guides, and the outputs of each set, each sce- 
nario, are examined in light of various perfor- 
mance indices 



Forest growth and yield research appears to be 
entering a new era where mathematical models 
are evaluated for generality of applicability as 
well as the traditional goodness of fit of equations 
to data. The basic model structure discussed in this 
paper is shown to be adaptable to single equation 
models as well as to models with large numbers of 
simultaneous equations. By using a relativistic 
approach to size class assignment the approach is 
made more computationally efficient than when 
fixed-diameter intervals are used. The concept of 
calibrating three resolution levels of the model 
and computing the needed coefficients in other 
resolution levels is presented, as is the notion that 
length of projection period should determine di- 
mensionality of the model. Variable resolution 
level, relativistic models will be more easily 
tailored to the individual needs and budgets of 
users than are today's models. 
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The need for growth projection systems that 
apply to all forest conditions of an ownership has 
been emphasized by Lundgren and Essex (1979), 
and a candidate projection system has been de- 
scribed by Leary (1979). Others have been de- 
scribed by Ek and Monserud (1974), Hegyi (1975), 
and Stage (1973). All the descriptions make it 
clear that a satisfactory model will not be devel- 
oped by analyzing a few growth plots. Rather, we 
need to have available a large amount of data. 
The data base for the current study covers 44 
different studies in Minnesota, Wisconsin, and 
Michigan. Included are 1,501 plots containing 
484,574 individual tree measurements on 92,649 
trees. 

The purpose of this paper is to describe the 
data base used to calibrate the growth projection 
system from several perspectives: space, time, 
age, site, and species. The data base used to val- 
idate the model is presented in Leary et al, (1979) 
<md the data needed to implement the model are 
escribed in Hahn and Brand (1979). 



SOURCES OF DATA 

Calibration of the growth projection system 
equires permanent growth plots with individual 
ree records at several points in time. Our sources 
vere plot records from (1) cutting experiments, 
i2) demonstration woodlots, (3) industrial contin- 
uous forest inventory, and (4) personal records of 
forest growth. In each case, the growth informa- 
tion came from permanent sample plots, rather 
than complete forest enumeration. The plots 
were all fixed size either circular or square. 
Square plots varied from 0.1 to 0.5 acre, and cir- 
cular plots varied from 0.1 to 0.2 acre. In some 
cases only trees greater than 9 inches d.b.h. were 



measured; for others all trees greater than 0.6 
inches d.b.h. were measured. Some plots were 
located randomly or systematically throughout 
the forest lands of an ownership. Other plots, 
coming largely from experimental forests, were 
all from similar quality sites. 

Measurement data from five properties were 
supplied to us on magnetic tape, and data from 
five other properties were supplied on cards. The 
data from the remaining 34 properties were sup- 
plied on the original field data collection forms. 
Inconsistent form layout between studies added 
considerably to the work of preparing field sheets 
for direct keypunching. Data transcribing was 
avoided as much as possible. 

During 1975 and 1976 each plot was visited and 
one or more of the following variables was mea- 
sured or estimated; d.b.h., crown class, crown 
ratio, and tree status. These data were added to 
the already existing data. 



DISTRIBUTION OF 
DATA SOURCES 

The spatial distribution of the data sources is 
shown in figure 1. Lower Peninsula of Michigan 
data are entirely from plantations of red and jack 
pine. Wisconsin data are mostly from natural 
stands of hardwoods oak-hickory type in the 
south and west, and northern hardwoods in the 
north and northern white-cedar. The two Wis- 
consin plantation studies are Star Lake red pine 
adjacent to Michigan and white spruce near the 
Menominee Reservation. The northern Minneso- 
ta studies include natural and plantation studies 
of red pine and jack pine, and natural stands of 
quaking aspen, white pine, and black spruce. 
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Figure 1. Location of counties with one or more calibration data sources in 
Minnesota, Wisconsin, and Michigan. 



The distribution of data sources by 1977 age 
(time) and site (matter/energy) is shown in table 
1. Each age-site table is for a different forest type. 
Although the total number of plots in table 1 is 
1,501, 1,529 plots were available for use; howev- 
er, the variable necessary for determining forest 
type was not available on 28 plots. The stand age 
at which plot measurements were used may be 
as much as 20 to 30 years less than the ages given 
in table 1. For example, the 26 jack pine plots in 
the 31-40 (age)-50 (site index) cell were 36 years 
old in 1977, They had, however, been measured 
since 1953 (table 2), so the measurements began 
when the trees were 12 years old. 

The distribution of the data sources by starting 
date and elapsed time between measurements is 
shown in table 2. The oldest measurements were 
made in 1938, giving nearly 40 years of observa- 
tions. Others are not far behind at 30+ years of 
records. Not every plot on a property was mea- 
sured where an "X" occurs in table 2. 



PROCESSING AND 
STORING DATA 

Information on field sheets was either key- 
punched directly or transcribed to 80-column 
sheets and then keypunched. The card data, as 
well as all data from magnetic tapes, were put 
into a standard format on disk files for editing. 
Following correction of all apparent data errors, 
the edited data were joined with the following 
information and placed in a master format: 

(1) Mode of stand origin, treatment, replication, 
property number, plot number. 

(2) Site index, site index species, and year of 
stand origin for up to five species per plot. 

(3) Soils information, 

(4) Forest type. 

(5) Denominator of plot size fraction. 

(6) Number of trees on plot. 

(7) Number of measurements in time, 
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Table 1. Distribution of 1,501 permanent forest growth plots used for model 
calibration by age and site within 13 forest types in the Lake States 
(Site class 30 ranges from 25 to 34, 40 from 35 to 44, etc.) 
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29 

157 




QUAKING ASF EH 


0-20 
21-30 
31-40 
41-50 
51-60 
61-70 
71-80 
Tolal 


1 1 2 
12211 7 
6251 1 15 
26631 18 
226 28 1 39 
167 14 
2 2 
97 




PAPER BIRCH 


21-30 
31-40 
41-50 
61-70 
91-100 
Total 


1 1 
1 1 
121 4 
1 1 

1 1 

a 











Table 2. Distribution of permanent growth plot measurements after the earliest 
measurement for calibration data base from Michigan, Wisconsin, 
and Minnesota. 
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Arjjoiine 'III!' 


oithem hnrduoodB 


47 


X 


> 


X 


y 


> 


X 


> 


X 






X 








X 


> 








y 


X 










X 






















(.airiin 'llil 


Hoi them hardwoods 


49 








j 








J 








> 










> 














J 
























f.llrivorth Till' 


Northern hardwoods 


18 










> 












X 




X 










> 
















X 






















Chil-icn Till 


lemlnck, Yellow 
birch, h'aple 


47 




X 


X 


X 


X 




X 




X 






X 








X 


X 


X 
















X 






















Ilundcp 'HIT 


Rod oaks, North- 
ern hardwoods 


17 






X 








X 






X 










X 
























X 




















llardii-H Ciuv\r THF 


Rud oak, White 


onk 


51 


X 












X 






X 


























X 

























Urcnlioldt Till' 


fliite pine 




52 














X 




X 


























X 






















HiuiBrtnkPi' IMF 


Jliitp !, Red pine, 
Aajien 


17 










X 




X 








X 








X 














X 










X 




















Chninptou Vnlley THF 


Northern hardwoods 


17 
















X 






X 














X 


















X 




















Stm I.nLe 


led i>lno 




11 




X 


X 




X 












X 














X 
















X. 












_x 










llutli-r'a Woods 


^nk, Dnnnwood, 
tlnple 




15 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 










X 




-g 




_ 




NcttLf'ton 


Jortliern liflidwo 


oils 


38 




M 








K 




J 


> 


X 










X 


X 






X 




X 








> 


X 








Stone 1 '] Hoods 


3nk, Hickory, 
Maple 




45 


X 


X 


_x 

X 


J 

X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 






















X 
















Hhim|'iH't!'jt-> 'n Hoods 


Ink, Hickory, 
Maple 




15 


X 


X 


X 


XI 


xl 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 






















X 
















HC-2f> 


Inlenn fir, Aopcn 


49 










Jl 


X 








X 




X 








X 




X 




X 


X 


X 


X 




X 
























Ull- 25 


lorthcrn liardwoodE 


51 












> 








X 










X 










X 






X 




























KH-11U (Arjjomip) 


Northern hardwoods 


56 












X 




X 








X 






X 






X 






































Nil-lift ftledfnrd) 


Northern hardwoods 


58 








X 






X 








X 










X 










































:l('-lfl 


Jhite epiuce 


19 












X 


X 
















X 




















y 
























St-li 


. Hliite ccdor 
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l An "X" denotes that at least one plot on the property was measuredat the number 
of years after the earliest measurement given by the column label. Only Butler's, 
Stone's, and Shaugnessey's Woods were measured in their entirety every year. 
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(8) Dates (month "and year) of each measure- 
ment, 

The master format consists of a maximum of 
512 word blocks of binary information. The block 
length is specified because of input/output buffer 
length on the Control Data Corporation's Cyber 
74 computer. Each word can contain up to 14 
integer values to code information. The master 
format and contents of each word are given in 
figure 2. An example output of a permanent 
growth plot record is shown in figure 3. Specific 
objectives of the master format were to (1) max- 
imize efficiency in reading/writing and storing the 
master files, and (2) have information in a stan- 
dardized format because five or six different per- 
sons may be getting information from the master 
at any one time. 



RETRIEVAL OF INFORMATION 
FROM DATA BASE 

Because the information is written in blocked 
nary form with 14-character words, one must 
se simple FORTRAN statements to extract the 
jsired values from a word in a block. For exam- 
ie, to analyze mortality, the master file is read 
.nd the desired information extracted to estab- 
lish an analysis-specific file. Most of the sort word 
contents are normally carried into the analysis- 
specific files. This information, together with 
other information, aids in tracing any problem 
observations and for organizing analyses by 
state, species, mode of stand origin, etc, The anal- 
yses then proceed species-by- species for the 26 



species groups for which functions are devel- 
oped. 



Word 
1 
2 



3-7 



10-9+meo 



10-hnea, 
ll+2mea , 



11-hnea- 
12+2raea 



Number of voids In record 
Sort word 

II Mode of stand origin 

11 Treatment 

II Replication 

13 Data source (property number) 

15 Plot number 

11 Record number (1, or 2, or 3, or '* 1C more 

than 500 words on preceding record 
(except for record 1)) 

Site index (up to 5 Sl-aperles combinations) 

IZ SI 

13 Speclea for SI 

13 Stand origin (last 3 digits of origin year) 

Plot invarlate parameters 

12 Foreat type 

12 Plot size (aa denominator of fraction, i.e., 

5=1/5 acre plot) 

13 Number of trees on plot 

12 Number of measurements = mea 

Soils Information 

Plot variable parameters 

14 Date of measurement (MMVY) 

13 BA/acre at time of measurement 

Tree invarlatG parameters 



13 Tree number 

13 Tree species 

12 Tree species group 

Tree variable parameters 

II Status 

II Condition 

13 DBH to last tenth Inch (F3.1) 
11 Crown class 

II Crown ratio 

11 Tree quality 

12 Form factor 

13 Height 



Figure 2. Lake States FREP binary master rec- 
ord format. 



92 1000100000011 A1L059LB Ar3/5941i 


10/01604 




J 059026 


1064037 


1 06 904 H 


I0/A0h3 


110501 


4009800000000 
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1011900000000 


10i:9:>6000064 
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J. 00 7200000000 


1008400000000 
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1014000000000 


10X4700000000 
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LOt LLOOOOOOOO 
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I01J/00000000 
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1011700000000 
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4005900000000 


J. 006900000000 


3007900000000 


10090^6000000 


837526 


4006000000000 


1007000000000 


1008000000000 


1009,1313000000 


909506 





4005,^00000000 


1006000000000 


1006/3HOOOOOO 


1037526 





4 005 MOO 0000 00 


1006300000000 


100/1 ''6000049 


1 1 09506 





4005600000000 


L006500000000 


I00/?2il000040 


1209506 





4005700000000 


1006900000000 


I 00/928000000 


131050J 








4006100000000 


LOOA034000000 


1437526 








4005900000000 


I007M6000000 


1&37526 








4oos:'oooooooo 


toor;H4Ei00oooo 


J 6 09506 











4005450000000 



Figure 3, Sample output listing of a single plot from the mas- 
ter file. Contents and meaning of each word are 
given in the master format in figure 2. 
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Jerold T. Hahn and Rolfe A. Leary, 

Principal Mensurationists 



When seeking an unknown value of a numerical 
variable it is common practice to bracket the un- 
known by two values one greater than the 
unknown and the other less, Mathematicians have 
used this practice for many years to find the solu- 
tion of nonlinear algebraic equations. A simple 
example is the problem of finding the value of x 
that satisfies the equation x =y, when y=142. 

The square root of a number can be estimated by 
guessing a number such that its square is greater 
than y and then guessing another with square less 
than y; e.g., 12 2 =144 and II 2 =121. The desired 
value is then known to be between 12 and 11. By 
judicious guessing one quickly gets several correct 
digits in the value of x. 

The idea of establishing an upper and a lower 
bound to the rate of a biological process is a very 
useful modeling device. In the case of the accretion 
processes it is reasonable to take zero as the nu- 
merical value of the lower bound. This paper is 
concerned with estimating the largest possible nu- 
merical values the potential for increase in di- 
ameter at breast height (d.b.h.) by all major forest 
species in the Lake States. Our concern is with 
both the potential diameter growth rates of indi- 
vidual trees and the diameter growth potential of 
all trees in a stand. 

The concept of a growth potential has been used 
in mathematical modeling of biological processes 
for many years (Verhulst 1965, Lotka 1956, Vol- 
terra 1931, Gause 1969). The potential is given by 
r in the generalized differential equation of popu- 
lation growth: dY_ 

dt 

The term r refers to the rate of self-renewal of the 
population in question. Its effect dominates the 
right-hand side of the equation as long as re- 
sources are not limiting. 



Foresters have only recently begun using the 
concept of potential in forest growth models. Ar- 
ney (1974), Hegyi (1974), Ek and Monserud (1974) 
have used potential height growth as the upper 
limit of growth. Potential height growth is usual ly 
estimated by taking the first forward difference, 
with respect to age, of equations to predict height 
based on site index and age (see Lundgren and 
Dolid 1970). The use of breast height diameter 
growth to characterize a tree's potential is less 
common than the use of height (Newnham 1964, 
Chen 1976). A disadvantage of this approacH is 
that there are no standard functions relating tree 
and environment variables to diameter growth , as 
there are to relate tree and environmental varia- 
bles to tree height. 

In theory, the following tree and environment 
variables affect the amount of lateral stem grow tli 
at breast height: (1) current d.b.h., (2) amount of 
foliage, (3) respiring surface area, (4) distance 
from breast height to center of live crown, and (5) 
abiotic flux as reflected by site index. The relation 
between each of these variables and diameter 
growth was extracted from the literature by 
Brand: 1 

(1) The closer the crown to breast height, trie 
more stem growth (Duff and Nolan 1953), 

(2) Le-s foliage is needed on a good site (Ass- 
mann 1970). 

(3) Amount of foliage is greater on a good site 
(Long and Turner 1975). 

(4) Respiring surface area is related allometri- 
cally to d.b.h. (Pienaar and Turnbull 1973). 



l Brand, G. 1976. Personal communication, St. 
Paul, Minnesota. 
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There are many ways of combining variables 
used to characterize each of these relations. Two 
broad categories are strictly additive association 
and strictly multiplicative association. Of course, 
one may have combined additive and multiplica- 
tive. 



Species 



Number of trees 



Our objective was to develop species-specific di- 
ameter growth equations that characterize the po- 
tential of Lake States tree species to increase in 
d.b.h. We have used a combined additive-multipli- 
cative model that requires knowledge of the fol- 
lowing variables: d.b.h., crown class, crown ratio, 
and site index (age 50 base). Crown ratio is the 
ratio of full live crown length to total tree height. 



Data Base 

Individual tree diameter measurements at two 
points in time were collected on dominant and 
codominant trees of 26 species in Michigan, Wis- 
consin and Minnesota. A total of 51,149 trees came 
from permanent growth plots or from forest survey 
plots in Wisconsin; of these, 36,726 were high- 
quality dominant and codominant trees usable for 
estimating potential growth. 

All trees had crown class and crown ratio esti- 
mates at the last measurement. Plot site index was 
also estimated. The geographical distribution of 
plots, and the distribution in the age-site matrix is 
given in Christensen et al. (1979). 

The diameter growth (AD) for each tree was 
determined over a period approximating 10 years. 
This was the dependent variable in a nonlinear 
regression analysis. The independent variables 
were initial diameter (D), final crown ratio (CR), 
and plot site index (SI). These data for all domi- 
nants and codominants were grouped by species, 
as shown in the tabulation above, by 1-inch diame- 
ter classes, by 10-foot site index classes, and by 10- 
percent crown ratio classes, If a cell contained just 
one tree, it was included in the most appropriate 
adjacent class. For each cell the means were com- 
puted for each variable, i.e., AD, D, SI, CR, as well 
as the standard deviation of diameter growth, SAD- 



Jack pine 


4,625 


Red pine 


11,663 


White pine 


845 


White spruce 


1,772 


Balsam fir 


1,277 


Black spruce 


1,006 


Tamarack 


165 


Northern white-cedar 1.070 


Hemlock 


140 


Black ash 


240 


Cottonwood 




Silver maple 


937 


Red maple 




Elm 


898 


Yellow birch 


343 


Basswood 


1,092 


Sugar maple 


3,753 


White ash 


866 


White oak 


517 


Select red oak 1,555 


Other red oak 290 


Hickory 


349 


Bigtooth aspen 261 


Quaking aspen 2,124 


Paper birch 


748 


Other hardwoods 190 


Total number 


of trees 36,726 



The mean diameter growth, AD, for each cell 
applies to dominants and codominants. Since it 
was our desire to estimate the potential of trees to 
increase in diameter at breast height, it was as- 
sumed that the diameter growth values were nor- 
mally distributed around this mean (see fig. 1), 
The potential is, of course, greater than the aver- 
age. We took the potential to be the mean diameter 
growth plus 1.65 standard deviations 
(AD-t-1.65s iD ). This diameter growth corresponds 
to the 95th percentile of dominant and codominant 
diameter growth. 



Mathematical Model 

The following four mathematical models were 
examined. All were nonlinear in form and com- 
bined additive-multiplicative with a core Rich- 
ards-type function that has anabolic and catabolic 
terms (Richards 1959): 
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Figure 1. The assumed normal distribution of 
diameter growth of dominant and codominant 
trees for a 10-foot site index class and 1 0-percent 
crown ratio class. The shaded area represents the 
fastest-growing 5 percent of the trees. 



AD 



AD 



AD 



AD 



+b 4 CRD 



D 



+b 4 SICRD 



+b,SICRD 



CD 
(2) 

(3) 
(4) 



The pattern of each is -vr= ["intercept" term] + 

[catabolicterm] + [anabolic term]. The signs onbi, 
b 3 , b.i, and b s are positive, and b 2 is negative. The 
"intercept" term gives the growth approached as 
d.b.h. approaches zero; in equations (3) and (4) 
lim AD 

D 



_ b 

'' 



At 
This limit is a function of SI in (1) and of CR in 

b 3 

(2). The catabolic term, b 2 D , has the form of an 
allometric relation, as required of a relation char- 
acterizing respiring surface area (Pienaar and 
Turnbull 1973). Several anabolic terms were 



tested. The terms in equations (3) and (4) make the 
most sense, because total tree height is multiplied 
by crown ratio, giving crown length. Stem diame- 
ter has often been related to crown diameter. The 
product of crown length and width gives an ap- 
proximation, uncorrected for crown shape, of 
crown volume. 

Fitting Procedures 

Each of the above models was fit to both mean 
and potential AD values for each cell of each spe- 
cies using Marquardt's algorithm modified for 
weighted nonlinear least squares regression (Mar- 
quardt 1963). Model evaluation criteria were: 
qualitative behavior of the solution to each equa- 
tion using the estimated parameters bi; pattern 
formed by residuals (predicted minus observed di- 
ameter growth) plotted against initial diameter, 
crown ratio, and site index; and goodness of fit as 
measured by R . For some species the qualitative 
behavior of the solution was unrealistic because 
an insufficient number of observations was availa- 
ble at large diameters. To compensate, the silvi- 
culture literature was consulted and from it we 
determined the upper limit of d.b.h. for the speciea. 
A very small AD value was assigned to the diame- 
ter, an average crown ratio and site index were 
attached, and the "add" values were inserted into 
the data base as needed to insure more realistic 
behavior of the function at large diameters. 

Results 

Models (3) and (4) gave better results than ( 1) or 
(2) as judged by model selection criteria. The re- 
sults of the regression analyses of models (3) and 
(4) for 26 species for both average growth and 
potential growth were examined and one or the 
other was selected as the best model. The coeffici- 
ents, number of cells, and R 2 values are given in 
the tabulation. When the b 5 parameter is exactly 
1 .0, the model selected was number (3). Otherwise, 
model (4) was selected. Figure 2 shows potential 
diameter growth values for species commonly 
found in the pine type; red pine, white pine, jack 
pine, and quaking aspen. The crown ratio was set 
at 4.0 and site index at 65. 

DISCUSSION 

To use the potential diameter growth function it 
is necessary to know the initial diameter, crown 
ratio, and site index. Crown ratio is not a com- 
monly collected tree variable, but the average 
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20 30 40 

D H H CLASS (INCHES) 



Figure 2. -Predicted potential diameter growth 
values for species commonly found in the pine 
type in the Lake States. Site index is 65 and 
crown, ratio is 4.0 for each species. 



crown ratio for a stand can be estimated by know- 
ing mean stand basal area, mean diameter, and 
stand age (Holdaway et al. 1979). 

The inferential population for equations (3) and 
(4), with coefficients in table 1, is the natural 
stands and plantations of the 26 species in Michi- 
gan, Minnesota, and Wisconsin. No attempt was 
made to separate trees in plantations from those in 
natural stands, because the individual tree is not 
thought to behave differently. 

The approach used for potential diameter 
growth is somewhat analogous to that used for 
models based on potential height growth. The lat- 
ter use height-age relations for dominants (Bite 
index curves). We used diameter growth-diameter 
relations for the fastest-growing dominants and 
codominants. By taking the 95th percentile of di- 
ameter growth, we intended to approximate diam- 
eter growth rates of forest trees not subject to sig- 
nificant competition from neighboring trees. We 
judged open-grown trees as not suitable for poten- 
tial diameter growth studies because of the signifi- 
cant difference in the distribution of increment 
among bole, branches, and limbs. 



Table 1. Number of cells and coefficients for average and potential diameter growth functions for 26 species in 
the Lake States l 



Number of 


Average 


Potential 


Species group 


groups 


b, 


b, 


b, 


b. 


b. 


R' 


b, 


b, 


b a 


b, 


b, 


R' 


1 Jack pine 


126 


0.09349 


-0.00005 


2.9144 


000003 


1.0 


044 


0.16062 


-0000009 


36245 


000004 


1.0 


035 


2 Red pine 


282 


.05309 


-00013 


20468 


,0001 97 


38337 


,64 


09446 


-0,00012 


2.0596 


00035 


.24225 


54 


3 White pirns 


161 


15242 


-.00030 


19660 


.00003 


1.0 


.34 


25578 


-.00058 


1.7263 


00004 


1.0 


28 


4 While apt UCB 


72 


18789 


-00776 


13904 


.00006 


1.0 


.51 


.17056 


- 014516 


1M60 


,00052 


27298 


74 


5 Balsam fir 


121 


037102 


-.00026 


19708 


.00016 


37750 


.64 


12200 


-.00080 


19390 


,00006 


1.0 


.50 


G Black spruce 


105 


.06667 


-00048 


21592 


00003 


.95715 


.26 


.10713 


-.00107 


20017 


.00006 


9(127 


,39 


7 Tamarack 


56 


05279 


-.000002 


34307 


.00003 


10 


.29 


11147 


-.00001 


30685 


00003 


10 


22 


8N white-cedar 


S3 


.06480 


- 00000002 


4.7468 


00012 


15411 


41 


.13403 


-.000001 


3,6880 


00002 


1.0 


17 


9 Hemlock 


59 


.11278 


-.000000009 


4.3699 


000007 


1.0 


25 


t6872 


-,0000003 


35738 


00001 


1,0 


25 


1 1 Black ash 


84 


.035027 


-000000001 


53434 


,00010 


46385. 


,41 


05BB07 


-.OOOOOOOOS 


50559 


.00024 


31553 


.26 


12 Cottonwood 




























13 Silver maple 


176 


.06793 


-00038 


1.5944 


00020 


.29971 


.44 


.10948 


-00004 


2226 


.00050 


06626 


30 


1 4 Red maple 




























15 Elm 


179 


148B2 


-.00016 


1,9910 


.00002 


1,0 


.12 


.28498 


-.00182 


1.5297 


00003 


1,0 


00 


16 Yellow birch 


96 


11333 


-00005 


2.3334 


000006 


10112 


.12 


.15155 


-.000003 


3,3104 


00007 


57302 


.16 


17 Basswood 


169 


.16782 


-00026 


1.8136 


000008 


1.0 


16 


.25402 


-.000004 


2,9396 


00001 


1.0 


09 


16 Sugar mapla 


272 


.13309 


-.00051 


17002 


.00002 


10 


29 


.1&772 


-.000007 


25839 


,00028 


08365 


,33 


19 While ash 


139 


.10829 


-.00025 


2.0926 


.00003 


1.0000 


.36 


21167 


-,000003 


33131 


.00001 


12430 


,27 


20 White oak 


152 


08114 


-0000004 


3.1795 


00004 


,61666 


.26 


,12654 


-000004, 


2.7538 


00005 


.62066 


21 


21 Select red oak 


288 


098793 


-0000002 


3.3873 


00008 


.41156 


.38 


.15535 


-00000010 


35367 


000t8 


25900 


22 


22 Other red oak 


92 


10142 


-0000002 


3.7728 


00012 


.20604 


,46 


17358 


-00007 


2.4451 


00003 


.95222 


30 


23 Hickory 


93 


.071906 


-.000008 


27884 


,00003 


69376 


44 


.16471 


-.02161 


13949 


00004 


.71628 


,26 


24 Big tooth aspen 


79 


,1-1120 


-.00029 


2.1217 


.00003 


10 


.28 


,23490 


-.00942 


11041 


00036 


,16368 


.33 


25 Quaking aspen 


201 


.11789 


-.00015 


2 361 8 


.000042 


1.0 


,39 


,21646 


-.000091 


26030 


.000044 


1.0 


,31 


26 Paper birch 


112 


.07565 


-.00059 


2.0175 


.00005 


10003 


.48 


1*971 


-.00032 


20236 


.00034 


.21268 


.41 


30 Other hardwoods 


78 


.10662 


-.00954 


10500 


.00002 


1.0 


.44 


.37510 


- 02479 


9,6268 


00004 


1.0 


.36 


Tolal 


3,284 



























bs 



AD a 

1 Crown ratio variable- used Is 3.0 rather than 0.3. General model ls-^-=b t 4-b 2 D +b<SI CR D 

At 



25 



To estimate stand potential growth, the mean 
stand d.b.h., mean crown ratio (either observed or 
predicted), and site index are substituted into 
equation 4, along with the appropriate coefficients 
from table 1; the potential growth of a single tree is 
then computed. This is the estimated rate of 
growth of a free-to-grow tree. Potential stand 
growth is simply the number of trees in the stand 
times the free-to-grow rate. 

CONCLUSION 

Potential and average diameter growth rates for 
26 species in the Lake States are adequately pre- 
dicted by combining initial tree diameter, crown 
ratio, and site index, into a nonlinear model con- 
taining three terms: intercept, catabolic, and ana- 
bolic. The intercept term is a species-specific 
constant, the catabolic is an allometric relation of 
tree diameter, and the anabolic term is a product of 
crown length and an allometric relation of tree 
diameter. Since the potential diameter growth 
function characterizes free-to-grow conditions, po- 
tential stand diameter growth is simply the prod- 
uct of tree potential and number of trees. 
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The amount of foliage in a stand of trees is an 
important indicator of the growth one may expect 
from the stand. Although this relation has been 
known for decades, the amount of foliage is not an 
easily determined variable that can be used in 
equations to predict stand growth. Field measure- 
ments on permanent growth plots have tradition- 
ally concentrated on stem dimensions (i.e., diame- 
ter and height), with supplemental correlation 
studies to associate crown diameter with stem di- 
ameter and other stand variables. An exception is 
the national forest survey, which has included 
measurements on tree crown ratio for many years. 
Recall that crown ratio is the fraction of total tree 
height in full live crown. Throughout this paper 
we express the "fraction" as an integer between 
and 10 i.e. , a crown ratio of 2 means a tree has 1 1- 
20 percent of total height in full live crown. 

Crown ratio was found to be an important varia- 
ble for predicting tree growth (Hahn and Leary 
1979), and was incorporated into an equation to 
predict potential diameter growth in a stand. Be- 
cause the potential growth equations require 
mean stand crown ratio as an independent varia- 
ble, and some prospective users of the generalized 
growth projection system (Hahn and Brand 1979) 
may not have included this variable in their plot 
data, it is necessary to develop equations to ap- 
proximate mean stand crown ratio. 



MODEL RATIONALE 

From a theoretical analysis it is suggested that 
the mathematical model should exhibit the follow- 
ing qualitative properties: 

(1) solving the mathematical equation must 
give a value between and 10, 



(2) crown ratio must not increase as BA (basal 
area per acre) increases, 

(3) crown ratio must ultimately decrease as age 
or average diameter increases, and 

(4) the rate of change of crown ratio with basal 
area must decrease as age or average diame- 
ter increases. 

The exact relations, rate coefficients, etc. appear 
to be species-specific. 



The core model tried initially was: 
b, 



CR= 



l+b 2 BA 



(1) 



Although this model satisfied condition (2) 
above, it did not satisfy (1) and, because only basal 
area is used, it did not satisfy the important condi- 
tions in (3) and (4). The additional variables of age 
and average diameter were used with equation (1) 
to produce models satisfying conditions (.1) 
through (4). 



METHODS 

The initial data base (Christensen et al. 1979) 
consisted of 1,629 permanent growth plots in the 
Lake States, Based on methods of treating plots of 
mixed species (Leary and Holdaway 1979), each 
plot was assigned as a pure species plot or as a 
combination of up to three species. In this study 
only plots fitted as "pure" or as containing two 
species were used. The individual tree crown ra- 
tios, measured on each plot in the last measure- 
ment year, were averaged to obtain a species mean 
crown ratio (CR). Plot basal area was averaged 
over the 10 years preceding the crown ratio mea- 
surement. When no measurement had been taken 
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10 years preceding, the measurement nearest to 
10 years was used. Only plots with no major distur- 
bances (i.e., losses to mortality and/or cutting) dur- 
ing the interval were included. If the final basal 
area was less than 90 percent of the initial basal 
area, a major disturbance was assumed to have 
occurred and the plot was not used. For all accept- 
able plots mean stand diameter (AD) was calcu- 
lated at the final measurement year and the stand 
age (AGE) was calculated at the midpoint of the 
interval. 



Because of insufficient or poorly distributed 
data for some species, it was necessary to group 
species as follows: 

Group number Species group 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 

12 



jack pine 

red pine 

white pine 

white spruce, balsam fir 

black spruce, tamarack 

northern white-cedar, hemlock 

cottonwood, silver maple, red maple 

black ash, white ash, elm, basswood 

yellow birch 

sugar maple 

white oaks, select red oaks, other 

red oaks, hickory 
bigtooth aspen, quaking aspen, paper 

birch 



The analysis showed two patterns of relations 
one for conifers and another for hardwoods. The 
conifers tended to have much larger crown ratios 
for low basal areas than did the hardwoods. Mean 
plot crown ratios for conifers were seldom greater 
than 8 or less than 2. The range for hardwoods was 
much smaller; crown ratios seldom exceeded 5 
even at very low basal areas. Fourteen different 
models, each containing the core relation, were fit 
to the data from 12 species groups using nonlinear 
regression methods (Marquardt 1963). The follow- 
ing model was selected for conifers: 
b.AGE 



, where 



(2) 



AD 



BA -basal area in square feet, 
AGE=years from stand establishment, and 
b,, b 2 - unknown numerical constants. 

Equation (2) coefficients for each conifer species 
group are as follows: 



Species group 

Jack pine 

Red pine 

White pine 

White spruce- 
balsam fir 

Black spruce- 
tamarack 

Northern white- 
cedar-hemlock 



b, 



Standard error Number 
ol estimate ol plots 



-0.0187 0849 38 35 123 

-,01270 .0504 .42 4,3 384 

- 0060 1399 42 2 1 62 

-,0058 .0815 71 37 bZ 

-.0039 .1045 .12 70 41 

-.0213 .0055 .41 37 42 



The amount and rate of change of crown ratio 
with age, average diameter, and basal area for 
hardwoods was quite different from that for coni- 
fers. In general, hardwood crown ratio changed 
very little over the range of stocking, size, and age 
available in the data base. The model for 
hardwoods was established as the following varia- 
tion of the core model: 



CR=b, AGE-f 



(3) 



AD 



Coefficients in equation (3) for the hardwood 
species groups are: 



Species group 


Stanford enar Number 
b, b, b, R 1 ol estimate ol plots 


Cottonwood, silver 


















maple, red maple 


-0 0186 


6 


.910 0. 


.0476 46 


2 


,2 


27 


Black ash, elm, bass- 


















wood, white ash 


0000 


3 


04 , 


,0000 


NA 


1 


4 


64 


Yellow birch 


-0222 


6 


.019 , 


,0264 


.37 


2 


6 


16 


Sugar maple 


-.0064 


4 


707 , 


,0282 


.26 


1 


9 


143 


White oaks, select 


















red oaks, other red 


















oaks, hickory 


-0080 


4 


.728 , 


,0381 


.32 


1, 


,2 


112 


Bigtooth aspen, 


















quaking aspen, 


















paper birch 


-.0070 


4 


.298 , 


,0230 


.28 


2 


,3 


112 



AD=arithmetic mean d.b.h., 



The crown ratio values for the 64 plots of black 
ash, elm, basswood and white ash showed so little 
variation that the mean value was assigned to the 
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tent b 2 (equation 3), and b i and b 3 were set at 
ornparing equations (2) and (3) we see that 
ents b 2 (equation 2) and b 3 (equation 3) op- 
Lo achieve the same effect; reducing the 
crown ratio as basal area increases. The 
& rate of change for conifers is twice that for 
jods. Another major difference between spe- 
>xips is based on the "intercept" values i.e., 
iwn ratio as basal area* age, and average 
er approach zero. For conifers, the "inter- 
3 10, a full live crown. For hardwoods the 
pt is b 2 , and ranges from 6.9 for cottonwood, 
pie, and silver maple to 3.04 for black ash, 
Lsswood, and white ash. Mean stand crown 
:>ver a range of average diameters, basal 
and ages are shown for red pine and sugar 
in figures 1 and 2. 









12 



20 



BASAL AREA PEn ACRE fsOUAHE FEET) 

1. Mean stand crown ratio for red pine 
s of selected ages and average diameters of 
and 20 inches. 






12 



so 120 leu 

BASAL AREA PER ACHE (SOUAItE FEET } 



Figure 2. Mean stand crown ratio for sugar ma 
pie stands of selected ages and mean diameters Oj 
4, 12, and 20 inches. 



DISCUSSION 

Adequate estimates of mean crown ratio for 
stand can be made if one knows the approxima 
mean basal area per acre over the last 10 years, th 
mean stand diameter, and the stand age. Bette 
estimates can be made if one knows, say, the rui 
ning average 10-year basal area. Simple poi 
estimates of stand basal area will produce p 
crown ratio estimates if one knows nothing abi 
prior stand disturbances. A stand disturbance 
less of a problem if the stand is old or the In 
large. For examplef if a red pine stand with i 
average d.b.h. of 6 inches, age of 40 years, ai 
basal area of 150 square feet is cut back to 1C 
square feet with no change in average stand diam 
eter, the predicted mean crown ratio would in 
crease from 0.30 to 0.37. However, if the mean 
diameter had been 28 inches, the age 200, and the 
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basal area 150 square feet, cutting back to 100 
square feet with no change in average stand diam- 
eter would have had no significant effect on the 
mean crown ratio. When applied in the projection 
program (Hahn and Brand 1979), the running 
average basal area for the past 10 years is used 
instead of an average of only two or three values. 

If cutting or mortality changes average stand 
diameter as well as stand basal area, the following 
relation holds for the models: a decrease in diame- 
ter as well as basal area (e.g., thinning from above) 
results in a lower crown ratio, whereas an increase 
in diameter and a decrease in basal area (e.g., 
thinning from below) results in a higher crown 
ratio. This relation is logical since larger trees of 
the same species in a stand (dominanats and 
codominants) tend to have higher crown ratios 
than smaller trees (intermediates and sup- 
pressed). 

A generalized relation between mean crown ra- 
tio, age, mean d.b.h., and basal area, applicable to 
both conifers and hardwoods, is given in equation 
(4): 

b a AGE 
b 2 e +b,i 



CR=b,AGE+ 



(4) 



1+ i6 BA 



For conifers, b ( is set at zero; for hardwoods b 2 is 
set at zero. 
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nise in the design logic for this gen- 
growth projection system is that the 
ag-gregate can be represented by a 
combination of a potential growth 
of the potential. The modifier takes 
mathematical function since another 
t the amount the potential is modi- 
i the amount of standing crop pres- 
oses of this paper are to describe 
results of studies to determine how 
3 interact within and between spe- 

e modifier function has in the basic 
m is shown in equation (1): 

Number 

of [Modifier]. (1) 

trees . 

n the analysis we assume that the 

leter growth function is known, 

difier function as the only function 

numerical constants. The relation- 

lifier function to the data base and 

is is shown in figure 1. Because the 

ion. is assumed known, estimation 

cannot begin until it has been com- 

,hown connected to the potential 

crown ratio function. Under ideal 

here would be no need for a crown 

in the calibration phase because 

ild have been observed at all mea- 

ice actual observations could be 

idy crown ratio was estimated only 

isurement, in 1975 or 1976; thus, if 

nt interval selected for estimating 

is from a prior time, it was neces- 

e the crown ratio used in the po- 

i. Because the modifier function 

pend on an estimated potential 

i turn may depend on an estimated 




MORTALITY 

ALLOCATIONS. _ X FUNCTION 

RULE 

(CROWN RATIO 
FUNCTION) 

I 

POTENTIAL 
FUNCTION 

MODIFIER 
FUNCTION 

Figure 1. Relation of modifier function to data 
base, potential function, and crown ratio func- 
tion. The latter is shown with parentheses and 
dashed line because when the crown ratio has 
been recently observed it is not necessary to use 
the crown ratio function in estimating numerical 
constants in the modifier function. 



crown ratio, the modifier function's sensitivity to 
variations in atand conditions may be less than it 
could be. 



REQUIREMENTS OP THE 
DATA BASE 

Unlike the potential and mortality functions, 
the modifier function is based on information 
about all the trees on a plot. Thus, the data re- 
quirements relate to plot attributes rather than 
individual tree attributes. For a plot to be useful 
for estimating modifier function coefficients it 
must: 

(1) have been measured three or more times, 

(2) have been measured over approximately 10 
years, 

(3) exceed a minimum number of trees, 

(4) exceed a minimum basal area density, 
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(5) be either a pure stand or a two-species mix- 
ture, and 

(6) have not had excessive cutting or mortality 
in the measurement period. 

The exact cutoff values for these six criteria are 
given in the methods section. Many plots were not 
used to estimate the modifier because they failed 
to meet one or more of the criteria. 

METHODS USED TO 

SUMMARIZE, STORE, 

AND RETRIEVE DATA BASE 

INFORMATION 

Before plot data are ready for analysis, the ana- 
lyst must: 

(1) Take the data, as found in the database, on the 
species present on each plot and transform 
them into a condensed form, recording the 
transformation. 

(2) Use this transformation information to consol- 
idate all species present into a maximum of 
four species categories and then summarize 
the needed plot information into one physical 
record to enable the mass of data to be effi- 
ciently stored on disk files. 

Design a system that enters a large file of 
information produced by step 2 (above) and 
picks out only those plots that match a speci- 
fied species combination, "unpeels" the data, 
reorders the species if necessary, and lists the 
unpeeled data on a new file in a format ready 
for analysis. 

Step 1 (above) needs to be discussed in more 
Letail. We found in a preliminary study that the 
,hree needed measurements could be taken over as 
short an interval as 4 years without affecting the 
results. Whenever possible, however, we tried to 
have at least 6 to 8 years in the interval from first 
to third measurement. A computer program was 
written to locate where the cuts, mortality, and 
ingrowth occurred. From these results we manu- 
ally selected an interval of at least 6 years with a 
minimum of disturbance due to death, cutting, or 
ingrowth, and a maximum number of live trees. In 
accommodating these constraints, the three-mea- 
surement interval tended to be as short as possi- 
ble. This enabled us to use more than one interval 
from plots with five or more measurements. Using 



the chosen interval(s) we determined the percent 
basal area for each species on the plot, ordered the 
species by basal area, retained the first three spe- 
cies, and combined the remaining species into a 
miscellaneous category. 

It was still necessary to reduce these species 
down to the principal ones. To accomplish this f a 
predominant species composition fitting rule was 
developed. The rule works approximately as fol- 
lows for a plot with only two species: the plot is 
located in a coordinate system at the point (x^Xa), 
where x f is the percent of plot basal area in species i 
(fig. 2). If the plot were a pure stand of species 1, 
the point would be located at coordinates (1,0); if it 
were a pure stand of species 2, it would be located 
at coordinates (0,1). All other mixtures would be 
located on a line between points (1,0) and (0,1) 
auch that Xi+x=l. For a two-species plot, we cal- 
culate the distance from (X U X K ) to three points, 
(1,0), (1/2,1/2), (0,1), and assign the plot to that 
species combination (pure species 1, two-species 
mixture, pure species 2) for which the distance ia 
smallest. 



100 



pure species 2 




(1,0) pure spocEei 1 



1QQ 

PERCENT SPECIES 1 

Figure 2. Schematic diagram of procedure for 
judging the composition of a plot containing two 
species. In this example d t is less than rf ia hence 
the plot will be considered a pure stand of species 
1 , Analogous methods are used in 3 and 4 dimen- 
sions for 3 and 4 species mixtures. 
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In our work it was possible to have a maximum 
of four species categories. Since the species are put 
in descending order by percent basal area, the 
reference points for determining species fits are 



(1,0,0,0) 

(1/2,1/2,0,0) 

(1/3,1/3,1/3,0) 

(1/4,1/4,1/4,1/4) 



corresponding 

to 



pure species 1 
two-species mixture 
three-species mixture 
four-species mixture 



The distance from the point (X|,x 2 ,x 3 ,x,i) to each of 
the four points was calculated and the species com- 
bination selected that had a minimum distance. 

There was still a need to handle minor species 
dropped by the fitting rule. That is, something had 
to be done to species 4 if the fitting rule showed 
that the point (xi,x 2> x 3 ,Xi) was closer to 
(1/3,1/3,1/3,0) than it was to (1/4,1/4,1/4,1/4). 
These species were added to a similar species 
based on a priority joining rule. Finally, for each 
plot we recorded the transformation that indicated 
the predominant species and what to do with all 
other species. For the purposes of this study only 
plots fit with one or two species were used and the 
others were skipped. 



ANALYSIS METHODS 



Pure Stands 

Once the plot information was summarized and 
stored in a retrievable manner, the pure stand 
analysis proceeded as follows: for a particular spe- 
cies, e.g., jack pine, the summarized data were 
used as boundary conditions in the nonlinear mul- 
tipoint boundary value problem (Leary 1970, 
Leary and Skog 1972). Solve 

?* [Potential diameter growth/tree] [NT] 1-e /BA 
At L J 

subject to the boundary conditions 

Y(t,) = k,, Y(t a ) = k 2f Y(t 3 ) = kg, where 
Y is the sum of tree diameters, 
NT is the number of trees, 



(2) 



ft is an unknown numerical constant (with units 
of basal area), 
ki,k2,k 3 are the boundary conditions, and 

AY = Y(t+l) - Y(t) 
At (t+1)- t 



(first forward difference). 



Other algebraic forms for the modifier were tested, 
this one most thoroughly 

(3) 



BA is the stand basal area (BA = .005454 I 



1 +-/BA 

but were found to be less desirable based on their 
pattern of residuals. 



Calculation of the potential diameter growth 
part of equation (1) requires use of mean stand 
crown ratio (Hahn and Leary 1979). Individual 
tree crown ratios on each plot had been recorded at 
the last measurement, but frequently the last 
measurement was not one of the times used to form 
the measurement set. Although most cases 
showed little difference between the ft values re- 
sulting from observed or predicted crown ratios, 
we used only observed crown ratios whenever 
there was a sufficient number of interval sets. This 
practice resulted in a large reduction in the sets of 
measurements available for use. For example, red 
pine had 603 sets of measurements spanning three 
measurements, but only 241 included the last 
measurement and an observed crown ratio. Be- 
cause we felt the observed crown ratio information 
was more precise, we used it unless doing so would 
have resulted in little or no information for a spe- 
cies. 

The ft values estimated for each measurement 
set judged to be pure were plotted against the 
mean stand diameter and their variation was 
studied. Based on these results, all measurement 
sets with 40 square feet or less of basal area stock- 
ing and with 50 percent (by number of trees) or 
more mortality and/or cutting between the begin- 
ning and final measurement were removed. Anal- 
ysis showed the ft values were independent of the 
number of trees, the site index, the interval length 
(5 to 29 years), and stand age to midpoint of the 
interval. There appeared to be a possible slight 
relationship between ft and average summer tem- 
perature over the interval, but data were insuffi- 
cient to warrant further study. The same held true 
for the effect of average summer rainfall on ft. 
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Some species showed a location (Minnesota, Wis- 
consin, Michigan) effect, but the limited range of 
mean stand diameters for many species made its 
study inconclusive. Some investigations were 
made on various weighting approaches, but the 
only one that offered any benefit was to weight 
proportional to plot size. 

Red pine, with 241 measurement sets, had the 
most data available to study any differences be- 
tween natural stands and plantations. The planta- 
tions tended to have smaller average diameters 
with the same basic pattern as the natural sets 
except for a scattering of higher values at about 8 
inches. A careful study was done to see if the cause 
of these higher values could be found, but it was 
unsuccessful. As a result, no separation as to mode 
of stand origin (plantation or natural) was made in 
the final analysis. 

In determining a function to fit to the /3's, we 
looked at those species with a fairly complete di- 
ameter range. Computing the average ft for 2-inch 
average diameter classes indicated a nearly con- 
stant function that decreases somewhat at the 
lower diameter values. From this relation the 
function 

fl-Ml-e 1 * 1 ) (4) 

was chosen. A few representative graphs are given 
in figure 3. 



Table 1 contains the coefficients bj and b 2 (equa- 
tion 4) for various species and the number of three- 
measurement sets on which each is based. Only 
490 measurement sets were usable after under- 
stocked plots, or those with high mortality or cut, 
were removed. Measurement sets were also lost 
because observed crown ratios were used, instead 
of predicted, whenever possible. 

Mixed Stands 

The mixed stand analysis proceeded similar to 
the pure stand analysis, but the model was ex- 
panded as follows: 

solve 



At 



(BA, 



1-e 



(6) 



(BA 2 -I- co a BA, 

-- 

At 

with boundary conditions 

Yl(t,) = k,, Yl(t) = k 12 Yl(ta) - k ia 
Y2(t,) = k Y2(t 2 ) - k,, Y2(t a ) - ku 




AVERAGE STAND DIAMETER (INCHESI 



Figure 3. Sample monomeokcular functions relating Rvalues 
to mean stand diameter. Conifer values tend to be considerably 
larger than hardwood values. The mathematical function is 
given in equation (4) and the numerical values are given in 
table 1. 



Table 1. -Parameter values for equation (4)' 



Species 

Jack pine 
Red pine 
White pine 
White spruce 
Balsam fir 
Black spruce 



Tamarack 



Northern white-cedar 
Hemlock 

Black ash 

Cottonwood 
Silver maple 
Red maple 
Elm 

Yellow birch 
Basswood 

Sugar maple 
White ash 

White oaks 

Select red oaks 
Ottier red oaks 
Hickory 



Bigtooth aspen 

Quaking aspen 
Paper birch 

Total 





Sets of 




b, bo 


measurement 


Comments 




intervals 




112 -.4638 


46 


Predicted crown ratio used 


121 -.4462 


241 


Observed crown ratio used 


117 -.1261 


77 


Predicted crown ratio used 


125 -.3652 


11 


Predicted crown ratio used 


86 -.4923 


15 


Predicted crown ratio used 




9 


data only for DBH 2"-4" (use 






white spruce coefficients) 




1 


Limited data (use jack pine 






coefficients) 


155 -.2230 


18 


Observed crown ratio used 




2 


Limited data (use white pine 






coefficients) 







Assigned sugar maple 






coefficients 


69 -.1434 


7 


Predicted crown ratio used 







Assigned sugar maple 






coefficients 




1 


Limited data (use coefficients 






for select led oaks) 




2 


Limited data (use coefficients 






for select red oaks) 


74 -.1902 


25 


Observed crown ratio used 







Assigned sugar maple 






coefficients 




4 


Limited data (use coefficients 






for red oaks) 


73 -.1303 


10 


Observed crown ratio used 







Assigned coefficients of 







select red oaks 






Assigned coefficients of 






select red oaks 







Assigned coefficients of 






quaking aspen 


117 -.1366 


21 


Predicted crown ratio used 







Assigned coefficients of 






quaking aspen 




490 





l For species in which the range ofaverage diameters was small, or only a few observations were available, the estimated beta values were 
used as a guide in selecting a stubstitute beta function. When no data were available, the species was assigned the beta function 
parameters for a similar species; for example, we had no pure bigtooth aspen or paper birch plots, so they were assigned the beta function 
coefficients for quaking aspen. 
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where 

Y1.Y2 are sums of diameters for species 1 and 2, 
respectively, 

PGi,PG z are potential diameter growth func- 
tions, 

NTi,NT 8 are numbers of trees, 

ADi,AD 2 are average diameters of trees in spe- 
cies 1 and 2, 

f i>f a are the monomolecular functions of average 
diameter (equation 4) for species 1 and 2, 

BA,,BA 2 are basal area (.005454 (Yj /NTt) 2 ) of 
species 1 and 2, and 

0)1,0)2 are the unknown interaction terms that 
tell how species 2 affects species 1 and how 
species 1 affects species 2, respectively. 

Only the plots judged to be two-species plots 
(figure 2) were used for this analysis. 

We tested the &> values resulting from using both 

a potential function containing the observed 

crown ratio and the predicted crown ratio. The 

esults were similar (actually slightly better for 

"edicted crown ratio), as judged by the deviations 

'observations from predictions. Since there were 

tore two-species plots available, we used the pre- 

icted crown ratio throughout the two-species 

inalysis. 

In the single species analysis we found that plots 
with less than 40 square feet of basal area had 
sharply decreased /3 values. The two-species anal- 
ysis showed that between 30 and 40 square feet 
they were nearly the same as above 40 square feet. 
Thus, we used 30 square feet as the cutoff value for 
understocked plots instead of 40 square feet to 
avoid losing additional data. The cutoff criterion 
for mortality, cut, or ingrowth was 50 percent of 
initial or surviving trees. Very few plots had to be 
discarded because of these criteria. 

From equation (5) it can be shown that if both o^ 
and tu 2 are exactly 1 and the ratio of average diam- 
eters is 1, we can say that species 2 affects species 1 
as species 1 affects itself, and species 1 affects 
species 2 as species 2 affects itself. In this analysis, 
the following relationships between coi and to 2 were 
thought to be most reasonable and were examined: 



In relationship [1] there is no fixed relation be- 
tween a) i and o) 2 , in [2] the sum of CD , and o> 2 is equal 
to 2, in [3] the product of w, and w a is 1, and in [4] the 
product of w, and w 2 is 2. 

Analysis of mixed plot data for each relationship 
in equation (6) showed that [11 is undesirable 
because w, and w 2 could vary widely with little 
difference in the goodness of fit of the solution of 
equation (5) to the observations. (This may have 
been caused by attempting to estimate four un- 
knowns, a*! and o). 2 and initial conditions for Yl and 
Y2 t based on only six observations of the system 
state.) Relationship [2] between w , and w 2 gave very 
poor results as judged by residual pattern and re- 
sidual magnitude. Relationships [31 and [4] both 
gave good results but Wi w 2 = 1 was slightly better 
and more logical. Thus, for each two-species mix- 
ture we had the constraint that w, co 2 = 1, For each 
plot of, for example, jack pine-red pine, we had an 
estimated value of i and w 2 . The mean values for 
each two-species mixture is given in table 2 along 
with the number of plots. Because we are dealing 
with the product of two averages, o>i and W K] the 
product does not necessarily equal 1, but is usually 
near 1. 



[1] w 2 and w, estimated separately 
[2] o 2 = 2- oj t 

[3] & 2 ~ Ifoi 
[4] w 2 - 2/wi. 



(6) 



Table 2. Interaction terms 
(5) 


oj, and 


w 2 for 


equation 


Measurement 


Species 


w, 


W 2 


sets 








available 


Jack pine-red pine 


2.757 


0483 


4 


Jack pine-quaking aspen 


.829 


1355 


3 


Red pine-white pine 


1.098 


1 192 


8 


Red pine-quaking aspen 


.455 


2,261 


2 


White pine-white oaks 


.776 


1.432 


2 


White pine-select red oaks 


3,500 


.286 


3 


White pine-quaking aspen 


1.783 


.662 


2 


Balsam fir-quaking aspen 


.581 


2.327 


16 


Black spruce-northern white-cedar 


.363 


3.066 


2 


Hemlock-sugar maple 


1.598 


.977 


5 


Red maple-yellow birch 


859 


1.205 


2 


Red maple-sugar maple 


1.571 


1.186 


4 


Red maple-select red oaks 


2.321 


.589 


5 


Red maple-paper birch 


.389 


2677 


2 


Elm-basswood 


.700 


2.010 


5 


Elm-sugar maple 


1.431 


.824 


11 


Elm-white oaks 


3.414 


293 


2 


Elm-select red oaks 


2,090 


.604 


7 


Elm-quaking aspen 


.880 


1.441 


2 


Basswood-sugar maple 


.976 


1.317 


24 


Sugar maple-white ash 


1.589 


.965 


9 


Sugar maple-select red oaks 


2.357 


.648 


22 


White ash-select red oaks 


.427 


2.346 


2 


White oaks-select red oaks 


.977 


1.601 


10 


White oaks-other red oaks 


.282 


3.546 


2 


Select red oaks-other red oaks 


1.120 


.893 


3 


Select red oaks-big tooth aspen 


.447 


2.497 


2 


Select red oaks-quaking aspen 


.499 


2.324 


2 
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The unit associated with /3 (equation 2) is square 
feet of basal area. Thus, the unit of bj in equation 
(4) is also square feet of basal area. The combina- 
tion of equations (2) and (4) makes this clear 



-b,(l-e b * DBH 



BA 



(7) 



Prom this relation we interpret the /3 function to 
give the basal area at which a stand with mean 
diameter, DBH, is growing at 63 percent of 
its potential. That is, when the exponent of e is -1, 
1 f-= 0.63. On this basis, one would expect toler- 
ant species to have higher function values than 
intolerant species. The softwoods in table 1 gener- 
ally follow this pattern although one would expect 
jack pine to have a lower bj value than balsam fir. 
The hardwood estimates deviate more from the 
theoretical ordering. One would expect quaking 
aspen to have the lowest b : value. Instead, it has 
the highest. 

Because of the way the model in equation 1 is 
assembled and the unknown numerical constants 
estimated, each subsequent step provides a check 
on previous steps. For example, the ft value esti- 
mates provide a check on the potential function in 
the following ways: (1) if the /3 values are ex- 
tremely large (recall the /3 unit is basal area) this 
means the potential function is predicting growth 
rates too low, and (2) if the /3 values exhibit a 
pattern that contradicts theory, such as quaking 
aspen having a higher b[ value than sugar maple, 
then the potential function for one of the species 
may need more work. Likewise, when estimating 
<Ui and ft> 2 , a check of the ft functions for species 1 
and 2, respectively, is made. For example, if two 
similar species have been observed growing to- 
gether and the observations are used to estimate 
MI and w 2 in the relation 

CUj Ct) 2 1, 

one would expect &>i to be very close to 1. A value 
significantly different from 1 would indicate that 
one of the /3 functions is not adequately character- 
izing pure stand growth of the species. 

Only pure stands and two-species mixtures were 
treated because the observational data simply will 
not support a thorough examination of three-spe- 
cies mixtures. Models for three- and four-species 
mixtures are built up from two-species models. For 



example, if we observe species 1 and 2 growing 
together and estimate <o 12 and oj 2 i for their interac- 
tion, if we observe species 2 and 3 growing together 
and estimate w 23 and <u 32 for their interaction, and 
likewise <a is and w 3 , for 1 and 3 growing together, 
then a three-species model is formed as follows; 



AYl 

At 



=[PG,J INT,] 



-/3,=f,(AD,) 



(BA, 



AD a /AD, + <a l3 BA 3 AD a /AD, 



1-e 



AY. 2 

At ' 



AY3 

At 



'[PG,]LNT 3 



'[PG 3 ][NT 3 ] 



(BA a + a,, BA, AD, /AD 2 + w 23 BA 3 AD 3 /AD 



1-e 



(8) 



-&=f 3 (AD 3 ) 



(BA 3 + u,, BA 2 ADz /AD 3 + (o 3 , BA, AD, /AD a 



1-e 



The explicit assumption here, for example, is that 
species 2 affects species 3 the same whether or not 
species 1 is present. 

For species mixtures where no information, or 
insufficient information, is available, a default 
value of 1 is given to the w values, pending acquir- 
ing additional information, 



CONCLUSION 

When the potential growth of a plot or stand of 
trees of a single species is characterized by the 
product of a potential function, number of trees 
and a modifier of the potential, the function 



satisfactorily reduces the potential growth to thai 
actually observed. The values of ft are species 
specific, being higher for conifers than foi 
hardwoods, and higher for tolerant species thar 
for intolerant species. For stands containing two' 
species mixtures, variation in the growth was sat 
isfactorily estimated by using two equations: 
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1-e 



1-e 



(BA, + a)! BA 2 



(9) 



(BA 2 + oj a BA, AD,/AD 2 ) 



We infer from this that mixtures of three or more 
species can be handled by expanding the number 
of equations. 
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THE PROBLEM 

Forest growth models may be categorized in 
many ways. An often-used dichotpmy is stand 
model vs. tree model. Munro (1974) presented a 
classification based on the use of inter-tree dis- 
tances and primary input parameters (variables). 
In this paper we have added a third dimension to 
Munro's: projection variables. We feel this is im- 
portant because one may choose not to project 
information at the level of detail potentially avail- 
able, but at an aggregated level. For example, 
even though a tree list containing d.b.h., species, 
stem maps, etc. is available as input, o'ne may 
choose to project the attributes of a class of trees. If, 
in going from the input form to the projection form, 
one aggregates tree attributes into class attri- 
butes, information is irretrievably lost unless a 
way can be found to disaggregate change in the 
projected variables. 

There are advantages to maintaining and up- 
dating a tree list through time. It allows more 
useful output because individual trees can be ag- 
gregated according to a criterion of interest for 
output. This is especially important when the out- 
put criterion is not the criterion used to aggregate 
input to the projected variables. For example, out- 
put information is often desired according to vari- 
ous merchantability standards, such as 4.5-inch 
d.b.h. (pulpwood) or 11.0-inch d.b.h. (saw log), 
There is, however, ample evidence that these size 
criteria have no biological significance useful for 
projection purposes. 

On the other hand, there are advantages to any 
approach that aggregates input information into a 
few classes for projection. Such approaches are 
more simple and economical, lending greater op- 
erational usefulness to the projection system. 

The purpose of this paper is to present a simple 
rule for disaggregating change in a projected vari- 



able i.e., of allocating change by a group of tret 
to the individual trees in the group. The followin 
discussion assumes that the basic input informa- 
tion is in the form of a tree list containing tree 
species and d.b.h., and that this information has 
been aggregated into a single class for projection 
purposes. 



APPROACH AND 
HYPOTHESIS 

Two fundamentally different ways to approach 
the allocation (disaggregation) problem may be 
termed absolute and relative. An example will 
show the nature of their differences. Assume we 
have a 12-inch d.b.h. sugar maple on our tree list, 
It is one of a group of sugar maple trees whose 
collective attribute (e.g., the sum of their diame- 
ters) is projected for 10 years. From values of the 
collective attributes at the initial and final time 
we calculate change. 

Implementing an absolute allocation rule to es- 
timate the diameter growth on the tree would nor- 
mally involve either consulting a table of mean 
diameter growth rates for sugar maple by size or 
size class, or solving a 10-year diameter growth 
equation for a 12-inch sugar maple tree. By apply- 
ing the same procedure to all other trees on the 
tree list and summing the projected growths, the 
sum could be compared with the amount of change 
to be allocated. If, for example, the predicted sum 
exceeds the amount to be allocated, one could use a 
reduction constant for all trees. 

Implementing a relative allocation rule would, 
on the other hand, consider the number of trees 
larger and smaller than the subject tree. Accord- 
ing to these numbers, the subject tree is allocated 
an amount of diameter growth, Because the social 
environment of a tree is an important determinant 
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of its growth, we have chosen the relative alloca- 
tion rule. 

The hypothesized relative allocation rule tested 
here is; growth of an individual tree in a stand is to 
the stand growth as growing stock of the individu- 
al tree is to the growing stock in the stand, or 

Adi . .JfSi w here the known quantities are 
Ad T GSr 

Ad T , GS h GS T (stand change given by the pro- 
jection equation, growing stock 
in the subject tree, growing stock 
total, respectively), and the un- 
known quantity is 
Ad, (diameter change of the subject tree). 

METHODS 

Only calibration data base sample plots judged 
to be "pure" (Leary and Holdaway 1979) were used 
for this analysis. The trees on each sample plot 
were sorted in ascending order by d.b.h. and then 
grouped into 10 cells of nearly equal frequency, 
corresponding to the 1-lOth percentile, the 11- 
20th percentile, etc., up to the 91-100th percentile. 
When the total number of trees, N, was not exactly 
divisible by 10, the extra trees were added, one per 
cell starting with the l-10th percentile cell, Thus, 
each cell contained approximately 10 percent of 
the trees on a plot with the smallest 10 percent in 
the first cell and the largest 10 percent in the last 
cell. 

The contents of each cell were analyzed for two 
attributes: (1) the growth of the trees in each cell 
relative to the total growth of all trees on the plot 
and (2) the growing stock represented by the trees 
in each cell relative to the total growing stock on 
the plot. Diameter growth per tree and sum of tree 
diameters were used as the respective variables. 



Proportion of stand growth due the 
computed as follows; 



Adbhu 



cell was 



k, 



Adbhy 



where k, is the number of trees in the ith cell, and 
Adbhy is the periodic d.b.h. growth for thejth tree 
in the ith cell. 

Proportion of growing stock in the ith cell was 
computed as follows: 

k. 



x, (2) 




where k, is the number of trees in the ith cell, and 
dbhy is the diameter breast height of theyth tree in 
the ith cell. 

These analyses gave two sets of 10 values for 
each plot having 10 or more trees: 

x, = the proportion of total growing stock in 

the ith 10-percentile cell, and 
y ( = the proportion or total growth due to trees 

in the ith 10-percentile cell. 

Proportion of growth and proportion of growing 
stock were computed for 10-percentile classes on 
each plot with 10 or more surviving trees in the 
time interval selected. Two intervals were selected 
on each plot: from the final measurement back 
approximately 10 years, and from the first mea- 
surement forward approximately 10 years. If over- 
Jap resulted, the plot was not used. A total of 1,599 
plot observations were available for analysis. 
Each of the 1,599 observations consisted of two 
sets of 10 proportions: growth and growing stock, 
as well as additional plot descriptors (mean plot 
basal area at midpoint of the interval, mean stand 
diameter, site index, property and plot number). 

Information on the 1,599 plots was sorted in 
order by species and mean stand basal area. Pre- 
liminary evaluations based on red pine plantation 
data showed a surprisingly linear relationship be- 
=yi (1) tween proportion of growth and proportion of 

growing stock. The same linear relation was evi- 
dent in the other plantation data available (jack 
pine and white spruce). Graphs of growth propor- 
tion against growing stock proportion for natural 
stands showed a marked difference from those for 
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plantations. This was especially evident for plots 
with greater mean d.b.h. values and/or lower 
stand basal areas. Figures la and Ib show the 
differences in trends as plot basal area changes. 
Figure Ib is for red pine plantations and figure la 
is for red pine natural stands. 

In both figures, there is a tendency for the rela- 
tion to rotate around or near the point with coordi- 
nates .10,. 10. This, of course, is to be expected since 
the_j>lot data were divided into 10-percentile 
classes. As basal area stocking decreases the fig- 
ures show three trends: the range in proportion of 
growing stock decrease, the range in proportion of 
growth decreases, and the slope of the relation 
changes from a very steep positive one at 179 
square feet in figure la to slightly negative at 66 
square feet in the same figure. Comparing these 
trends with those in figure Ib shows that the 
ranges also decrease, but the slope of the relation 
is still positive at 50 square feet. 



Several approaches were used to develop mathe- 
matical equations to characterize the types of 
changing relationships between proportion of 
growth and proportion of growing stock. The at- 
tempt was basically to estimate parameters for a 
general equation fit to each set of observations 
based on basal area classes. For example, a four- 
parameter model fit to observations for each basal 
area class would yield one set of parameters 
(a,b,c,d): 

GlUn = f(GS prop , a,b, C> d). (3) 

Given a set of parameters (a,b,c,d) for each basal 
area class, the attempt was to develop an expres- 
sion as follows: 

a = f(BA), where 
a is a parameter in equation 3, 
BA is mean stand basal area in square feet, and 

f is an unknown function. 
The general function, then, would be 

GR prop =flGS prop ,a=f 1 (BA),b-f 2 (BA), 

c=f a (BA),d=f 4 (BA)) (4) 
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Figure I. Proportion of growth compared with proportion of growing stock for 
(a) red pine natural stands and (b) plantations in the Lake States. 



where fj (i - 1,2,3,4) are functions relating basal 
area (BA) to parameters a,b,c, and d. 
This approach resulted in sets of parameter val- 
uesfor example, a set of six b values that varied 
so widely that no relation between the b values 
and basal area could be established. As a result, 
this approach had to be discarded. 

The approach finally used was to take the vari- 
ous graphic forms of relations between growth pro- 
portion and growing stock proportion shown in 
figure 1 and align a series of them on a light table 
so that they produced, as nearly as possible, a 
single smooth graph. To this single graph was 
added a new set of axes. The new coordinates of 17 
points were read for segments of the composite 
graph and used as observations to estimate the 
parameters in an equation to characterize the 
composite graph. The function used was: 
log(GRp rop ) = a + b(log GS mp ) + 



c(logGS prap ) 2 + d(logGS prop ) 3 



(6) 



where 

log is logarithm to the base 10, 

GRp rop is growth proportion, 

GSp r0 p is growing stock proportion, and 

a 0) b,c,d are unknown numerical constants. 

Estimated values of the numerical constants are 
a = -2.379 
b = -5.077 
c = -3.972 
d = -.2964 

t 2 = 0.986, Standard error = 0.0204). A graph of 
nis function is shown in figure 2, 

The strategy in putting together a composite 

relation is to have a mathematical function that 

possesses all the required types of shapes. The 

desired shape for a particular situation can be 

obtained simply by translating the axes so that the 

ibserved data coincide with a segment of the com- 

losite function graph. Because of the analogy with 

he giant tree concept (Grosenbaugh 1954), equa- 

ion (6) is called the "giant function". 



GR=k 1 +10 (a +b(log(GS+k 2 ))+ 
c(log(GS+k 2 )) 2 + 
d(log(GS+k 2 )) 3 ) 



(6) 



where a ,b,c,d have the numerical values given 
above, 

kj is the vertical translation of the axes, and 
k 2 is the horizontal translation of the axes. 



Observations of the growth proportions and 
growing stock proportions for basal area classes, 
such as shown in figure 1, were used to determine 
values of ki and k 2 for each basal area class for 
each species. Table 1 shows ki and k 2 values for 
eight red pine basal area classes for natural 
stands. The equation used to characterize the rela- 
tion between k : and mean basal area is a simple 
parabola: 

k^aj+a-jBA+agBA 2 (7) 

where 

ai,a 2 ,a 3 are unknown numerical constants, and 
the function for k 2 is 



where 

b 1( b 2 ,b 3 and b 4 are unknown numerical con- 
stants. 

Constants a i} i ~ 1,2,3 and b i( i 1,...4 were deter- 
mined by linear and nonlinear regression meth- 
ods, respectively, to give the following generalized 
form to the giant function: 

GR=k,+10 ( a + b ( lo g (GS+k 2 ))+ 
c(log(GS+k 2 ))Hd(log(GS+k 2 )) 3 ) (9) 

where a ,b,c,d are known constants for the gianl 
function," and k! and k 2 are computed from plot 
basal area using equations (7) and (8). 
In most cases it was necessary to examine the 
observations and assign numerical values to bj 
and b 3 in equation (8). Table 2 gives the coefficients 
of equations (7) and (8) for all species with suffi- 
cient data. 



Table 1. Vertical (k s ) and horizontal (k s ) transla- 
tion values for giant function to account 
for variation in relation between growth 
proportion and growing stock proportion 
for red pine natural stands (groups are 
by basal area classes) 



Basal area 


k, 


k 2 


R 2 


No. of 
plots 


Mean Range 


65.6 (20-69) 


-.1060 


.1346 


.42 


8 


79.5 (70-88) 


-.1111 


.1174 


.44 


21 


101.5 (90-109) 


-.1099 


.0961 


.53 


52 


119.1 (110-129) 


-.1043 


.0826 


.93 


151 


137.9 (130-149) 


-.1009 


.0734 


.97 


86 


156.0 (150-166) 


-.0985 


.0673 


.96 


39 


178.9 (170-189) 


-.0565 


.0315 


.99 


19 


197.6 (191-207) 


-.0689 


.0396 


.71 


5 
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0.32 



Figure 2. Graph of composite function in equation 5 showing location of 17 
observations from separate graphs. 



RESULTS AND DISCUSSION 

The functions to predict translation values, k t 
and k 2s show a considerable difference in the pat- 
tern of translations as basal area increases. Figure 
3 shows the patterns for Lake States species with 
sufficient plots judged to be pure. Species differ- 
ences relate primarily to the magnitude of ki and 
k 2) and the rate at which k, and k 2 approach zero. 
In general, tolerant species' values approach zero 
very slowly while intolerants rapidly approach 
aero. Quaking aspen appears to be an exception to 
this statement. 

Small differences between species' k t and k 2 
values can make a sizable difference in the rela- 



tion between growth proportion and growing stock 
proportion. For example, the kj and k 2 values for 
red pine natural stands and plantations (figs. 3g 
and 3h, respectively) are very similar, yet the pat- 
tern for 3g produces the relations in figure la and 
that for 3h produces the relations in figure Ib. 



The total sequence of operations in predicting 
growth proportion is given schematically in figure 
4. Four input quantities are required for a given 
tree and stand to determine the tree's d.b.h. 
growth; total stand growth (from the stand compo- 
nent projection equation), d.b.h. of the ith tree, 
sum of tree diameters on the plot, and the total 
stand basal area. 
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Figure 3. Translation values k, and k 2 for Lake States species. The solid lines 
denote the range of basal area for which data were available. The dotted lines 
are extrapolations based on equations (7) and (8) using coefficients in table 2. 
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Table 2. Coefficients of equations (7) and (8) used to predict vertical and horizontal translations of axes so 
that appropriate section of giant function is used to characterize growth proportion-growing stock 
proportion 



Species 






k, (vertical) 


h 2 (horizontal) 


i 




a 2 


a 3 


b. 


b 2 




b 3 


b 


4 


Jack pine (plantation) 


-.630 


E-01 


-.730 E-03 


.770 


E-05 


.10 00 


-.4000 


E-02 


.150 E+03 


.7400 


E 00 


Red pine (plantation) 


-.844 


E-01 


-.540 E-03 


.380 


E-05 


.15 EOO 


-.5585 


E-02 


.300 E+03 


.2139 


E+01 


(natural) 


-.643 


E-01 


-.968 E-03 


.500 


E-05 


.20 E 00 


-.9400 


E-02 


.300 E+03 


.4102 


b+U1 


White pine (natural) 


-.736 


E-01 


-.350 E-03 


,140 


E-05 


.10 EOO 


-.7300 


E-03 


,250 E+03 


.1759 


E 00 


White spruce (plantation) 


-.111 


EOO 


-.100 E-03 


.339 


E-05 


.10 E 00 


-.1799 


E-01 


.250 E+03 


.8264 


E+01 


Balsam fir (natural) 


-.871 


E-01 


-.242 E-03 


.283 


E-05 


.10 EOO 


-.1720 


E-02 


.250 E+03 


.3479 


E 00 


Black spruce (natural) 


-.100 


EOO 


.620 E-04 


.800 


E-07 


.10 E 00 


-.3190 


E-02 


.250 E+03 


.3724 


E 00 


No. white-cedar (natural) 


-.154 


EOO 


.583 E-03 


-.120 


E-05 


,10 EOO 


-.2018 


E-02 


.400 E+03 


.3724 


E 00 


Sugar maple (natural) 


-.840 


E-01 


-.240 E-03 


.240 


E-05 


.10 EOO 


-.6200 


E-02 


.250 E+03 


.8322 


E 00 


Quaking aspen (natural) 


-.819 


E-01 


-.247 E-03 


.248 


E-05 


.10 E 00 


-.3470 


E-02 


.250 E+03 


.7540 


E DO 


Other hardwoods (natural) 1 


-.394 


E-01 


-.196 E-02 


.150 


E-04 


.10 EOO 


-.2297 


E-01 


.200 E+03 


.6384 


E+01 



'Red maple, elm, basswood, yellow birch, white oak, select red oak 



f 

(TOTAL GROWTH) (PROPORTION GROWTH) =AD b h 



PROPORTION GROWING STOCK (OS) 




Figure 4, Sequence of steps in the estimation of diameter growth for an individ- 
ual tree in a stand. Steps a, b, and c establish the relation between growth 
proportion and growing stock proportion, Steps d, e, and f apply the relation to 
all trees in a stand. Total growth at step f comes from the projection equation. 
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In the projection procedure each tree is treated 
individually rather than being grouped into 10- 
percetitile classes. Considering a single tree on a 
plot of N trees is equivalent to dividing the plot 
into N 100-percentile classes. In this case the re- 
lation between proportion of growth and propor- 
tion of growing stock should have the same shape 
as the 10-percentile graph but should be centered 
at -^-^-instead of (.10, ,10). This can be accom- 
plished by transforming the per-tree growing 
stock value to a 10-percentile class value. The lat- 
ter value can be used as in figure 4 f step e, and then 
be transformed from a 10-percentile growth pro- 
portion value to a ((-jj-)lOO) percentile growth pro- 
portion value. The conversion relations are based 
on the following proportions: 



1 



:h>_._ _JJL_ nv Y Y 

x " ,10 01 X " X(I 



where 

x,, is the per tree growing stock proportion, 

N is the number of trees, and 

X is the 10-percentile growing 

stock proportion, 

The value X is used as in figure 4e and a 10- 
percentile growth proportion is computed, This 
value is transformed back to a (( ]r)10Qj percentile 
basis using the following relation: 



1 

_N_ 
.10 



, or Yo= 



Y 



where 

y is the per tree growth proportion, 
N is the number of trees, and 
Y is the 10-percentile growth proportion, 
Because of the rather complicated procedures 
used to estimate all the giant function parameters, 
we used observed growing stock proportion and 
observed basal area, and computed a predicted 
growth proportion. These values were compared 



with the observed growth proportion. It appeared 
that no systematic bias had been introduced into 
our procedures because the deviations of predicted 
from observed were as a rule normally distributed 
about zero for both the percentile class and stand 
basal area. This test does not, of course, constitute 
an independent test; rather, it simply checks for 
bias introduced by our procedures. 



The growth of an individual tree in a stand is a 
function of its position in the stand and of the stand 
growth. The initial hypothesis 

tree growth . , tree growing stock 
stand growth" * stand growing stock 
is overly simplistic. The growth proportion is re- 
lated to growing stock proportion in a nonlinear 
manner that changes with basal area stocking lev- 
els. Furthermore, the nonlinear relation differs 
between natural and plantation stands of the same 
species. 
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Knowing an individual tree's likelihood of dy- 
ing within the next year is essential to a tree- 
model forest growth simulator. A mechanism for 
selecting trees to be removed from a stand as a 
consequence of mortality is an integral feature of 
the North Central Forest Experiment Station's 
generalized forest growth projection system 
CLeary 1979). 

A tree's probability of dying depends upon its 
ability to cope with its environment. The tree it- 
self, its neighboring trees, and numerous ad- 
ditional influences including light, moisture, 
nutrients, insects, and fungi combine to form the 
unique niche for the individual tree. 

How might we assess the tree's response to the 
conditions of its neighborhood? One of the simplest 
indicators is the tree's vigor as reflected in its 
growth. Diameter growth rate is a readily obtain- 
able, direct measure of the tree's ability to with- 
stand the stresses of life. 

Diameter growth rate is the key tree character- 
istic in the species-specific mortality functions 
presented in this paper. We attempted to select a 
common functional form for all species, with 
unique parameters for each species. 



DATA ORIGIN 

Permanent growth plot data were available for 
the major forest tree species in Minnesota, Wiscon- 
sin, and Michigan (Christensen et al. 1979). Most 
individual tree records contained measurements 
for at least three points in time, permitting calcu- 
lation of an annual mortality rate for a given 
growth rate. Determining diameter growth rate 
required two diameter measurements; determin- 
ing the tree's status (live or dead) and subsequently, 
the mortality rate, required a third measurement. 



The time between measurements ranged from 1 
year to more than 10, with the most frequent inter- 
val being 5 years, 

Intervals of at least 2 years but no more than 6 
were chosen to establish a reliable measure of di- 
ameter growth rate. A comparable time interval (2 
to 6 years) following growth rate determination 
was chosen for recording tree status. 

These criteria reduced the size of the data base 
eligible for use in developing the mortality func- 
tion. Even so, more than 47,000 trees (table 1) had 
the required three measurements. The tree count 
ranged from 18,208 for red pine and 7,350 for sug- 
ar maple down to about 100 for white spruce, tam- 
arack, and bigtooth aspen. 

Furthermore, many of the trees had more than 
three measurements that met the time require- 
ments. However, only one set of measurement da- 
ta one diameter growth rate followed by a status 
observation was used from each tree when deriv- 
ing the mortality function. This set was chosen at 
random for each tree where two or more sets were 
available. 

The number of trees by species in the final 
mortality data base is presented for each diameter 
growth rate class (table 1) and for each diameter 
class (table 2). 



MORTALITY RATE 

Individual trees either lived, (mortality zero), 01 
died, (mortality 1.0). Diameter growth class mor- 
tality rates, however, depended upon the compos- 
ite of the individual rates of those trees forming 
the class. 

Growth class limits chosen were 0.0 inch (no 
growth), 0.02 inch, 0,04 inch, and so on up to 0.30 
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Table 1. Number of trees used 
classes 



in estimating mortality rates and parameters, by 5-year diameter growl 



Five-year diameter qrowth class (inches) 




.00- 
.09 


.10- 
.19 


.20- 
.29 


.30- 
.49 


.50- 
.69 


.70- 
.89 


.90- 
1.0 


Over 
I.O 


Tot 


Jack pine 


104 


165 


178 


492 


833 


763 


718 


528 


3,7* 


Red pine 


528 


605 


1,033 


3,987 


4,571 


3,663 


2,040 


1,781 


18 h 2C 


White pine 


49 


44 


55 


130 


139 


189 


170 


407 


1,1 


White spruce 


3 


2 


6 


10 


12 


25 


9 


39 


1C 


Balsam fir 


118 


207 


234 


562 


648 


431 


271 


331 


2,8C 


Black spruce 


12 


34 


59 


123 


96 


43 


21 


11 


3S 


Tamarack 


1 


5 


6 


16 


19 


8 


3 


3 


6 


N. white-cedar 





4 


11 


42 


53 


30 


5 


5 


15 


Hemlock 


17 


23 


34 


92 


68 


61 


39 


22 


35 


Black ash 


8 


10 


21 


45 


42 


23 


10 


5 


16 


Red mapfe 


36 


35 


52 


193 


246 


221 


166 


185 


1,13 


Elm 


168 


65 


141 


258 


241 


204 


234 


374 


1,68 


Yellow birch 


38 


32 


51 


143 


129 


100 


51 


53 


59 


Basswood 


103 


59 


110 


227 


262 


286 


328 


348 


1,72 


Sugar maple 


402 


216 


538 


1,145 


1,447 


1,347 


1,214 


1,041 


7,35 


White ash 


20 


22 


60 


129 


131 


182 


237 


319 


1,10 


White oak 


33 


29 


50 


117 


103 


66 


30 


18 


44 


Select red oaks 


19 


37 


50 


216 


323 


263 


250 


235 


1,39 


Other red oaks 


11 


4 


21 


28 


32 


20 


13 


11 


14 


Hickory 


22 


24 


60 


114 


67 


29 


13 


7 


33' 


Bigtooth aspen 





1 


4 


20 


20 


27 


20 


16 


10i 


Quaking aspen 


105 


213 


198 


398 


382 


389 


443 


791 


2,91! 


Paper birch 


64 


86 


91 


165 


167 


99 


78 


55 


80! 


Other hardwoods 


16 


20 


48 


54 


29 


36 


20 


26 


24! 


Noncommercial 


101 


22 


70 


77 


24 


10 


3 


2 


30! 


Total 


















.47,50: 



inch per year. These limits correspond to 5-year 
growth remeasurements (as mentioned earlier, 5- 
year remeasurements prevailed in the data base) 
of zero diameter growth, 0.1 inch, 0.2 inch, etc. 

The mortality rate (R), the annual mortality for 
a growth class, was calculated from the more read- 
ily determined annual survival rate (SR), with 
R := l"~fc)R. 

If we assume the SR for a species to be the same 
each year of the 2- to 6-year status measurement 
interval, while allowing a tree to die but once with 
its demise being independent of others of the spe- 
cies, we have 



= 2N,SR' 



(1) 



In (1) i is the length of the status observation 
interval, Si is the number of trees surviving the i th 



interval, and NI is the number of live trees at thf 
beginning of the i" 1 interval. 



These terms become more meaningful if we 
examine real status data, for example, the status 
of 528 "no growth" red pine trees: 



Number of 


Number of 


Number of 


live trees 


years (I) to 


live trees 


(Ni) at second 


third measurement 


(Si) at status 


measurement 


(status) 


measurement 





2 





16 


3 


13 


31 


4 


29 


417 


5 


347 


64 


6 


44 



48 



Table 2. Number of trees used in estimating mortality rates and parameters, by d.b.h. class 



Species 










n.M 


i. In Inches al 


second 


.measurement forlor to status observation) 


to 
1.5 


1.6- 
2.5 


2.6- 
3.5 


3.6- 
4.5 


4.6- 
5.5 


5.6- 
6.5 


6.6- 
7.5 


7.6- 
8.5 


8.6- 
9.5 


9.6- 
10.5 


10.6- 
11.5 


11.6- 
12.5 


12.6- 

13.5 


13,6- 

14.5 


14,6- 

15.S 


15.6- 
16.5 


16.6& 
larger Total 


Jack pine 


62 


426 


851 


482 


477 


509 


362 


199 


159 


103 


51 


42 


20 


16 


14 


4 


4 3,781 


Red plna 


12 


91 


234 


928 


1,937 


2,889 


2,851 


2,221 


1,811 


1,491 


1,124 


853 


628 


432 


320 


183 


20318,208 


While pine 








23 


72 


64 


70 


70 


76 


75 


81 


81 


78 


82 


79 


71 


26* 1,183 


While spruce 


2 


11 


12 


10 


7 


9 


19 


12 


10 


5 


3 


1 


1 


2 


1 




1 106 


Balsam fir 


92 


596 


631 


463 


298 


303 


207 


107 


61 


22 


12 


5 


4 


1 






2,802 


Black spruce 


3 


12 


35 


44 


68 


132 


80 


16 


6 


3 














399 


Tamarack 




3 




1 


9 


20 


11 


10 


3 


2 


2 












61 


N. white-cedar 






1 


1 


3 


36 


34 


21 


27 


15 


6 


4 




1 






1 150 


Hemlock 










30 


36 


44 


33 


35 


30 


33 


24 


23 


26 


16 


12 


14 356 


Black ash 








9 


17 


47 


33 


20 


7 


14 


8 


5 


1 


1 


2 




164 


Red maple 




7 


3 


19 


97 


237 


185 


156 


115 


92 


74 


50 


22 


28 


15 


10 


24 1 f 134 


Elm 




18 


69 


119 


156 


178 


180 


175 


172 


155 


136 


85 


65 


46 


35 


33 


63 1,685 


Yellow birch 




1 




6 


62 


140 


114 


61 


56 


54 


28 


12 


22 


15 


7 


8 


11 597 


Basswood 




20 


33 


68 


148 


225 


220 


195 


191 


170 


134 


72 


54 


33 


32 


27 


10t 1,723 


Sugar maple 


23 


252 


269 


743 


1,123 


1,324 


1 r 028 


705 


486 


336 


234 


187 


135 


113 


101 


86 


205 7,350 


White ash 




10 


20 


51 


86 


150 


139 


154 


142 


121 


94 


61 


29 


16 


11 


5 


11 1,100 


White oak 




2 






10 


5 


28 


39 


49 


60 


45 


49 


26 


29 


29 


17 


58 446 


S red oak 








1 


26 


72 


87 


90 


85 


109 


134 


120 


118 


119 


79 


83 


270 1,393 


Other red oaks 


3 


7 


10 


13 


12 


6 


6 


3 


6 


12 


17 


9 


6 


9 


6 


3 


12 140 


Hickory 




2 


4 


7 


11 


14 


50 


82 


70 


34 


29 


17 


6 


5 


3 




2 336 


Blgtooth aspen 










3 


20 


18 


15 


6 


12 


19 


11 


3 


1 






103 


Quaking aspen 




21 


98 


313 


304 


439 


513 


448 


300 


234 


147 


57 


29 


7 


4 


4 


1 2,919 


Papar birch 




169 


128 


72 


77 


132 


113 


55 


29 


14 


6 


4 


2 


1 


2 




1 805 


Other hardwoods 




10 


3 


11 


25 


30 


38 


25 


26 


1fl 


23 


18 


6 


7 


3 




6 249 


Noncommercial 




28 


39 


51 


85 


49 


30 


15 


7 


5 
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If the status observation follows the second di- 
ameter measurement by 2 or more years and the 
status interval is of common length for all trees 
(i.e.> one level of i), then (1) reduces to S = N SR YR , 
where YR corresponds to the length, in years, of 
the observation interval. This yields SR 
(S/N) 1/YR and can be solved by substituting the ob- 
served N, S, and YR in the function. Based on the 
second line in the tabulation (i = 3), the annual 
survival rate, SR, is (13/16) 1/3 or 0.933. Calculat- 
ing further, SR = 0.983 for i = 4, 0.964 for i = 5, 
and 0.939 for i = 6. 

When several different status interval lengths 
are encountered, as with i = 3, 4, 5, and 6, applica- 
tion of equation (1) will yield a composite SR for all 
intervals. Substituting observed values in S| = 
SNiSR 1 yields 433 - 16 SR 3 + 31 SR 4 + 417 SR 5 + 
64 SR 6 . Solving, we find SR to be 0.9611 and R to 
be 0.0389. 

We then have IS, - 433 and 2N|SR' = 433.05 
with the red pine data, both a check on our calcula- 
tions and the completion of the demonstration. In 
summary, we have; 



NSR 






2 





0.00 


16 


3 


13 


14.20 


31 


4 


29 


26.45 


417 


5 


347 


341.96 


64 


6 


44 


50.44 



Inspection of each species' mortality data re- 
vealed that decreasing mortality was associated 
with increasing diameter growth. When diameter 
growth exceeded 0.5 inch in 5 years, the mortality 
rate essentially leveled off without reaching zero. 
Tamarack and northern white-cedar were excep- 
tions; these species had consistently low mortality 
at all growth levels. 

There are four critical points where reliable 
mortality data are needed. Foremost is the mortal- 
ity rate for "no growth" trees. If the data at this 
point need to be bolstered by including trees show- 
ing slight growth, then the estimate of mortality 
at no growth loses some reliability. 
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A second critical data point is the mortality rate 
for the fast-growing trees. For each species the 
most ample data appeared here. 

At least two reliable data points between the 
nongrowers and the fast-growers are needed to 
establish the form of the mortality curve as mor- 
tality declines. 

Mortality rates for at least the four critical 
diameter growth rates were established for all spe- 
cies except other red oaks, tamarack and northern 
white-cedar. With tamarack and northern white- 
cedar the relationship of mortality to diameter 
growth was not established because too few trees 
were available, especially too few slow or no 
growth trees. Therefore, a uniform mortality rate 
was assigned to all tamarack and another to north- 
ern white-cedar. 



Three variations from this basic functioj 
were also evaluated for usefulness in eating 
mortality: 

R = [1+EXP (a+/3*DGR)]' 1 + S, 
R = [1+EXP (a+/3*DGinr, and 
R - [1+EXP (a-t-/>DGR y )]- 1 + 8. 
Still another model was considered, 
R = [1+EXP Ca+/3DGR v +ojDBH)]- 1 + 

to accommodate those species in which the eff< 
DGR depended upon tree size. 



0* 



* , 



MODEL SELECTION 

A mathematical expression for relating mortal- 
ity to growth rate was sought for each species. The 
logistic function Y = [l+EXP(a+2/3X)]' 1 has been 
used effectively for modeling tree mortality (for 
example, Hamilton and Edwards (1976) and Mon- 
serud (1976) and for modeling human health risk 
(for example, Truett et al (1967)). 

Writing the logistic in terms of our application 
we have 

with DGR being the annual diameter growth rate. 
Equivalent forms are R' =&' EXP (DGR) and 
R' = EXP (a+^DGR), with R' = (1 - R)/R. 



In these models, the parameter, a, 
mortality at the no-diameter growth level with 
very minor contribution of 8. DBH and its m 
plier, &>, provide an adjustment based on tree 
especially for trees having little or no growth, 
elusion of a background death rate, S, provides 
the occasional removal of a fast-growing tree, 
progression of the mortality estimates from 
growth trees to that of fast growers is controlled 
andy. Adding the exponent, % greatly 
the function's flexibility. 

Model (6) was selected for seven of the 
(table 3) and model (5) for the remainder. Esti- 
mates of the parameters under (5) are given for 
species (table 4) to facilitate parameter compe 
aon, even though (6) is the recommended model f or 
seven species. 

A preliminary screening of models (2) through 
(5), the first step in the selection process, was 
based on the ability of each to reproduce the mor- 
tality observed for each diameter-growth class. In- 
clusion of the background death rate (8), models (3) 
and (5), greatly improved the estimates for annual 
growth rates in excess of 0.1 inch where a low, but 
non-zero, mortality rate was observed. 

Table 3. Mortality function l [or selected Lake States tree species (estimates based on diameter and diameter" 
growthrute^ ., 

Parameters _,- 



Species 



ft 



OJ 



Jack pine 


0,3678 


60.178 


1.0727 


0.12220 


0.00502 


Red pine 


1.0843 


11.012 


.5887 


.48500 


.00020 


Balsam fir 


2.6976 


30.550 


.8996 


-.23069 


.00400 


Basswood 


-.5861 


101.870 


1.6493 


.68398 


.00489 


Sugar maple 


.8150 


66.937 


1.5934 


.75000 


.00435 


Other red oaks 


-1.1492 


9.109 


2,4552 


,35000 


.00000 


Quaking aspen 


-.0786 


213.630 


1.7126 


.25812 


.02276 



= diameter breast high. 
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The variable exponent for diameter growth rate 
Cy) yielded exponents less than 0.75 for three spe- 
cies, and more than 1.25 for 10 species (table 4) a 
substantial departure from 1. In these instances 
model (5) more realistically portrayed the decline 
in mortality for trees with no growth to trees with 
0.1-inch annual growth than did model (3). 

Finally, model (5) could well be superior to the 
others eliminated at this point, yet still be defi- 
cient. Numerous questions come to mind. Is it 
applicable for all sites? Is it applicable over the 
geographic range of the species, for all stand densi- 
ties, and for all diameters? And more. 

One major evaluation was completed; estimated 
tree mortality was compared to mortality observed 
within each diameter class. To accomplish this, a 
two-way table was constructed having 1.0-inch di- 
ameter classes as one dimension and 0.02-inch 
annual growth classes as the other. Observed and 



estimated tree deaths along with total trees were 
calculated for each cell as well as for each diameter 
class and each diameter growth class. 

Disagreement between observed and estimated 
tree deaths could be accepted, and was expected, 
within a diameter class. However, if these discrep- 
ancies formed a pattern whereby the estimated 
deaths were consistently too low for small trees 
and too high for large trees, or if the reverse pat- 
tern were true, then model (5) was deemed inade- 
quate. 

For six species, including red pine, balsam fir, 
basswood, sugar maple, other red oaks, and quak- 
ing aspen, model (5) overestimated mortality on 
the small trees and underestimated on the large. 
Jack pine revealed a modest reversal from this 
pattern. Model (6) adequately adjusted the esti- 
mated diameter class mortality patterns. Hence 
we used model (6) for these species, 



Table 4. Mortality functions 1 for Lake States tree species (estimates based on diameter growth rate) 

Parameters 



a 





7 


1.2178 


149.3669 


1.3699 


3.2185 


173086 


.6165 


2.3507 


33.6094 


1.0088 


2.8806 


30.7737 


1.2309 


2,0754 


29.5772 


.9541 


2.8774 


432.9522 


1.8576 


13.0000 2 


1.0000 


1.0000 


13.0000 2 


1.0000 


1.0000 


2.8233 


89 5603 


1.0966 


1.4875 


13.3384 


.5429 


2.2405 


53.5938 


1.0392 


1.9889 


44.5414 


1.0965 


2.6705 


63.1477 


1.3819 


2.2779 


25.3602 


.9027 


3.1831 


82.5642 


1.3605 


1.7821 


232.3075 


1.7267 


3.6946 


95.6369 


1.4517 


2.3092 


6,999.8220 


2.6615 


1.0044 


69.2301 


1.8048 


3.0626 


63.4946 


1.0498 


.8805 


56.5925 


1.2154 


1.2199 


693.5903 


2.0753 


2.1808 


11,6263 


.4750 


2.2339 


29.6290 


1.3427 


2.3581 


113.0644 


.9993 



Specjes^ 

Jack pine 

Red pine 

White pine 

White spruce 

Balsam fir 

Biack spruce 

Tamarack 

N. white-cedar 

Hemlock 

Black ash 

Red maple 

Elm 

Yellow birch 

Basswood 

Sugar maple 

White ash 

White oak 

S. red oak 

Other red oaks 

Hickory 

Bigtooth aspen 

Quaking aspen 

Paper birch 

Other hardwoods 

^commercial 

'Estimated annual probability of a tree dying Is (1/[1 +EXP(+/3-DGR T ) ]} + 
2 Mortallty is constant for all growth rates, relation to DGR not established, no 



0.00856 
.00042 
.00201 
.00396 
.00265 
.00235 
.01334 
.00402 
.00320 
.00018 
.00494 
00899 
.00382 
.00317 
.00405 
.00222 
,00188 
.00468 
.00000 
.00134 
,01288 
.02384 
.0002' 
0074 
.0192 



8 with DGR = diameter growth rats (Inches per year), 
slow-growing trees In sample. 
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VALIDATION 

To evaluate model performance the mortality 
models were applied to two data sets that were not 
utilized in model development. One set consisted 
of 266 pei-manent plots from the Cloquet Experi- 
mental Forest, Carlton County, Minnesota; the 
Forest is managed by the College of Forestry, Uni- 
versity of Minnesota. These plots provided test 
data primarily for aspen, birch, and northern coni- 
fers (table 5). An initial measurement and three 
re measurements, spanning a period of 17 years, 
were available. 

The second set consisted of 38 permanent plots 
from the Dukes Experimental Forest, Marquette 
County, Upper Michigan. The Dukes Forest is 
managed by the North Central Forest Experiment 
Station. These plots provided a small set of test 
data for northern hardwoods (table 5). An initial 
measurement and two remeasurements, spanning 
a period of 10 years, were available. 



One validation consisted of comparing the num- 
ber of trees, by species, expected to die during the 
status observation interval with those actually dy- 
ing in the interval. With this, the probability of 
death was determined for each tree of the species 
using that species' mortality model. As in the mod- 
el derivation process, diameter growth rate was 
based on the first two diameter measurements for 
a tree, and the status was observed at the third 
measurement. The count of trees observed to die 
was compared to expected deaths based on the- 
model (table 5). 



The largest discrepancy between actual tree 
deaths and model estimates was for jack pine. Ex- 
cessive deaths were yielded by the model, especial- 
ly for trees growing less than 0.08 inch per year, 
and for trees less than 8.6 inches in diameter. The 
tamarack and quaking aspen models also overesti- 
mated mortality, though the aspen differences 
were not large. 



Table 5. Observed and estimated tree mortality on validation plots by tree species 
_ _ CLOQUET PLOTS 



Initial measurement through 
second remeasurement 1 



First remeasurement through 
third remeasurement 2 



Jpecies 


Trees in 


Mortality 


Trees in 


Mortality 




test 


Observed 


Estimated 


test 


Observed 


Estimated 


Jack pine 


1,674 


111 


211 


1,597 


207 


278 


Red pine 


899 


3 


5 


1,032 


6 


17 


White pine 


170 


11 


6 


185 


20 


9 


White spruce 


69 


6 


4 


104 





6 


Balsam fir 


513 


73 


65 


613 


130 


94 


Black spruce 


778 


54 


55 


832 


118 


81 


Tamarack 


300 


6 


19 


318 


16 


29 


N. white-cedar 


354 


7 


7 


394 


16 


11 


Ash 


61 


1 


5 


76 


13 


8 


Red maple 


75 





3 


93 


8 


5 


Bigtooth aspen 


53 


8 


10 


50 


10 


14 


Quaking aspen 


611 


89 


103 


653 


166 


174 


Paper birch 


800 


8 


12 


889 


39 


24 








DUKES 


PLOTS 3 






Elm 


71 


1 


4 





_ 


__ 


Basswood 


67 


1 


2 











Sugar maple 


400 


5 


6 












'Diameter growth rate based on 5 years' growth, status observed 5 years later. 
2 Diatneter growth rate based on 5 years' growth, status observed 7 years later. 
3 Diameter growth rate based on 7 years' growth, status observed 3 years later. 
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Numerous species, including white spruce, 
northern white-cedar, bigtooth aspen, elm, bass- 
wood, and sugar maple, showed good agreement 
between actual and estimated mortality. With the 
exception of 10 predicted deaths among the 148 
1-inch trees (none were found), the red pine esti- 
matfisjwere good. Ash and red maple showed-varia- 
tioftl^rresurfiEfbetween observation periods,"but on 
balance yielded good results. 

Mortality was underestimated for white pine, 
primarily with the 6-inch trees; for balsam fir and 
black spruce, primarily with the 5-, 6-, and 7-inch 
trees; and for paper birch. 

The second validation consisted of simulating 
mortality for each plot, paralleling the use of mor- 
tality models in stand growth simulation. Based 
on random selection in proportion to each tree's 
likelihood of death, each tree was assigned either a 
mortality or a survivor status. A count of these 
mortality trees provided an estimate of plot mor- 
tality through simulation for comparison with ob- 
served plot mortality. 

To arrive at the probability associated with 
each mortality count, each tree on a plot was given 
100 assignments of 1 (mortality tree) or (survivor 
tree). This was done through generating 100 ran- 
dom numbers between and 1 for each tree and 
comparing each number with the tree's likelihood 



of death. If the random number did not exceed this 
likelihood, the tree was assigned a status of 1, 
death. Otherwise the tree was scored 0, a survivor 
tree. 

The first scores were added for all trees on the 
plot to arrive at one mortality _count for the plot. 
This process was repeated for each of the 99 re- 
maining scores. A distribution of these 100 counts 
was made for each plot, providing a relative fre- 
quency or probability for each mortality count. 

Consider the results from Cloquet plot 82 (table 
6). This plot had 30 trees with two diameter mea- 
surements. When status was recorded 5 years lat- 
er, one dead tree was found. Expected mortality 
using the mortality functions for the species pres- 
ent was 0.94. On this plot the probability of select- 
ing exactly one tree to die in the 5-year interval 
based on 100 applications of the selection process 
was 0.30, i.e., 30 of the 100 simulation runs pro- 
duced exactly one dead tree. The probability of 
selecting either zero, one, or two trees, a count 
within one of the observed mortality, was 0,91. 

Ten Cloquet plots demonstrate the discrepan- 
cies that might be encountered in stand growth 
simulation (table 6). These plots were selected be- 
cause of their close agreement between observed 
and expected mortality. Still, the highest propor- 
tion of hits, simulation results equal to observed 
mortality, was the 0.47 (47 of the 100 simulation 



Table 6. Relative frequencies (probabilities) of mortality tree counts for 10 Cloquet plots determined by 
random selections in proportion to each tree's probability of dying in 5 years (basis: 100 sets of tree 

counts per plot) _____ 

Mortality trees 



Mortality tree counts 



Plol Trees Observ. Exp. 1 



1 



8 



10 11 12 13 



82 30 



-Number- 



0.94 .43 .30 .18 .( 



.01 ----- - 



68 


19 


1 


1.01 


.36 


.42 


.17 


.05 























237 


30 


1 


1.00 


.34 


.47 


.16 


.03 





. 





. 











4 


36 


2 


1.55 


.26 


.28 


,27 


.15 


.02 


.01 


.01 














398 


18 


2 


2.43 


.01 


.18 


.40 


.18 


.16 


.05 


.02 











, 


147 


53 


3 


3.48 


.02 


.10 


.12 


.25 


.17 


.15 


.10 


.07 


.02 








111 


33 


4 


4.63 


.01 


.01 


.11 


.13 


.21 


.24 


.17 


.08 


.02 


09 

. Ut. 


, 


285 


29 


5 


4.65 





.02 


,11 


.17 


.23 


.21 


.13 


.06 


.05 


.02 





177 


36 


5 


4.92 





.01 


.05 


.12 


.21 


.23 


.21 


.14 


.02 


n-j 
.u i 





204 


45 


7 


7,88 


. 





.01 





,02 


.13 


.18 


.17 


.16 


.17 .10 .03 


.02 .01 



'Expected mortality is the sum of the tree probabilities of dying over 5 years. 
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runs produced one dead tree each) found for plot 
237. Plot 177, with observed mortality at five trees 
and expected at 4.92, showed only 23 hits in 100 
applications. 

Clearly, the mortality counts from stand growth 
simulation will often differ from the actual, even 
when the model is appropriate. However, the dem- 
onstration material (table 6) indicates we can ex- 
pect to be within one tree of the true count most of 
the time. 

The Cloquet and Dukes plots having at most one 
mortality tree were summarized on the basis of 
conformity between each plot's observed mortality 
and the results from mortality simulation (table 
7). On four plots (Cloquet plots having 11 or more 
trees, one with no mortality and three with one 
mortality tree each) the probability of getting the 
correct mortality count through simulation was 
zero. At the other extreme, 26 of the Cloquet plots 
having 10 or fewer trees and no mortality showed 
probability of 0.91 to 1.0 of getting the correct 
mortality through simulation. 



Except for plots having 10 or fewer trees and no 
mortality, most plots did not have a high probabili- 
ty of hitting the observed mortality through mor- 
tality simulation. However, most plots met the 
relaxed criterion of being within one tree of the 
correct number. 

Only 31 of the 87 Cloquet plots with two or more 
mortality trees had a probability exceeding 0.20 of 
producing the correct mortality count through 
simulation (table 8). Furthermore, the criterion of 
being within one tree of the correct count wag not 
met consistently by plots having three or more 
mortality trees. 

Thirteen plots containing 11 or more trees each 
had zero probability of producing the correct mor- 
tality count through simulation. Of these, 11 had 
mortality counts in excess of the observed mortali- 
ty, and two were less than the observed mortality. 
Jack pine was the prime contributor to nine of the 
plots having excessive counts. 



Table 7. Probability of conformity by plot between mortality tree counts determined by random selections in 
proportion to each tree's probability of dying (simulation) and observed counts for each plot having 
at most one mortality tree (basis; 100 sets of simulation tree counts per plot) 

____ (Number of plots) ________ ___ 

Dukes plots 2 



Cloquet plots 



Plots: 10 or 


fewer trees Plots: 11 or 


more trees 


Probability 


Mortality (O) 1 


Mortality (1) Mortality (0) 


Mortality (1) 


Mortality (0) 






Within 




Within 


Within 




Within 




Within 




HIH 


1 tree 


Hit 


1 tree Hit 


1 tree 


Hit 


1 tree 


Hit 


1 tree 


.00 








1 




3 


2 






.01-.10 






3 


13 


2 


5 


2 


1 




.11-.20 






2 


8 




5 


1 






.21-.30 






2 


6 


5 


7 


2 


1 




.31-.40 






1 


5 


5 


11 - 




2 




.41-.50 


3 




2 


5 


5 


5 


4 


2 


1 


.51-.60 


1 






5 


4 




2 


4 




.61-70 


2 






6 


5 




3 


7 


1 


.71-.80 


6 


1 




10 


5 




2 


5 


A 


.81-.90 


16 


3 


1 


4 


11 




5 


8 


3 


.91-1.0 


26 


50 




11 8 


29 




13 


2 


23 


Total plots 


64 


54 


11 


11 71 


71 


3fi 


36 


3? 


32 



'One plot with two mortality trees not included. 

'The simulation count for the plot Is the same as the observed mortality. 



Mortality (1) 
Within 
Jit 1 

1 

2 
2 
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Table 8, Probability of conformity by plot between mortality tree counts determined by random selections in 
proportion to each tree's probability of dying (simulation) and observed counts for each Cloquetplot 
having two or more mortality trees (basis: 100 sets of simulation tree counts per plot) 



(Number of plots) 


Plots: 10 or fewer trees 


Plots: 11 or more trees 


Mortality (2 or 3) 1 


Mortality (2) 


Mortality (3) 


Mortality (4 or 


More) 


Probability 




Within Within 




Within 


Within 




Within 


Within 




Within 


Within 




Hit 2 


1 tree 2 trees 


Hit 


1 tree 


2 trees 


Hit 


1 tree 


2 trees 


Hit 


1 tree 


2 trees 


.00 


2 




2 






4 


2 


1 


3 


1 


1 


.01-. 10 


5 


1 


3 


2 


1 


6 


3 


2 


15 


5 




.11-.20 


1 


1 


9 


1 


1 


3 


1 




3 


5 


3 


.21-.30 


3 


3 


10 


1 




5 


2 


1 


7 


5 


4 


.31-.40 




2 


5 


2 


1 


1 


3 


1 




4 


1 


.41-. 50 








1 






2 


2 




1 


4 


.51-.60 




1 




3 


1 




2 


1 




5 


3 


.61-70 








12 


1 




3 


2 




2 


5 


71-.80 




2 




5 






1 


3 






1 


.81-.90 








2 


5 






4 






4 


.91-1.0 




1 11 






19 






2 






2 


Total plots 


11 


11 11 


29 


29 


29 


19 


19 


19 


28 


28 


28 



'Either two or three mortality trees observed on these plots. 

2 The simulation count for the plot Is the same as the observed mortality 
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Gary J. Brand, Research Forester 



Previous papers have described a forest growth 
projection system capable of "growing" and "kill- 
ing" trees, Because it does not select trees for re- 
moval, the system applies only to unmanaged 
stands. To accurately project managed stands a 
procedure is needed that will remove trees (simu- 
late cutting) according to prescribed guidelines. 
This paper describes the procedure (algorithm) 
used to simulate tree removals recommended in 
management guides, and indicates some possible 
uses for the projection system with the addition of 
the algorithm. FORTRAN subroutines called by 
the projection system computer program (Hahn 
and Brand 1979) implement the tree removal 
algorithm, 



dominant and codominant stems, basal area, basa! 
area of dominants and codominants, or some other 
method. The Society of American Foresters (1964) 
favors using the number of dominant and codomi- 
nant stems. The USDA Forest Service (1974) 
states that cover type should be based on species 
most fully occupying the area. The latter point 
indicates that tree size as well as number is impor- 
tant. Since crown class cannot be projected by the 
model, a classification based on crown position is 
impractical. To incorporate size into the proce- 
dure, basal area rather than total stems was used 
to determine cover type. An algorithm used by 
North Central Forest Experiment Station Forest 
Survey 2 is the basis for determining the type. 



DESCRIPTION OF THE 

SILVICULTURAL 

ALGORITHM 

The algorithm has three major parts. Because 
ifferent species produce different products, ma- 
-ure at different ages, and grow under different 
conditions, the guides used by forest managers are 
species dependent, and are therefore developed for 
specific cover types. In the projection system a tree 
list characterizes the portion of the forest of inter- 
est (Leary 1979). The first part of the algorithm 
determines, on the basis of this tree list, the cover 
type and thus, the most appropriate guide to use. 
The selected guide, along with plot values, deter- 
mine a silvicultural action to be carried out on the 
tree list. The final part of the algorithm performs 
the suggested action. 

Cover type may be determined by the species 
with the greatest number of stems, number of 



work was conducted under a cooperative 
agreement between the North Central Forest Exper- 
iment Station and the College of Forestry, Universi- 
ty of Minnesota, St. Paul, Minnesota, 



SELECTING THE ACTION 



In most cases, guides prepared by North Central 
Forest Experiment Station silviculturiats were 
used to develop the algorithms. Each guide was 
translated to a diagram that systematically evalu- 
ates stand characteristics to reach a final recom- 
mendation. The sources for the diagrams are as 
follows: 



Cover type 
Jack pine 
Red pine 
White pine 

Spruce-fir 

Black spruce 
Northern white-cedar 
Mixed swamp conifers 
Lowland hardwoods 



Source 

Benzie (1977a) 

Benzie (1977b) 

USDA Forest Service 

(1958) 

USDA Forest Service 

(1967a) 

Johnston (1977a) 

Johnston (1977b) 

Johnston (1977b) 

Godman (personal 

communication) 



2 Hahn, J. T. Personal communication, 1976. 



fifi 



Northern hardwoods 



Oak-hickory 
Aspen -paper birch 



Tubbs (1977), Godman 
(personal 
communication) 
Sander (1977) 
Perala (personal 
communication) USDA 
forest Service (1967b) 
Understocked stands of ^Benzie (personal 
above cover types communication) 



Cover Typing 

The purpose of the typing algorithm is to deter- 
mine the most appropriate cover type and there- 
fore the correct guide for evaluating the tree list. 
Most tree lists will have a cover type assigned to 
them that was either determined in the field or by 
computer analysis. In either case, there are advan- 
tages to determining cover type directly from the 
tree list using a computer algorithm. One is that it 
will ensure that identical tree lists have the same 
cover type regardless of their source. It also allows 
the cover type to change with projected composi- 
tion changes. 

These guides have several advantages over 
other possible guides. They enjoy rather wide- 
spread use in the Lake States and thus it is more 
likely that region-wide projections will be 
accurate. They exist for the major timber types in 
the Lake States, so there is some consistency. And 
finally, they do not recommend actions so 
intensive that the model cannot respond to them. 
For example, fertilization could not be adequately 
represented by the model because coefficients are 
not available to predict growth responses to 
fertilization. 

Since management guides are interpreted by 
foresters in the field, professional judgment can 
influence the choice of a final action. To use com- 
puter methods in determining an action requires 
decisions using a tree list that does not contain all 
the information available in the field. As a result, 
some simplifying assumptions had to be made in 
constructing algorithms for the various types. 
These assumptions were checked with silvicultur- 
ists and appear to be reasonable. 

Various actions are suggested by the guides 
such as "do nothing," "clearcut", "thin", "shelter- 
wood cut", and "remove overtopping trees". The 
recommended actions are based on site quality, 



stand characteristics, and cover type. For exam- 
ple, clearcutting is prescribed in a jack pine stand 
if the site index is 60 and stand age is 65 years, 
Clearcutting in an aspen stand, however, is pre- 
scribed if the stand age is 36 years. 

Decision diagrams containing these criteria 
were produced for each type. Stand and site char- 
acteristics are used to guide decisions in the dia- 
grams, which are constructed along the lines of a 
taxonomic key (fig. 1). For example, assume a jack 
pine tree list has the following stand characteris- 
tics: average stand diameter 6 inches; basal 
area 100 square feet/acre; site index 50 feet; 
stand age 35 years; number of trees 500/acre; 
and BAMIN 55.6 square feet/acre (BAMIN is the 
amount of basal area needed to achieve full stock- 
ing.). The dashed line in figure 1 indicates how a 
final recommendation would be reached. Since av- 
erage stand diameter is 5 inches, the bottom 
branch is followed. Again the bottom branch is 
taken since the basal area is greater than 55.6 
square feet per acre. The site index and age 
indicate that the top branches are now 
appropriate, leading to the final recommendation, 
"do nothing". 

A disadvantage of this approach is that it lacks 
flexibility only this decision diagram can be used 
for jack pine. To allow more flexibility, the critical 
values are read as input variables and can there- 
fore be changed. Figure 2 shows the more general 
form of the decision diagram for jack pine. Now, by 
selecting JP(7) as 50 years and the rest equal to 
zero, jack pine can be managed on a 50-year rota- 
tion with no intermediate cutting. The guides can 
be manipulated to test the effect of different maxi- 
mum and residual basal areas, different rotation 
ages, or less intensive management. 



Marking 

After an action has been selected according to 
the guide, it must be carried out. The marking part 
of the algorithm selects trees for removal accord- 
ing to tree characteristics. In the simple case of 
clearcutting, all live trees are "marked" as cut. For 
more complicated actions, such as thinning, char- 
acteristics must be selected that allow ranking of 
trees for desirability of removal. Low-preference 
trees are marked for cutting until the required 
residual stand is achieved. 
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Figure I. Specific jack pine decision diagram. 
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The various removal procedures used to carry 
out specified actions are listed below: 



Removal procedure 

1. rough and rotten 

2. small trees, 
unspecified species 

3. large trees, 
unspecified species 

4. small trees, 
specified species 

5. large trees, 
specified species 

6. shelterwood cut 

7. reduce number of trees 

8. all overtopping 

9. all trees of a certain 
size and species 



Input needed 
size class to remove 
size class to remove 
four species to leave 
size class to remove 
four species to leave 
size class to remove 
four species to remove 
size class to remove 
four species to remove 
size class to remove 
four species to favor 
number of trees 
to leave 

lowest d.b.h. of trees 
to remove 
size class to remove 
four species to remove 



The removal procedures have been selected to 
accomplish all the final actions recommended by 
the management guides. For instance, jack pine 
stands less than 60 years old with 120 square feet 
of basal area per acre and a site index greater than 
GO should be thinned from below to 80 square feet 
per acre. The first, second, and fourth procedures 
will accomplish this action. Rough and rotten trees 
are removed first, and if the basal area is still too 
high, the smallest tree of an undesirable species is 
removed. When all undesirable species are re- 
moved and the residual stand has still not been 
achieved, the smallest trees of desirable species 
are removed. Marking is discontinued when the 
residual basal area is achieved. 



As with action selection, it is desirable to have 
the ability to change the marking rules. To permit 
this flexibility, different hierarchies of removal 
procedures may be selected. This option permits 
up to five of the removal procedures listed above to 
be selected for each recommended action. Return- 
ing to the jack pine example by selecting the 
first, third, and fifth removal procedures the stand 
can be thinned in such a way that the largest trees 
are removed first. 



RELATIONSHIP TO THE REST 

OF THE PROJECTION 

SYSTEM 

Because effective use of the management guide 
algorithm depends on the rest of the projection 
system, it is important to present a simplified de- 
scription of how the pieces interrelate. A tree list 
(representing a stand) is brought into the system 
and its cover type and silvicultural action are de- 
termined as previously described. The action is 
carried out and the appropriate trees are coded as 
cut. One year's growth is determined for the re- 
maining trees and the mortality function applied. 
Trees are coded as dead if selected by the mortality 
function. The revised tree list is cycled through the 
growth and mortality functions to determine an- 
nual changes until the desired number of years 
between cutting cycles has occurred. Cover type, 
silvicultural action, and trees to mark are again 
determined using the current tree list. This proce- 
dure continues until the tree list has been pro- 
jected for the desired number of years. (If another 
stand (tree list) is to be projected, the new tree list 
is brought into the system and the entire cycle is 
repeated.) 



POSSIBLE USES 

Incorporation of the management algorithms 
into the projection system makes possible several 
additional uses of the system. Survey information 
in the form of tree lists, is available on a statewide 
basis. Applying the projection system to these tret 
lists would produce timber yields expected frorr 
the state if all lands were managed according t< 
the guides. Comparison with present yields wouk 
indicate the intensity of management as well ai 
opportunities for increasing timber yields. 

Another use for the system would be testing 
various management guides and marking rules 
Actual or contrived tree lists could be projected 
using different guides, and the yields then com- 
pared. This would provide a much faster means of 
comparing alternatives and help identify those 
most promising for field testing. By experimenting 
with marking rules, the effect on timber yields of 
leaving rough and rotten trees as den trees or 
favoring different species or size classes could be 
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determined. Comparing these yields with those 
obtained from the original marking rules would be 
a rneaure of the "cost" of the new marking rule. 
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and Gary J, Brand, Research Forester 



The various parts of the generalized forest 
growth projection system have been described 
elsewhere (Buchman 1979, ChristenseneiaJ. 1979, 
Leary, et aL 1979, Hahn and Leary 1979, Leary 
and Holdaway 1979, Brand 1979). TRES (Timber 
Resource Evaluation System) is the computer pro- 
gram that combines the parts into a package capa- 
ble of simulating the growth of trees entered into 
the system. This paper is not intended to be a 
detailed user's guide, but rather a general descrip- 
tion of the input needed, how the parts are put 
together, and the output produced. 

INPUT 

The input consists of a parameter deck and a 
tree list. A tree list consists of tree species, d.b.h., 
crown ratio, and quality class for a group of trees. 
Also required are site index, stand age, and plot 
size. 

Parameter Deck 

The parameter deck specifies two groups of in- 
put values: user-selected options, and numerical 
constants needed by each portion of the growth 
processor. 

Numerical constants are used for: 

(1) assigning species group codes to individual 
trees, 

(2) crown ratio function coefficients, 

(3) silvicultural treatment specifications, 

(4) volume equation coefficients, 
(6} mortality equation coefficients, 

(6) potential growth function coefficients, 

(7) growth modifier coefficients, 

(8) allocation rule coefficients, and 

(9) a species grouping priority matrix (used for 
combining species when more than three 
species occur on a plot). 



These values are determined from calibrating the 
model on a specific data base of permanent (re- 
measured) plots. 

Projection related options that may be selected 
by the user are: 

(1) length of projection cycle (years), 

(2) number of projection cycles, and 

(3) type of growth projection desired (valida- 
tion, projection, or update). 

The number of years the tree list is projected is 
determined by the number and length of the cy- 
cles. The projection is limited to a maximum of 100 
years, but it may be in the form of 100 cycles of 1 
year each or one cycle of 100 years. The length of 
the cycle determines how often summaries are 
made and the tree list examined for management 
opportunities. 

Three types of growth projections can be se- 
lected: validation, projection, and update. A val- 
idation run uses a tree list with measurements 
recorded for several different years. The tree list is 
projected from the initial measurement to the final 
measurement. Projected and actual measure- 
ments are compared to determine the accuracy of 
the model. A projection run starts with initial con- 
ditions and grows the trees for a specified length of 
time to estimate future conditions. Finally, an up- 
date run starts with initial conditions from the 
past and grows the trees to the present. Using 
known product removals, the trees that produced 
those products are removed from the tree list. This 
update section of the program is currently being 
developed. An example of an echo of the input for a 
validation run appears in Appendix A. 

Tree List 

Regardless of whether the tree list comes from 
sample plots or is derived (for example, from a 
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stand table), the tree list can be in one of three 
forms. The first is a list in the format described by 
Christensen et al. (1979). The second is tree infor- 
mation in the unit record format used by the 
USDA Forest Service Renewable Resources Eval- 
uation (forest survey) projects. And finally, a tree 
list in a computer-readable format specified by a 
user may be selected, in which case the user is 
responsible for supplying the necessary computer 
code for input and output routines. 



PROGRAM ORGANIZATION 

Besides providing input/output capability, the 
program controls when the different parts of the 
system are used according to the options selected, 
Table 1 gives a listing of all the subroutine and 
function names, their calling routines, and a brief 
statement of their purpose. The program has seven 
parts: 

(1) the main calling program, 

(2) a function, VALUE, where all mathemati- 
cal equations are located, 

(3) input subroutines, 

(4) a subroutine for aggregating the input tree 
list, 

(5) stand treatment subroutines, 

(6) stand growth subroutines, and 

(7) output subroutines. 

Figures 1 to 7 contain flow charts for each of the 
seven parts. Figure 1 shows the structure of the 
main calling routine and figure 2 the structure of 
the function VALUE. 

Figure 3 shows the set of subroutines for reading 
the input parameters and the tree list. The struc- 
ture of subroutine READ is shown in figure 3a. It 
reads input parameters and echo-prints all values. 
Figure 3b shows the subroutine structure for read- 
ing tree lists in master format (Christensen et al. 
1979); figures 3c and 3d show subroutine structure 
for reading tree lists in Renewable Resources 
Evaluation (forest survey) unit record format. The 
possibility of a user-written subroutine is indi- 
cated by figure 3e, 

When all input (parameters and tree list) has 
been read, the tree list is partitioned into species 
groups and size classes by subroutine FIT (fig. 4), 
called from the main program TRES (fig. 1). 




Figure 1. Main program calling TRES. 
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Figure 2. Function VALUE. 
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Table I. Subsystem TRES (Timber Resource Evaluation System) of program FREP (Forest Resources 
Evaluation Program) 



'rogtam TRES 

Main calling routine Sets up and initializes arrays, performs 
housekeeping functions, organizes calling of subroutines 
Function VALUE 

Called from RMASTR, RSURTP, GROW, ISUM. OUTPTI, and OUTPUT2 
Computes values for crown ratio, potential growth, modifier, 
mortality, volume, and growth allocation functions. 
Subroutine READ 
Called from TRES 

Reads user-selected options and numerical constants and 
echo-prints all values. 
Subroutine RMASTR 
Called from TRES 

Reads data from binary flies in MASTER format into arrays TRS and ITRS 
Initializes array parameters 
Subroutine RSURTP 
Called from TRES, 

Reads data from files in forest survey unit record format Into arrays TRS and ITRS, 
Initializes array parameters. 
Subroutine OPENRD 
Called from RSURTP 

Unblocks forest survey unit record file and supplies one record at a time to RSURTP 
Has entry point READIT. 
Subroutine RUSER 
Called from TRES 

Supplied by user to read user dat file Into arrays TRS and ITRS. 
Subroutine FIT 

Called from TRES, SCREEN and REMOVE. 
Classifies trees in arrays TRS-ITRS into projection groups 
Subroutine SCREEN 
Called from TRES. 

Screens plots for possible treatment and/or sets up treatment schedule 
Subroutine CQVTYPE 
Called from SCREEN 

Determines cover type using live trees in tree list. 
Subroutine TREAT 
Called from SCREEN 

Performs silvlcullural treatments prescribed In SCREEN. 
Subroutine GROW 
Called from TRES. 

Grows trees, selects trees formortallty, allocates growth to individual trees, updates TRS and ITRS array. 
Subroutine REMOVE 
Called from GROW. 
Removes trees to meet specified annual removals for updates or trees recorded as having been 

removed by remeasurements. 
Subroutine ISUM 

Called from GROW, OUTPTI, TRES 

Summarizes trees Into classification groups and computes group totals. 
Subroutine OUTPTI 
Called from TRES. 

Outputs summaries for specified points in time. 
Subroutine TFILE 
Called from TRES. 

Outputs updated forest survey unit record file. 
Subroutine OPENWR 
Called from TFILE. 
Writes blocked unit record tape 
Has entry points WRITIT and CLOSWR. 
Subroutine OUTPTS 
Called from TRES. 
Outputs file for validation tests. 
Subroutine QUTUR 
Called from TRES. 
Supplied by user to output updated tree list to user specifications. 



3 Input 
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Figure 3 .Subroutines (a} READ, (b) RMASTR, (c) RSURTP, (d) OPENRD, and (e) RUSER. 
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Figure 4. Subroutine FIT. 



Selecting the option to manage the stand as it is 
being projected, following the management guides 
of Brand (1979), causes subroutines SCREEN, 
COVTYPE and TREAT to be called. Figure 6a 
shows the flow chart of SCREEN, from which 
COVTYPE and TREAT are called. COVTYPE 
(fig, 5b), which is also called at the end of a pro- 
jection by OUTPUT2, determines the cover type; 
TREAT (fig. 5c) performs the silvicultural treat- 
ments selected in SCREEN. Trees selected for re- 
moval are designated as cut, and the live trees are 
again divided into species and diameter groups by 
a recall of subroutine FIT. This allows the changes 
in species composition caused by cutting to modify 
the numerical constants used to project the plot- 

The tree list is projected 1 year at a time for the 
number of years in a cycle and for the number o* 
cycles in the run by subroutine GROW, REMOVE 
and ISUM (fig. 6). A plot is then projected through 
the end of projection period; it cannot be revisited 
until all plots have been processed. The croitfj 
ratio function (Holdaway et al. 1979), potent***' 
growth function (Hahn and Leary 1979), and mo^' 
ifier function (Leary and Holdaway 1979) are use* 
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to determine total growth of each diameter/species 
projection group (fig. 6a). Individual tree growth is 
then determined from the growth allocation func- 
tion (Leary et al. 1979). Mortality trees are deter- 
mined vising the mortality function (Buchman 
1979) and paeudo- random number generator. 



Four different output types are possible: 

(1) "Yield table" output containing such varia- 
bles as basal area and number of mortality 
trees, surviving trees, basal area and sum of 
diameters on a per plot and per acre basis, 
and cut tree volumes on a per acre basis. 
This output is given by subroutine OUTPT1 



(fig 7a) and is available at the end of each 
cycle as well as at the beginning and end of 
the projection period. Appendix B shows an 
example of this form of output. 

(2) A new Renewable Resources Evaluation 
(forest survey) unit record file (subroutine 
TFILE and OPENED, figs. 7b and 7c respec- 
tively). 

(3) A file for validating the projection system. 
This is produced by subroutine OUTPT2 
(fig. 7d) and has meaning only if two or more 
measurements have been made on the plot. 
The file produced by OUTPT2 is ordinarily 
stored on disk and analyzed by a separate 
user- written program. 

(4) A user-designed output file (subroutine 
OUTUR, fig. 7c). 






Figure 5,-^ubroutines (a) SCREEN, (b) COVTYPE, and (c) TEEAT. 
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Figure 6. Subroutines (a) GROW, (b) REMOVE, and (c) ISUM. 
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TEST OF UPDATE CHANGES USING CLOflUET FILE 
tt**tttt_ ,PJ,OI_SUHMARIES. AT THE END OF 1959 PLOT IDENT. 1 ***** 



NUMBER OF MORTALITY TREES PER ACRE 


SPECIES 
GROUP ALL 


DIAMETER 
TOTAL 


CLASS 


JACK PINK! 
PAP. BIRCH 








BLK. SPRUCE o ' 

TOTAL 


6 





DIAMETER SUMS PER ACRE 


SPECIES 
GROUP ALL 


DIAMEIEK 
TOTAL 


CLASS 


JACK PINE 37?. 4 
PAP, BIRCH 125.3 


372.4 
135.3 




ELK. SPRUCE 
TOTAL 497.7 




497.7 




BASAL AREA OF MORTALITY TREES PER ACRE 


SPECIES 
GROUP ALL 


DIAMETER 
TOTAL 


CLASS 


JACK PINE 
PAP. BIRCH 








BLK. SPRUCE 
TOTAL 








BASAL AREA ON PLOT 






SPECIES 
GROUP ALL 


DIAMETER 
TOTAL 


CLASS 


JACK PINE 3.1 
PAP. BIRCH .6 


3.1 

.6 




BLK, SPRUCE 
TOTAL 3.7 



3.7 




NUMBER OF TREES PER ACRE 






sptcity ' 

GROUP ALL 


UJAHtltH 
TOTAL 


CLAbb 


JACK PINE 35.0' 
PAP. BIRCH 21.0 


35.0 
21.0 




BLK. SPRUCE 
TOTAL 56.0 


111 o 1 

56.0 




BASAL AREA PER ACRE 






SPECIES 

GROUP ALL 


1 DIAMETER 

TOTAL 


CLASS' 


JACK PINE 2?.0 
PAP. BIRCH 4.1 


22. 'o 

4.1 




BLK. SPRUCE 
TOTAL ' 26.0 



26.0 




NUMBER OF TREES ON PLOT 






SPECIES 
GROUP ALL 


DIAMETER 
TOTAL 


CLASS. 


JACK PINE 5.0 
PAP, BIRCH 3.0 


5,0 
3,0 




BLK, SPRUCE 
TOTAL B.O 



8.0 




PRESENT STAND 


VOLUMES 




SPP.GP. CUFT/A. 
JACK PINE 389.50 


CORDS/A, BDFT/A. 
4,93 1866,35 


PAP. BIRCH 68.91 
TOTAL 458.41 


.87 
5.80 1860.35 
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Rolfe A. Leary and Jerold T. Hahn, 

Principal Mensurationists, 
and Roland G. Buchman, Principal Biometrician 



The purposes of this paper are (1) to report on 
tests of the design for assembling projection sys- 
tem parts, and (2) to report on tests of projection 
system reliability in various use domains. The use 
domains are many. For example, the forest type 
domain covers all Lake States types; we are inter- 
ested in the projection system reliability for each 
type. Another domain is that of stand origin mode 
plantation or natural. A third relates to geo- 
graphical regions within the Lake States; does the 
system work better in one State than in others? 
For which types in which State? A fourth and 
particularly important domain is time. What is 
the time decay rate for projection system reliabili- 
ty? Other domains include stocking levels, sites, 
ages, species, sample plot types (fixed or variable 
radius) and mean stand diameters. Each domain 
or combination of domains warrants investiga- 
tion. In this paper we report on item 1 above and 
the domains of (1) cover type, (2) geographical 
area, (3) time, and (4) sample plot type. 

TESTING MODELS OF 
SYSTEM DYNAMICS 

Mathematicians recognize two possible sources 
of error in predicting the future state of a dynamic 
system, such as a forest stand: 

(1) errors in estimating the initial state (initial 
conditions) of the system, and 

(2) errors in the mathematical model that char- 
acterizes system dynamics. 

Errors in predicted future states may be caused by 
either of these sources or both may act in concert. 
Source (2) may be further broken down into: 

(a) the mathematical equations are inappropri- 
ate for describing the process, 

(b) the manner of combining the equations is 



incorrect, and 

(c) the numerical constants obtained from cali- 
bration are not those governing the system. 

An example of error type 1 for the projection 
system is when a person does not have an observed 
tree list but only a dot tally of trees on the plot by t 
say, 4-inch diameter classes. 

An example of error type 2a is if the potential of 
trees to grow in diameter at breast height is not 
satisfactorily characterized by our potential func- 
tion (Hahn and Leary 1979). Error type 2b is exem- 
plified by a finding that the realized stand growth 
cannot be characterized by a product of a potential 
and modifier functions. Error type 2c is obvious; 
the numerical constants in use do not characterize 
the specific process under test. 

Considerable effort has gone into developing 
theoretically sound mathematical equations for 
describing the process of forest growth. Hence er- 
ror type 2a will not be evaluated in this paper. 
Such evaluations belong more properly in the sep- 
arate papers decribing each equation. We will, 
however, evaluate the projection system to assess 
the magnitude of type 2b and 2c errors. We are not 
concerned with assessing type 1 errors this is an 
inventory problem. We have prevented type 1 er- 
rors by using only permanent growth plot informa- 
tion. To assess type 2b errors we have used the 
permanent growth plot information available to 
calibrate the mathematical equations (Christen- 
sen et al 1979). If the projection system performs 
well when checked against the calibration data 
base, it may be tested against an independent data 
base (i.e., one not used in calibration). By doing so, 
and by comparing predicted future states with ob- 
served future states (starting from a known past 
condition) we are able to assess the magnitude of 
type 2c errors. The validation data base for this 
comes from five sources: 
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(1) Cloquet Experimental Forest of the College 
of Forestry, 1 University of Minnesota; 

(2) Chequamegon National Forest, Wisconsin; 

(3) Nicolet National Forest, Wisconsin; 

(4) A second-growth northern hardwood study 
(NH-2) on an experimental forest near Mar- 
quette, Michigan; and 

(5) Hiawatha National Forest, Upper Michi- 
gan. 

Figure 1 shows the geographical locations of 
sources. The fifth source, the Hiawatha National 
Forest, has an eastern and a western part. The 



l We thank Dietmar Rose and the staff of the Cloquet 
forestry Center for making past measurement data 
available and for assistance in making a measure- 
ment in 1976. 



distribution of plots by age class and site index 
class is shown in table 1. 



To test both the method of assembling system 
parts and the applicability of the system to other 
Lake States forest conditions, the validation op- 
tion was selected in the computer program that 
implements the projection system (Hahii and 
Brand 1979). Selecting this option causes the pro- 
gram to: 

(1) read up to four sets of observed diameter 
measurements for each tree on a plot and 
store them for comparison with sets of pre- 
dicted diameters, 




Figure 1, Geographical location of the five 
sources of permanent growth plot data used for 
independent tests of the projection system. 



80 






40 : 50 : 60 ; ?0 "0 



NORTHERN WHITE-CEDAR 



VHITE SPRUCE-BALSAM FIR 



3 1 

NORTHERN HAftDWOQD 



21 -SO 

31-40 

41-50 

51-60 

61-70 

71-80 

91-100 

100* 



LOWLAND HARDWOOD 



41 -SO 

51-60 

61-70 

71-80 

81-90 

91-100 

100* 

Total 



0-20 
21-30 
31-40 
41-50 
51-60 
61-70 
71-80 
81-90 
91-100 

loot 

Total 



6-20 
21-30 
31-10 
41-50 

51-60 

61-70 

71-80 

91-100 

lOOt 

Total 



5T4TJ 
11-50 
51-60 
Total 




(2) use the first set of diameter measurements 
to start the projection, 

(3) remove a tree from the projection tree list 
when it has been observed to have been cut 

(4) add a tree to the list of trees being projected 
when it is observed to be an ingrowth tree, 

and , , , 

(5) write an output file that contains observed 
and predicted diameter and status for each 
tree present on a plot at the first measure- 
ment. This is done for up to three points in 
time. 

Ingrowth trees are not written on the output file 
(step 5) for comparing observations and predic- 
tions. They are included in the projection process, 
however, because in some cases ingrowth trees 
comprise a major part of the stand by the time ot 
the final observation. Ignoring them could cause 
serious overestimation. 



Analysis of the information produced at step 5 
proceeds in two directions: (1) analysis for accura- 
cy, and (2) analysis for precision. 

The accuracy analysis is based on all the plots 
that occur in a forest type. We co rnpare the ^pre- 
dted and observed mean stand diameter trees 
per acre, and basal area per acre in a forest type by 
forming the ratio 




Clearly, ratioa consistently greater than one m- 
dicat overprediction (positive bias), and those 
condstently less than one indicate underpredic- 

observed and predicted are expressed as integers. 
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The precision analysis is needed to assess the 
repeatability of the estimated ratios for each type. 
Here we compute for each plot in a forest type the 
ratio of predicted to observed variable (mean stand 
diameter, trees per acre, basal area per acre), the 
mean ratio over all plots, and a standard deviation 
of the ratios. 

Conceivably the projection system could, at 10 
years, have a positive bias of exactly 6 percent for 
every plot in a forest type. This would not be an 
accurate predictor but would be precise. Lack of 
accuracy can, of course, be corrected. Lack of preci- 
sion cannot be easily corrected. Rather, it indi- 
cates a general inability of the underlying model 
to capture the essence of the processes underlying 
stand dynamics. 

To test the method of assembling the parts of the 
system (error type 2b) we used selected forest prop- 
erties from the calibration data base, and did an 
accuracy analysis only. To test the applicability of 
the numerical constants to other forests in the 
Lake States region (error type 2c) we used the five 
properties in the validation data base (fig. 1), and 
did both an accuracy and precision analysis, In 
both cases, we based our analysis on all plots in a 
property regardless of their stocking level; plots 
having as few as three or four trees were used as 
well as fully stocked plots. Clearly, failure to satis- 
factorily pass the test on method of combining 
system parts makes the other test largely academ- 
ic. Finally, it should be noted that the plots are 
typed by the projection program (Hahn and Brand 
1979), and that trees of several species may be 
present on a plot, not just trees of the species the 
type name bears. 



TEST RESULTS 
AND DISCUSSION 

This section is broken into two parts; the first 
gives the results of testing the manner of combin- 
ing system parts, and the second tests the applica- 
bility of the complete system, including values of 
numerical constants, to other forests of the Lake 
States. 

Combining the System Parts 

The results of tests are shown in figure 2. The 
model seems to predict mean stand diameter bet- 



ter than either number of trees or basal area. This 
is not surprising, since an error in predicting num- 
ber of trees even if diameter is predicted exact- 
lywill usually cause an error in predicting total 
stand basal area. 

Even so, the three ratios for most species are 
gratifyingly close to one at most points in time. 
Several points can be made based on the contents 
of figure 2: 

(1) Predictions for Minnesota jack pine are de- 
cidedly better than those for Lower Michi- 
gan when the evaluation is based on mean 
stand diameter and total stand basal area, 
However, the model estimates number of 
trees better for Lower Michigan stands. 

(2) Predictions for Minnesota and Wisconsin 
red pine plantations and natural stands are 
quite good, although nearly all show ratios 
somewhat greater than one. 

(3) White pine mean stand diameter is pre- 
dicted fairly well. The mean d.b.h. for the 
Minnesota source is underestimated, while 
that for the Wisconsin source is somewhat 
overestimated. Urenholdt trees per acre af- 
ter 22 years is off some (23 percent), but 
trees per acre at the other locations are quite 
good. 

(4) Only one small plantation was available to 
calibrate the system for white spruce. Even 
after 25 years the predictions are good, 

(5) Predicted values for the balsam fir type are 
less than the observed values. The average 
stand diameter prediction gets progres- 
sively worse for the first 18 years, but then 
improves over the remaining 7 yeai-s of 
observations. 

(6) Both average diameter and number of trees 
are underestimated for the black spruce and 
northern white-cedar types. Combined, 
these greatly underestimate basal area. In- 
terestingly, most of the prediction error 
occurs in the first 10- to 13-year period, 
Thereafter the prediction essentially paral- 
lels the observation. 

(7) Average diameter for test plots in the north- 
ern hardwood type is, for the most part, 
within 4 percent of the observed. An ex- 
ception is the Dundee timber harvest forest, 
which overestimates average diameter by 
12 percent after 27 years. Perhaps the sever- 
est test of the projection system occurred 
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JACK PINE TYPE 



10 



YEARS SINCE INITIAL MEASUREMENT 
15 20 25 30 

i_J t I. ..I I I I I J i~J I l_J I l~4_ 



35 
.. . '- 



1 NC-57 (P/22) AD 

TR 
BA 

2 NC-40 (P/25) AD 

TR 
BA 



.90 
.95 
.94 



.96 

.95 

I 

1.05 

.97 

1.12 



.99 
1.09 
1.04 



1.13 

.99 

1.29 



1.14 
1.00 
1.31 




RED PINE TYPE 



1 NC-84 (N/11) AD 

TR 
BA 

2 NC-155 <P/54) AD 

TR 
BA 



1.02 
1.00 
1.34 



1.03 
.39 

106 



1.03 1.03 

.99 99 

1 OS 1.06 



1.01 
1.00 
1.02 



1.01 

1.01 
1.03 




RED PINE TYPE 



3 STAR (P/B) AD 

LAKE TR 

BA 

4 NC-431 (P/27) AD 

TR 
BA 



10 
.til 



15 



20 



30 



35 



1.03 
1.04 
1,07 



.97 

1.00 

.92 



.99 

1.01 
.95 



1.03 
I.OB 
1.09 



.95 

1.01 

.90 



1,03 
1.09 
1.09 




WHITE PINE TYPE 



1 NC-56 (N/27)AD 

TR 
BA 

2 IN/10) AD 
URENHOLDT TR 

BA 

3WASAUKEE (N/18JAD 
TR 
BA 



,66 
1,01 

,91 



1.00 
.92 
.97 



.96 

1.02 

.91 



1.00 

1.04 

.99 



.96 
1.03 

.94 



1.03 
.77 

.as 
1.00 

1.03 
.98 



1.01 

1.04 

.99 




WHITE SPRUCE TYPE 



1 NC-18 (P/12) AD ,98 

TR 1,00 
BA .93 



BALSAM FIR TYPE 



2 NC-26 (N/2Q) AD 
TR 
BA 



99 

1.04 
.94 



.94 

.99 

1.01 



.92 

.99 
.98 



1.00 

1.10 
.94 



.97 
.91 
.95 



gure 2. Test results of accuracy analysis of 
method of combining system parts. Each forest 
property is identified by name or number, planta- 
tion (P) or natural (N) stand and number of plots 
used in the comparison. Criteria for comparison 
ITQ mean stand diameter for the type (AD), mean 



trees per acre for the type (TH), and mean stand 
basal area (BA) per acre for the type. Maps show 
approximate location of each property in Lake 
States. Ratios are computed as 
predicted- mean stand attribute for aforesMype 
observed mean stand attribute for a forest type 
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BLACK SPRUCE TYPE 

1 NC-59 (N/12) AD 

TR 
BA 

NORTHERN WHITE 
CEDAR TYPE 

2 NC-32 (N/19) AD 

TR 
BA 

NORTHERN HARDWOOD TYPE 

1 CAIRNS (N/16) AD 

TR 
BA 

2 (N/12) AD 
NETTLETON TR 

BA 

3 NH-25 (N/105) AD 

TR 
BA 



.93 
.95 
.84 



.95 
.93 
.87 



.94 
.94 
.SB 



,95 
.98 
.67 



.96 
.94 
-B7 



.98 
.89 
.98 



.99 

.93 

1.11 



.97 

.95 

1.08 



.96 

1.07 

.89 



.95 

1.09 

.68 



.98 

1.00 

.96 



,98 

1.00 

.95 



.98 

1.04 

.99 



.96 

1,05 

.OS 



4 NH-138M (N/10) AD 1.00 

TR 10 
BA 1.06 

5 STONE'S (N/19) AD 

WOODS TR 

BA 

6 DUNDEE (N/12) AD 1.04 

TR .95 

BA 1.06 



1 

.99 
1.10 
1.05 



1.01 
1.07 
1.06 

.99 
1.00 
1.05 



.96 
1.14 
1.19 



1.06 

.95 

1.07 



1.12 
1.02 
1.23 



OAK-HICKORY TYPE 



Jt 



15 



20 



25 



30 



1 STONE'S (N/65) AD 

WOODS TR 

BA 

2 DUNDEE (N/18J AD 

TR 
BA 



.99 
.98 
.99 



.96 
1.02 
1.04 



1.02 
.91 
.97 



1.01 
.88 
.94 



1.03 
.83 
.96 



35 

..>.... -I.-, 




ASPEN TYPE - 

2 NC-96 (N/28) 



1 NC-26 (N/25) 



AD .96 
TR 1.01 
BA .87 

AD 
TR 
BA 



.97 
.97 
.85 



.96 

1.01 

.82 



.92 
.82 
.86 



.91 
.78 
.85 



.95 
.73 
.88 



with the Nettleton plots. The average stand 
diameter of the first measurement was be- 
tween 1 and 2 inches. Even though pro- 
jections were made for 36 years, the 
estimated mean stand diameter was within 
4 percent of the observed. 

(8) Predictions for the oak-hickory type, as 
- found in southern Wisconsin, are reason- 
ably good, although there is a tendency to 
overestimate mortality on the Dundee tim- 
ber harvest forest. 

(9) The Minnesota aspen type is reasonably 
well predicted for both average diameter 
and number of trees. The fact that basal area 
is underestimated so much is probably due 
to allocating too much of the diameter 



growth to small trees. The Wisconsin aspen 
comes from the same study as the balsam fir 
test. The 25 plots under test here probably 
had a large component of balsam fir. As with 
the 20 plots in the balsam fir type, mortality 
is overestimated and the average diameter 
is underestimated. 

Tests on Independent Data Sets 
Accuracy analysis 

Five different data sets form the basis for these 
tests (fig. 1). The 266 Cloquet plots have been 
measured four times. The ratios of predicted to 
actual stand attributes is shown in table 2 for 5, 10, 
and 17 years after the initial measurement. Ratios 



Table 2. Summary of Cloquet Experimental Forest test showing ratios of projected to actual mean plot 
diameter, trees per acre, and basal area per acre for the type totals converted to a per-acre basis 
(Cloquet plots are all VT acre, fixed-radius) 



Forest 
type 


No. of 




Years after 


Initial measurement 


plots 


5 10 


17 5 


10 


17 


5 


10 


17 



(ratio of predicted 
to actual mean plot 
diameter) 



(ratio of predicted 

to actual number of 

trees per acre) 



(ratio of predicted 

to actual basal area 

per acre) 



Jack pine 


43 


1,02 


1.03 


1.04 


1.02 


1.03 


1.06 


1.05 


1.09 


1.13 


Red pine 


53 


.99 


.39 


.98 


1.00 


.99 


.99 


.99 


.97 


.96 


White pine 


6 


1,00 


.99 


.96 


1.04 


1.11 


1.20 


1.02 


1.03 


1.01 


W.spruce-fir 


37 


1.02 


1.03 


1.04 


1.08 


1.12 


1.26 


1.12 


1.19 


1.37 


Black spruce 


29 


1.01 


1.03 


1.05 


1.05 


1.09 


1.18 


1.07 


1.14 


1.32 


N,w.~cedar 


15 


1,01 


1.02 


1.04 


.99 


1.01 


1.03 


1.01 


1.06 


1.12 


Tamarack 


15 


1.02 


1.03 


1.06 


.97 


.97 


.97 


1.01 


1.04 


1.09 


Lowland hdwd. 


2 


1.00 


1.01 


1.01 


1.13 


1.07 


1.39 


1.09 


1.00 


1,26 


N, hardwood 


2 


1.01 


1.05 


.92 


1,00 


1.00 


1.33 


1.02 


1.10 


1.20 


Aspen 


32 


1.02 


1.05 


1.07 


.93 


.93 


1.00 


.97 


.99 


1.11 


Paper birch 


26 


1.02 


1.03 


1.04 


1.03 


1.03 


1.01 


1.06 


1.09 


1.10 



Table 3. Summary ofNH-2 test projection, for the northern hardwood type (NH-2 plots are Vie acre, fixed- 
radius) 



Number 
of 
plots 


Years after 


initial measurement 


7 


8 


11 : 


7 


8 


11 : 


7 


8 


11 



37 
62 



(ratio of predicted 
to actual mean plot 

diameter) 
1.00 ,99 



(ratio of predicted 

to actual number of 

trees per acre) 

.99 .99 



(ratio of predicted 

to actual basal area 

per acre) 

1.02 1,00 
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Table ^.-Summary of 10-year national forest test projections (all comparisons are based on repeated measure- 
ments of variable radius 10-point cluster plots) 



Cheouameaon 


Nicolet 


Hiawatha 


Forest 
type 


No. 
of 
plots 


D.b.h. 


Trees 


Basal 
area 


No. 
of 
plots 


D.b.h. 


Trees 


Basal 
area 


No. 
Of 
plots 


D.b.h. 


Trees 


Basal 
area 


Jack pine 


7 


1.04 


1.02 


1.09 


5 


1.06 


1.18 


1.23 


20 


1.09 


1.03 


1.16 


Red pine 


6 


.99 


1.11 


1.00 


6 


1.02 


.85 


.96 


12 


.99 


.94 


.92 


White pine 


3 


.96 


1.09 


.93 


2 


.93 


.84 


.87 


3 


.97 


.99 


1.00 


W. spruce-fir 


15 


1.05 


1.09 


1.10 


14 


1.03 


.91 


1.05 


8 


1.02 


1.fa3 


1.16 


Black spruce 


2 


1.13 


.65 


1.02 


3 


1.04 


1.15 


1.26 


5 


1.04 


1.14 


1.16 


N.w. -cedar 


6 


1.00 


1.00 


1.02 


12 


1.02 


.91 


.96 


14 


1.04 


M 


\m 


Tamarack 


2 


1.00 


.75 


.85 


1 


1.01 


.95 


1.01 











Oak-hickory 


4 


1.00 


.70 


.83 


1 


1.05 


1.24 


1.18 


2 


.95 


.86 


.87 


Lowland hdwd. 


2 


1.00 


.71 


.99 




















N. hardwood 


45 


1.01 


.87 


.95 


54 


1.00 


.88 


.94 


36 


1.01 


.93 


.96 


Aspen 


23 


.98 


.88 


.98 


44 


1.02 


.70 


.97 


17 


1.04 


.94 


1,03 


Paper birch 


7 


.96 


1.18 


1.09 


6 


.96 


.72 


.77 


3 


1.04 


1.16 


1.31 


All types 


122 


1.01 


.74 


.95 


148 


1.01 


.94 


1.00 


120 


1.01 


1.00 


1.04 



for the second growth northern hardwood study 
are shown in table 3. Thirty-seven plots had three 
measurements and 62 had only two measure- 
ments. The three national forests had two mea- 
surements only, so only one ratio is presented for 
them in table 4. 

The Cloquet test shows S of 11 types have 
overestimates of mean d.b.h., 9 of 11 have over- 
estimates of trees per acre and 10 of 11 have 
werestimates of basal area. Although the mean 
iameter overestimates are not serious, they are 
otable for their consistency and regular progres- 
ion; for example, 1 or 2 percent after 5 years, 2 or 3 
ercent after 10 years, and 4, 5, or 6 percent after 
.7 years for white spruce-fir, black spruce, north- 
am white-cedar, tamarack, and paper birch types. 
In fact, all types except red pine, white pine, and 
northern hardwoods have nearly the same pro- 
gression of accuracy decay in predicting mean di- 
ameter. There is close agreement in both the type 
if error (over- or underestimate) and progression 
f accuracy decay when evaluating the system on 
lumber of trees. Exceptions appear to be the white 
-ine, tamarack forest types and northern 
lardwoods. The ratios for predicted and observed 
>asal area follow closely the estimates of mean 
liameter and number of trees. 

The Cloquet results confirm earlier forecasts 



that the projection system would overestimate ob- 
served variables. This forecast was based on per- 
sonal observations of the field remeasurement 
crew on the quantity of large shrubs and small 
trees present on Cloquet plots. Any competing veg- 
etation that has affected growth of the trees on the 
plot but is not a part of the recorded observations 
will result in overestimates. 

The major independent teat for northern 
hardwoods was based on the repeated measure- 
ments of a study (NH-2) in Michigan's Upper Pen- 
ninsula (table 3). The ratios are close to one in 
almost every case, and hence warrant little 
discussion. 

The Chequamegon, Nicolet, and Hiawatha Na- 
tional Forest plot data cover a period of 10 years 
(table 4). These plots differ from the Cloquet and 
NH-2 plots in three ways: the plots are 10-point 
clusters of variable radius plots, the trees are mea- 
sured down to 1.0 inch diameter, and the re- 
measurement occurred 10 years after the initial 
measurement. Most notable of the Chequamegon 
results are the relatively large overestimates of 
mortality in the black spruce, tamarack, oak-hick- 
ory, lowland hardwood, northern hardwood, and 
aspen types. One might discount overestimates in 
the first four types as chance occurrences because 
of the low plot numbers. We feel it is more likely 
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because the mean plot diameters tend to be small, 
especially on the black spruce type, and because 
the number of trees represented by a small tree is 
larger than the number represented by a large tree 
(recall these are variable-radius plots). Thus, if the 
projection should assign a mortality status to a 
very small tree, it may mean that 100 trees per 
acre (its tree factor) are considered to have died. 
This points up one of the major differences in the 
way fixed- and variable-radius plots are projected. 

The Nicolet also shows low estimates of num- 
bers of trees for the aspen and paper birch forest 
types. The Hiawatha shows a 53-percent overesti- 
mate of trees per acre for the white spruce-fir type. 
Most of this difference is attributable to three of 
the eight plots. Because we are underestimating 
mortality on some plots by as much as 242 trees 
per acre in only 10 years, it is apparent that some 
variable other than competition is dominating the 
mortality process. 



Precision analysis 

The precision analysis, based on the Cloquet and 
NH-2 data sources, differs from the accuracy anal- 
ysis in two respects: 

(1) we consider only the ratio of predicted to 
actual mean stand diameter (ratios for num- 
ber of trees and basal area have been com- 
puted but are not included here), and 

(2) we consider the ratio for each individual plot 
in computing the mean and standard devia- 
tion of the ratios over time for for each mea- 
surement following the first. 

Again, the analysis is organized by forest type. 
Figure 3 contains graphs for nine of the Cloquet 
types (lowland hardwood has been omitted), and 
figure 4 shows the mean and standard deviation at 
the two points in time for the northern hardwood 
type from NH-2. 

Our concern with this analysis is with the pat- 
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Figure 4. Pattern of mean ratios and standard 
deviation of ratios of projected to actual mean 
plot diameter for the northern hardwood type. 
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tern created by plotting the mean 1 standard 
deviation over time. Specifically, we look at (1) the 
closeness of the mean ratio to 1.00 at each time, 
and (2) the magnitude of the standard deviation 
and its relation to time since the initial measure- 
ment. Concern number 2 is particularly impor- 
tant. 

Figures 3 and 4 show that all mean ratios are 
within 7 percent of 1 after 17 and 11 years, re- 
spectively. For most types the ratios begin near 1 
at 5 years and increase linearly to the value at 17 
years. Jack pine, white spruce-fir, quaking aspen, 
and paper birch have noticeable concave patterns 
to the mean ratio. Tamarack and northern white- 
cedar show signs of a convex pattern. In all cases 
except the red and white pine types, extrapolation 
backward to the initial measurement would yield 
a ratio near 1.0. 

Comparing the standard deviation magnitude 
shows that northern white-cedar is smallest, with 
black spruce, white spruce-fir, and jack pine also 
low. Quaking aspen has the largest standard devi- 
ation, followed by paper birch, red pine, and white 
pine. 

The trend in standard deviation magnitude over 
me shows three patterns; (1) linear, (2) concave, 
nd (3) convex. We consider jack pine, paper birch, 
nd black spruce to be linear. Red pine, white pine, 
aspen, and white spruce-fir are judged concave. 
Tamarack and northern white-cedar appear to ex- 
hibit more of a convex pattern. With only two 



points, a trend cannot be assigned to the northern 
hardwood type (fig. 4). 

The least desirable standard deviation trend 
type would be convex; linear would be better, and 
the most desirable would be concave. Trend type 
alone, however, is not sufficient to rate the good- 
ness of projections by forest type, because it is 
confounded with standard deviation magnitude 
and the ratio trend itself. For example, we rate 
white pine, red pine, aspen, and white spruce-fir as 
having the best standard deviation trend type 
(concave), but none of these rates low on standard 
deviation magnitude. The northern white-cedar 
type shows a linear trend type and it has the lowest 
standard deviation. Because of these conflicting 
relations (fig. 3) it is not possible to say which type 
is projected the best (as measured by predicted vs. 
actual mean diameter). It is, however, possible to 
identify three candidates for "poorest" rating; 
quaking aspen, paper birch, and tamarack. A more 
definitive analysis could probably be made if num- 
bers of trees were also included. 

Table 5 summarizes the results of the precision 
analysis based on individual plot ratios of pre- 
dicted to actual mean diameter for each forest 
type. 



Table 5. Summary of precision analysis findings 
for projections of mean stand diameter 
(all results are based on Cloquet Forest 
plot data except the northern hardwood 
type, which is based on NH-2 from Up- 
per Michigan) 







Compar- 






Trend In 


ative 


Time 




ratio of 


overall 


trend In 


Forest type 


projected to 


magnitude 


magnitude 




actual mean 


of standard 


of standard 




plot diameter 


deviation 


deviation 


Jack pine 


linear 


small 


linear 


Red pine 


concave 


medium 


concave 


White pine 


concave 


medium 


concave 


White spruce-fir 


concave 


medium 


concave 


Black spruce 


linear 


small 


linear 


N. white-cedar 


convex 


small 


convex 


Tamarack 


convex 


medium 


convex 


Northern hardwood 


NA 1 


medium 


NA 1 


Quaking aspen 


concave 


large 


concave 


Paper birch 


linear 


large 


linear 



1 Not applicable because only two points in time are available. 
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SUMMARY 

Tests of the projection system were baaed on 
comparisons of predicted and actual mean diame- 
ters, number of trees per acre, and basal area per 
acre for two data sets: (1) selected properties from 
the calibration data base, and (2) an independent 
data base. A summary of the important findings, 
organized around two central questions, is given 
below. 

I. Is the method of combining system parts sat- 
isfactory? 

1 . The method of combining system parts 
(potential, modifier, allocation rule, mor- 
tality, etc.) is satisfactory. 

2. Many of the calibration data base proper- 
ties are projected very well for up to 30 
years (fig. 2). 

3. Properties with small trees are generally 
less well projected (e.g., NC-26, NC-32, 
NC-18, Nettleton, NC-59, fig. 2) than 
those with large trees. 

4. Sizable growth rate differences for jack 
and red pine appear between Lower 
Michigan and Wisconsin-Minnesota. If 
validated by further checks, this indi- 
cates the need for separate coefficients for 
at least Lower Michigan. 

5. Projections of average diameter are bet- 
ter than numbers of trees, which are bet- 
ter than basal areas. 

II. Is the projection system applicable to other 
areas in the Lake States? 

1, The projection system did a reasonably 
good job of projecting plots from Minneso- 
ta, Wisconsin, and Upper Michigan. 
These plots had not been used to calibrate 
the system. 

2. Consistent positive bias (as much as 6 
percent in 17 years) was found in pre- 
dictions of mean d.b.h. for several cover 
types on the Cloquet Forest in 
Minnesota. 



Variable-radius plots (10-point cluster 
design of Forest Survey) can be satisfac- 
torily projected for 10 years, although 
mortality is sometimes overestimated for 
plots with many small trees. 
No discernible deficiencies were found in 
projection capabilities for any species in 
the independent data sets as geographi- 
cal location changed. It should be noted 
that we had no independent data set from 
Lower Michigan. 

The standard deviation of the ratio of 
projected to actual mean plot diameter 
showed three types of trends when plot- 
ted over years since the original mea- 
surement: concave, linear, and convex. 
Northern white-cedar and tamarack for- 
est types showed the least desired (con- 
vex) trend. Red pine, white pine, white 
spruce-fir, and quaking aspen showed the 
most desired (concave) trend. 
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The purpose of this report is to summarize 
trends of precipitation and temperature at se- 
lected locations in the Lake States and discuss 
these trends and their possible impact on the cali- 
bration of a forest growth projection study. Sub- 
stantial differences in the rainfall or temperature 
between the period used for calibration and future 
periods could have an effect on the accuracy of the 
projections. 

Previous work on the impact of climate on tree 
growth has indicated a close correlation between 
precipitation and radial stem growth (Dils and 
Day 1952). The transition from wet years to dry 
years induces a gradual, rather than sharp, de- 
crease in tree growth, reflecting a lag in soil mois- 
ture response to decreased rainfall. The reverse 
transition from dry to wet periods produces imme- 
diate growth response because soil moisture is 
replenished more rapidly than it is lost by evapo- 
transpiration. Temperature also affects tree 
growth, primarily because it regulates growth ini- 
tiation in the spring (Dils and Day 1952). Howev- 
er, we have made no attempt to correlate climate 
data and tree growth in this report. 



DATA SOURCES 

Most of the data in this report were obtained 
from records provided by the Minnesota State Cli- 
matology Office, 1903 Hendon Avenue, St. Paul, 
Minnesota 55108. Additional data were obtained 
through correspondence with Forest Service per- 
sonnel in Michigan. 

Precipitation and temperature data for the 6- 
month period, April through September, were 
compiled by month for 11 weather stations 
throughout Minnesota, Michigan, and Wisconsin 
from the year of each station's establishment 



through 1976 (fig. 1). The criterion for weather 
station selection was the station's proximity to 
permanent studies of forest growth. 

In all but two cases, both precipitation and 
temperature were recorded at each station. At the 
Cass Lake, Minnesota station, only precipitation 
was recorded; at the Virginia, Minnesota station, 
temperature was the only variable recorded. 

Precipitation was recorded by month at each 
station from April through September from the 
year of station establishment through 1976. Tem- 
perature was not recorded in the annals monthly 
until 1931; consequently, the values prior to 1931 
are the April through September averages from 
the date of station establishment through 1930. 

Minnesota and Wisconsin temperature values 
were recorded by month from 1931 through 1976. 
Michigan monthly temperature values were not 
recorded until 1941 at Cadillac, and 1951 at Alle- 
gan, Grayling, and East Tawas. No temperature 
information was available for Grayling or East 
Tawas until 1931; from 1931-1951, only the April- 
through-September average was available. April- 
through-September temperature averages were 
recorded for Cadillac and Allegan from station 
establishment to 1941 and 1951, respectively. 



GROWING SEASON 

Prior to discussing the data analysis, the term 
"growing season" as used in this text should be 
defined. Usually the growing season is described 
as the period between the last killing frost in the 
spring and the first killing frost in the fall. Some 
factors that determine whether a plant will be 
killed by frost are species, age, time of year, and 
duration and rate of freezing (Moran and Morgan 
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, , | East Tawai (1897) 
Cadillac (1909) 




Figure 1. Location and date of establishment of meteorological stations. 
Monthly (April to September) precipitation and temperature infor- 
mation was recorded from each station, except Cass Lake, where only 
precipitation was recorded, and Virginia, where only temperature 
was recorded. 



1976). The Weather Bureau commonly designates 
the growing season as the "freeze-free" period, 
when minimum temperatures are greater than 
32F. (0C.). Another measure employs a 28F. 
(~2C.) free period, which accommodates for local 
temperature variations. 

In the Lake States the frost-free period gener- 
ally extends from April through September. 
Moran and Morgan (1976) cite the three Wisconsin 
stations used in this study as having 28F. free 
periods that range from 110 to 195 days. The Min- 
nesota stations, Virginia, Cloquet, and Leech 
Lake, have 28F, free seasons that extend from 104 
to 159 days (Baker and Strub 1963). While most 
tree growth occurs in June and July, it does begin 
as early as April and can continue until Septem- 
ber. Therefore, in this study we have defined the 
growing season as the 180-day period from April 
through September, to include all tree growth oc- 
curring in any year. 



DATA ANALYSIS METHODS 
AND RESULTS 

To prepare the information for analysis, several 
steps were taken. The monthly precipitation data 
were added to obtain growing season rainfall to- 
tals. Average temperatures for each growing sea- 
son were also calculated. 

Five-year averages for both precipitation and 
temperature were calculated for each station from 
its date of establishment through 1976 (figs. 2 and 
3). "Even" 5-year intervals were used for ease of 
manipulation (i.e., 1901-1905, 1906-1910, etc.). 
The first average for a station may include more 
than five yearly measurements to accommodate 
this scheme. For example, Allegan, Michigan pre- 
cipitation data began in 1909, but the first "even" 
interval is 1911-1915; consequently, the first in- 
terval is 7 years, from 1909-1915, All subsequent 
averages are 5 years, except for the last interval, 



91 



15 



MINNESOTA 




O Cass Lake 
D Cloquet 
A Leech Lake 



-A 



WISCONSIN 



20 



16 



25 




O Long Lake 
Q River Falls 
A Waukesha 



16 - 



MICHIGAN 



O Allegan 
P Cadillac 
A Grayling 
EastTawas 

I 



,o 




1000 



1910 



1920 



1830 



1940 
YEAR 



1950 



1960 



1970 



Figure 2. April-to-September precipitation trends in the Lake States, Five year 
averages from selected weather stations were plotted at the endpoint 
of each interval. 
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Figure 3.April-to-September temperature trends in the Lake States. Five year 
averages from selected weather stations were plotted at the endpoint 
of each interval. Dashed lines indicate the average for periods when 
monthly or yearly averages were missing. 



1971-1976. The average for an interval is plotted 
at the endpoint of that interval. 

Trends in precipitation and temperature were 
analyzed for three periods: (1) station establish- 
ment through 1976; (2) station establishment 
through 1935; and (3) 1936-1976. The period 1936- 
1976 encompasses the interval when field mea- 
surements for the forest growth projection study 
were taken, and is referred to as the calibration 
period (Christensen et al. 1979). In addition, the 
months of June and July were analyzed separately 
for the calibration period. Soil moisture is critical 
during these months, because it is then that soil 
moisture deficits are first apparent, and it is also 
the time when most growth occurs (Fraser 1962). 



Yearly fluctuations in precipitation for the vari- 
ous weather stations are large, yet there are some 
notable trends. For example, in Minnesota from 
1911 through 1976, growing season precipitatioi 
at the Leech Lake, Cass Lake, and Cloquet sta 
tions increased slightly overall, as shown by linea 
regression analysis. Until 1940 a decrease in rain- 
fall was recorded for all three stations. Cass Lake 
showed a marked increase in precipitation from 
1941 to 1976 (fig. 2a). Leech Lake showed an in- 
crease in rainfall from 1940 to 1945, but a gradual 
decline from 1946 to 1976. 

Cloquet showed a decrease in precipitation from 
1940 to 1976. From 1936 to 1976 for June-July 
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rainfall gradually increased at all three Minneso- 
ta stations. 

In Wisconsin, precipitation increased from 1918 
to 1976. From station establishment to 1976, the 
Long Lake and River Falls stations showed in- 
creasing rainfall amounts, but the Waukesha sta- 
tion recorded a decrease (fig. 2b). The period from 
station establishment to 1935 was generally drier 
than the period from 1936-1976. The 1936-40 in- 
terval marked an increase in rainfall for a decade. 
Since 1940, all Wisconsin weather stations have 
recorded an overall decrease in precipitation. For 
the June-July analysis, Long Lake and River Falls 
showed a decrease in the amount of rainfall, while 
Waukesha showed an increase. 

In Michigan the overall trend from 1909 to 1976 
was one of decreasing rainfall (fig. 2c). From sta- 
tion establishment through 1976, Allegan and 
Cadillac exhibited this downward trend; however, 
Grayling and East Tawas showed an increase in 
rainfall for the same period. From station estab- 
lishment to 1935, East Tawas experienced pro- 
gressively more rainfall, while the other Michigan 
stations experienced a drying trend. During the 
calibration period these trends continued. Howev- 
er, the June-July analysis showed decreasing pre- 
cipitation at Grayling and increased precipitation 
at the other stations. 

In the Lake States overall, the period from the 
various stations' establishment to 1935 was gener- 
ally drier than the calibration period (1936-1976). 
Following 1935, which marked the end of the 
"Dust Bowl" era, precipitation increased for a dec- 
ade throughout the Lake States. This was followed 
by a general drying trend through the 1976 grow- 
ing season at all stations in the study, except at 
East Tawas and Cass Lake. 

The average precipitation for the Lake States 
from station establishment through 1976 was 19.5 
inches. From station establishment through 1935 
the average was 18.7, and for the calibration peri- 
od, 1936-1976, the average was 20.0 inches. Wis- 
consin had the highest averages: 21.0, 20.0, and 

21.4 inches for the periods from station establish- 
ment-1976, station establishment- 193 5, and 1936- 
1976, respectively, Minnesota had an overall aver- 
age of 18.7 inches, a station establishment- 193 5 
average of 17.4, and a calibration period average of 

19.5 inches. Michigan averaged 18.9 inches over- 
all, 18.6 inches from station establishment- 193 5, 
and 19.1 inches during the calibration period. 



Temperature trends for the study period did. not 
fluctuate as much as precipitation patterns. Sta- 
tion temperature values were fairly constant 
through time, but varied by location. Allegan had 
the highest calibration period average, 63.2F* 
(17.3C.), followed by Waukesha and River Falls 
with 62.0F. (16.6C.) and 61.5F., (16.4C.), re- 
spectively. Cloquet had the lowest average, 
56.3F. (13.4C.). 

Because 5-year averages in Minnesota and Wis- 
consin were not available prior to 1931, and be- 
tween 1941 and 1951 in Michigan, the information 
on temperature trends is limited to the calibration 
period. 

Each State showed a slight cooling trend for the 
calibration period when linear regression analysis 
techniques were applied. On a station-by-station 
basis, however, only six stations exhibited a cool- 
ing trend. The other four stations Leech Lalce, 
River Falls, Waukesha, and Cadillac indicated 
slight warming trends (fig. 3). 

When June-July temperatures were examined., 
they were similar to the growing season trends at 
all stations except Cloquet, which showed gradu- 
ally increasing June-July temperatures. 

The average Lake States growing-season tem- 
perature for the period from station establishment 
to 1976 and for the calibration period was 58.9K. 
(14.9C.). The corresponding averages for these 
two time periods by State are as follows: Wiscon- 
sin 60.1F. (15.6C.), and 60.0F. (15.5C.); Mich- 
igan 59. 5F. (15.3C.J for both periods; and 
Minnesota57.0F. U3.9C.) and 57.1F. 
(13.9C.). 

Although the data indicate downward trends 
for both precipitation and temperature during the 
growing season and an increase in precipitation 
for the months of June and July (table 1), these 
trends may not continue into the future. 



IMPLICATIONS OF FINDINGS 

Precipitation and temperature are external fa.c- 
tors affecting tree growth. They fluctuate both 
yearly and over longer time periods. 

It is logical to assume that short-term fluctua- 
tions of both precipitation and temperature are 
normal and that their overall effect on forest 
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Table I. Comparative growing season precipitation and temperature trends from Station establishment to 
1976 and 1936-1976 (for the period 1936-1976, both the April to September and June and July 
trends are summarized) 

PRECIPITATION TRENDS 





Minnesota 


Wisconsin 


Michigan 


Station establ.- 
1976 
April-Sept. 
1936-1976 


Cass 
Lake 


Cloquet 


Leech 
Lake 


Long 

Lake 


River 

Falls 


Waukesha 


Allegan 


Cadillac 


Grayling 


East 
Tawas 


increase 


increase 


decrease 


Increase 


increase 


decrease 


decrease 


decrease 


increase 


increase 


April-Sept. 

June-July 


increase 
increase 


decrease 

Increase 


decrease 
increase 


decrease 

decrease 


decrease 
decrease 


decrease 
increase 


decrease 
increase 


decrease 
increase 


decrease 
decrease 


increase 
Increase 



TEMPERATURE TRENDS' 



Station establ.- 






















1976 






















April-Sept. 


increase 2 


Increase 


Increase 


decrease 


Increase 


increase 


decrease 


decrease 


increase 


increase 


1936-1976 






















April-Sept. 


decrease 2 


decrease 


increase 


decrease 


increase 


increase 


decrease 


increase 


decrease 


decrease 


June-July 


decrease 1 


Increase 


Increase 


decrease 


Increase 


increase 


decrease 


increase 


decrease 


decrease 



'Michigan temperature trends from 1941-1976 at Cadillac, and from 1951-1976 at other stations 
Virginia recording station. 



growth is minimal. These fluctuations are re- 
flected in seasonal growth variation. During years 
of abundant rainfall and warm temperatures, in- 
dividual trees or stands of trees put on more 
growth. Conversely, sparse rainfall and cold or 
extremely hot temperatures inhibit growth (Diller 
1935). If these periods of highs or lows are pro- 
longed, a marked difference in growth rates would 
be expected to follow. 

Since the projection system is calibrated on 
measurements taken between 1938 and 1976, a 
very favorable period for tree growth in the Lake 
States, its use for projections into periods of un- 
known precipitation and temperature must be 
made with caution. If rainfall and temperature are 
higher in the future than during the calibration 
period, the probable result would be underestima- 
tion of tree growth by the models. If, on the other 
hand, the period following calibration is dry and 
cool or dry and exceptionally warm for an extended 
period of time, the models would probably overes- 
timate actual growth. 
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